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The  purpose  of  this  study  was  to  determine  the  rela¬ 
tionship  between  the  timing  o-f  intelligence  updates  on  the 
location  o-f  strategic  relocatable  targets  (SRT's)  and  the 
probability  o-f  detection.  A  simulation  was  designed  to 
represent  a  bomber  mission  attacking  an  SRT.  Only  notional 
numbers  were  used  in  building  the  model,  but  the  model 
should  still  produce  valid  results  with  real  world  data. 

Most  research  involved  classi-fied  information,  and 
imposed  some  limitations.  Many  generalities  were  used  to 
describe  the  problem  to  avoid  using  classified  information. 
Also  verification  and  validation  checks  were  completed  using 
notional  numbers.  However,  the  final  model  is  still  as 
valuable  as  the  accuracy  of  the  input  allows. 

Strategic  relocatable  targets  are  a  growing  concern  of 
strategic  planners.  New  technology  makes  war  fighting  more 
difficult  and  new  procedures  have  to  keep  up  with  the  trends 
that  are  present.  This  project  is  meant  to  remove  some  of 
the  uncertainty  dealing  with  strategic  relocatable  targets. 

By  applying  this  simulation  to  real  world  situations, 
the  probability  of  mission  success  of  the  strategic  bomber 
can  be  improved.  Results  of  simulation  runs  will  indicate 
the  expected  detection  probability  for  each  of  a  series  of 
update  times.  A  decision  maker  can  then  decide  what  level 
of  mission  success  he  wants  as  approximated  by  the  detection 
probab i 1 i t 1 es  and  apply  his  resources  accordingly. 
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Abstract 

This  project's  primary  goal  was  to  determine  the  rela¬ 
tionship  between  the  timing  of  the  last  intelligence  update 
with  the  probability  o-f  detecting  a  strategic  relocatable 
target  (SRT ) .  A  computer  simulation  of  a  bomber  attacking 
an  SRT  was  built  to  develop  the  rel at i onshi p .  Many  charac¬ 
teristics  o-f  a  bomber  aircra-ft  and  an  SRT  were  investigated 
and  included  in  the  model  to  ensure  a  reasonable  representa¬ 
tion  o-f  the  actual  system.  A-fter  integrating  the  signifi¬ 
cant  -factors  into  a  comprehensive  model,  experimental  runs 
were  made. 

Notional  numbers  were  used  throughout  the  project. 
After  completing  the  verification  and  validation  stages, 
experimental  data  was  run  through  the  model  to  demonstrate 
the  format  of  the  output  and  usefulness  of  the  model.  From 
the  example,  inferences  were  made  about  the  true  relation¬ 
ship  of  the  update  timing  and  detection  probabi 1 1 ti es. 

Experimental  results  indicate  that  a  significant  de¬ 
crease  in  detection  probab i 1 i t i es  occurred  when  the  last 
update  on  the  target's  location  is  received  by  the  attacking 
aircrew  at  a  point  where  the  time  remaining  until  entering 
the  search  area  is  less  than  half  the  dwell  time  of  the  SRT. 


Therefore,  this  model  not  only  provides  the  detection  proba- 


bility  a  decison  maker  can  expect  at  each  update  time  but 
also  the  point  where  significant  changes  in  the  probability 
of  detection  occurs. - ' 


HOW  THE  TIMING  OF  UPDATES  ON  THE  LOCATION  CF  STRATEGIC 


RELOCATABLE  TARGETS  CHANGES  THE 
PROBABILITY  OF  DETECTION 


I .  Introduction 


Background 

Fighting  a  war  involving  strategic  relocatable  targets 
(SRT's)  complicates  the  role  a-f  the  strategic  bomber.  Until 
recently,  attacking  strategic  targets  was  straight  -forward. 
A  bomber,  knowing  the  target's  exact  location,  would  -fly  a 
planned  -flight  path,  identity  the  target,  and  release  the 
appropriate  weapon.  However,  times  have  changed.  Today,  a 
new  class  of  targets  has  become  increasingly  important;  a 
class  of  targets  that  can  change  location  throughout  the 
planning  and  execution  phases  of  a  war.  These  strategic 
relocatable  targets  pose  a  new  threat  to  strategic  bombing 
because  of  the  uncertainty  about  their  location.  To  combat 
this  new  threat,  new  procedures  must  be  developed  to  improve 
a  bomber's  chance  of  finding  and  destroying  an  SRT.  But 
before  old  procedures  can  be  discarded  and  new  procedures 
implemented,  information  regarding  SRT's  must  be  gathered. 

Since  the  ultimate  goal  is  to  ensure  bomber  aircraft 
have  the  best  chance  of  finding  arid  destroying  assigned 
targets,  factors  that  significantly  affect  the  bomber's 
chance  must  be  examined.  Attacking  SRT's  involves  many  fac- 


tors  which  can  be  grouped  into  -four  primary  areas: 


target  , 


terrain,  aircraft,  and  updates  (13). 

T arqet .  Target  movement  is  the  main  difference  between 
strategic  relocatable  targets  and  static  targets.  Until  the 
development  of  SRT's,  bomber  aircrews  only  had  to  rely  on 
the  aircraft  navigation  systems  to  guide  them  to  the  target 
area  since  the  target  location  was  known  with  certainty. 
Now,  locating  the  SRT  has  become  an  added  oroblem.  While  the 
SRT's  location  may  be  known  at  takeoff  time,  the  target  may 
move  before  the  aircraft  arrives.  In  trying  to  predict 
where  the  SRT  will  move,  several  questions  come  to  mind. 
How  fast  does  the  target  travel  from  one  location  to 
another?  How  far  does  the  target  travel  before  another  stop 
is  made?  Once  a  new  location  is  found,  how  long  does  the 
target  stay  in  the  same  location  before  another  move  is 
made7  Answers  to  these  questions  will  help  planners  esti¬ 
mate  a  new  SRT  location  for  the  aircraft  arrival  time. 
However,  while  character  1 st 1 cs  of  the  target  itself  may 
provide  some  answers,  information  about  the  target  environ¬ 
ment  will  supply  other  answers. 

Terr ai n .  Terrain  around  the  SRT  can  affect  the  speed 
of  the  target  as  well  as  limit  the  number  of  possible  relo¬ 
cation  sites.  Terrain  features  such  as  dense  forest  or 
mountainous  terrain  can  male  relocation  efforts  very 
difficult.  Average  speed  and  distance  traveled  each  move  in 
this  type  of  terrain  can  be  considerably  lower  than  for 


terrain  that  allows  good  mobility.  Even  more,  other  terrain 
■features  such  as  lakes  or  mountain  peaks  may  prove  to  be 
in-feasible  locations  for  deployment  of  an  SRT.  But  in  either 
case,  information  about  the  terrain  can  be  very  useful  in 
forecasting  the  movement  of  the  target. 

Predictions  about  future  locations  of  an  SRT  are 
easier  if  the  possible  movement  area  around  the  target's 
present  location  is  reduced.  Reductions  in  the  area's  size 
can  be  made  due  to  slower  SRT  speeds  and  shorter  moves  by 
the  SRT  because  of  certain  terrain  features.  In  these 
cases,  some  segments  of  the  movement  area  can  be  discarded 
from  consideration.  A  smaller  movement  area  means  a  smaller 
search  area  for  the  aircraft.  Then,  because  of  other  air¬ 
craft  constraints,  a  smaller  search  area  can  improve  the 
probability  of  mission  success. 

Aircraft.  During  mission  planning,  tradeoffs  among  the 
aircraft  factors  have  to  be  made  to  ensure  mission  success. 
Aircraft  are  limited  by  allowable  search  time  and  fuel 
constraints  when  searching  for  an  SRT  because  of  other 
mission  requirements  such  as  more  targets  to  attack  or 
recovery  procedures.  The  amount  of  area  covered  during  the 
allotted  search  time  will  depend  mostly  on  the  aircraft's 
speed  and  altitude. 

A  faster  speed  will  allow  the  aircraft  to  cover  more 
area  and  a  higher  altitude  will  provide  a  larger  search 
radius  (distance  the  aircraft  can  see) .  However,  flying 
faster  burns  fuel  at  a  higher  rate  and  therefore  could  force 


the  aircrew  to  end  the  search  early  due  to  low  -fuel. 


A1  so , 


exposure  to  enemy  defenses,  aircraft  equipment  limitations, 
and  the  target  detection  factor  (how  well  the  target  can  be 
seen)  are  other  aircraft  factors  that  need  to  be  consid¬ 
ered.  When  making  tradeoffs  between  the  probability  of 
locating  the  target  versus  the  aircraft  probability  of  sur¬ 
vival,  these  three  factors  contribute  to  the  probability  of 
mission  success.  To  a  certain  extent  the  aircraft  routing, 
speed,  etc.,  can  be  adjusted  as  the  mission  requirements 
change.  But  in  order  to  increase  the  likelihood  of  mission 
success,  accurate  intelligence  information  on  the  SRI  must 
also  be  provided. 

Updates.  Updated  information  on  an  SRT's  location 
prior  to  the  aircraft  arrival  could  reduce  the  mission  to 
one  of  an  attack  on  a  static  target.  If  the  final  update  is 
received  "close"  enough  to  the  aircraft  arrival  time,  the 
target  will  not  have  a  chance  to  move  and  the  aircraft  can 
attack  as  it  would  against  a  static  target.  The  "farther" 
(earlier)  away  from  the  aircraft  arrival  time  that  the  final 
update  is  received,  the  more  the  target  has  an  opportunity 
to  move.  This  situation  reduces  the  likelihood  of  target 
detection  because  more  target  movement  translates  into  a 
larger  area  tor  the  aircraft  to  search.  As  the  search  area 
increases,  the  probability  of  finding  the  target  decreases 
due  to  the  extra  time  required  for  the  search  and  increased 
area  for  the  target  to  hide.  Similarly,  as  the  time  of  the 
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last  update  becomes  closer  to  the  aircraft's  arrival  time 
into  the  target  area,  the  probability  of  target  detection 
decreases.  But  how  "close"  should  the  update  be  to  produce 
a  satisfactory  detection  probability? 

The  best  time  to  pass  updated  information  to  the  air¬ 
crew  is  when  the  aircraft  enters  the  target  area.  However, 
more  flexibility  in  allocating  intelligence  resources  would 
be  available  if  updated  information  can  be  given  to  the 
aircrew  at  an  earlier  time  but  not  significantly  reduce  the 
probability  (chance)  of  finding  the  target.  In  other  words, 
instead  of  trying  to  get  updated  information  to  all  aircraft 
at  their  search  area  arrival  time,  the  information  could  be 
passed  to  the  aircrews  at  any  time  during  a  specific  time 
interval.  If  an  aircrew  receives  their  update  within  the 
specified  time  interval,  decision  makers  would  be  confident 
that  target  detection  probabilities  would  not  be  signifi¬ 
cantly  reduced  from  that  of  a  "closer"  update.  By  computing 
and  displaying  the  relationship  between  the  timing  of  the 
update  and  probability  of  detection,  the  dec i si  on-maker 
would  be  able  to  accurately  assess  the  likelihood  of  finding 
and  destroying  the  SRT  and  allocate  intelligence  resources 
appropr 1  at el y . 

Pr  ob 1 em  Statement 

To  attack  strategic  relocatable  targets,  bomber  crews 
must  receive  intelligence  updates  on  the  location  of  the 
target  in  a  timely  manner.  However,  not  much  is  known  about 


how  the  time  of  the  update  affects  the  probability  of  target 
detection. 

Ob jecti ves  of  the  Study 

Three  primary  objectives  highlight  the  overall  goal  of 
this  research  project.  The  first  objective  is  to  model  the 
movement  of  a  strategic  relocatable  target.  Character i sti cs 
of  SRT's  are  investigated  to  identify  important  factors 
affecting  the  target  movement.  A  second  objective  is  to 
write  a  computer  program  to  simulate  an  aircraft  mission 
against  an  SRT.  One  computer  run  will  simulate  a  series  of 
aircraft  missions  and  a  detection  probability  will  be  com¬ 
puted  for  that  set  of  missions.  More  computer  runs  will  be 
done  for  a  specified  number  of  update  times  with  a  detection 
probability  computed  for  each  time.  Finally,  the  third 
objective  is  to  plat  the  detection  probabilities  against  the 
correspond! ng  update  times.  Confidence  bands  will  also  be 
computed  to  specify  the  range  of  accuracy  of  the  estimate. 
Completion  of  these  objectives  will  accurately  describe  the 
relationship  between  update  times  and  detection  probabil¬ 
ities. 

Li terature  Search 

Class! f i ed  Documents.  Much  is  still  unknown  about 
strategic  relocatable  targets  and  most  information  published 
is  classified.  Headquarters  Strategic  Air  Command  (SAC)  is 
very  interested  in  any  aspect  of  SRT's  for  future  mission 
planning  and  has  written  documents  concerning  SRT's.  Clas- 


sified  documents  outlining  SAC's  views  on  SRT's  were  re¬ 
viewed  for  background  information.  While  classified  infor¬ 
mation  will  not  be  discussed  in  this  paper,  the  documents 
supplied  valuable  information  about  which  variables  concern¬ 
ing  strategic  relocatable  targets  need  to  be  considered  in 
the  simulation  model.  But  just  understanding  the  variables 
involved  does  not  answer  all  the  questions.  Moving  from 
areas  of  certainty  (static  targets)  to  areas  of  uncertainty 
(SRT's)  also  requires  some  background  in  probability  theory. 

Probabi 1 i ty  Theory.  Basic  concepts  about  probability 
distributions  were  reviewed  to  discover  how  an  SRT's  future 
location  could  be  estimated.  Discussions  on  stochastic 
processes  provided  insight  about  how  probability  distribu¬ 
tions  change  from  the  present  time  to  future  times  (4:376). 
This  information  helped  devise  an  initial  unconditional 
probability  distribution  of  the  target's  location.  Proce¬ 
dures  were  then  investigated  to  develop  the  probability 
distribution  of  an  SRT's  movement  and,  together  with  the 
initial  conditions,  to  estimate  the  probability  of  future 
target  locations.  From  this  set  of  possible  locations  a 
search  pattern  can  be  developed  for  the  aircraft  to  ensure 
the  "best"  areas  are  given  first  priority  in  the  search. 

Sear ch  and  Detecti on  Theory.  To  estimate  the  proba¬ 
bility  of  SRT  detection  for  any  update  time,  the  model  must 
include  a  method  for  the  aircraft  to  search  for  the  target. 
Research  uncovered  several  sources  to  help  decide  on  the 
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appropriate  pattern  to  ensure  a  high  probability  of  mission 


success. 


Many  journal  articles  have  been  written  on  search  and 
detection  theory.  In  his  book,  Search  and  Detect i on  Theory , 
Alan  R.  Washburn  summarizes  the  past  -findings  and  claims 
that  up  to  now  search  theory  is  not  really  concerned  with 
optimal  search  "paths'1  which  are  required  for  a  realistic 
search  pattern  for  an  aircraft  (16:1-1).  He  states  two 
reasons  as  the  primary  causes.  First,  an  optimal  path  is 
not  useful  unless  it  is  easy  to  follow.  Second,  he  asserts 
that  a  "path"  is  not  a  convenient  mathematical  object  (16:1- 
1).  Based  on  these  reasons,  most  research,  according  to 
Washburn,  involves  a  distribution  of  effort  over  an  entire 
search  area  to  maximize  the  probability  of  detecting  the 
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target  (16: 1-1 ) . 

A  review  of  recent  literature  concerning  search  and 


detection  theory  confirmed  Washburn's  observation.  S.  S. 
Brown  in  one  of  his  recent  articles  elaborated  on  an  optimal 
search  plan  that  allocated  a  search  effort  in  each  of  a 
finite  number  of  time  intervals  to  maximize  the  overall 
probability  of  detecting  a  moving  target  (3:1275).  A  book 
by  L.  D.  Stone  stated  the  basic  problem  in  search  theory  is 
that  of  maximizing  the  probability  of  detecting  the  target 
within  a  limited  resource  such  as  time  or  fuel  (12:32). 
Finally,  another  article  by  Luke  Tierney  and  Joseph  Kadane 
defines  a  "search  strategy"  as  "an  allocation  of  the  avail¬ 
able  search  effort"  (15:720).  In  these  and  other  examples, 


while  a  search  path  is  not  directly  ignored,  emphasis  is 
placed  on  an  allocation  o-f  search  effort  instead  of  the 
search  pattern  used  to  find  the  target. 

A  myopic  search  algorthim  presented  by  Washburn  can  be 
used  to  develop  a  search  pattern  for  an  aircraft  <16:6-4  - 
6-8).  Patterns  resulting  from  a  myopic  search  may  prove 
infeasible  for  an  aircraft  to  fly  but  by  modifying  the 
procedure  slightly  the  probabilities  derived  will  initially 
put  an  upper  bound  on  the  probability  of  finding  the  target. 

Another  source  of  information  concerning  SRT  '  s  was  a 
segment  of  HO  Air  Force  Studies  and  Analyses  Advanced  Pene¬ 
tration  Model  (APM)  .  The  APM  simulates  an  entire  bomber 
mission  (14:slide  3).  In  the  model,  movement  of  an  SRT  is 
disregarded  after  the  war  begins.  But,  because  some  aircaft 
missions  may  be  quite  long,  SRT  movement  can  almost  be 
assured  before  the  aircraft  arrives.  Also,  another  problem 
with  the  APM  segment  is  that  an  exhaustive  search  (a  pattern 
that  searches  the  entire  area)  is  used.  (14:slide  10).  An 
exhaustive  search  may  prove  infeasible  because  required 
search  time  is  too  long  and  aircraft  fuel  reserves  could 
become  too  low  to  fulfill  other  mission  requirements  or  to 
recover  safely.  Strategic  planners  today  envision  only  one 
trip  through  the  search  area  with  at  the  most  only  one  turn 
not  to  exceed  90  degrees  (Fig.  1.1)  (5,13). 
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Fig.  1.1  Expected  Flight  Path  of  Aircraft  Through 

a  Search  Area 


Scope  and  Limitations 

Target  T ype.  Knowledge  of  the  type  of  SRT  being  at¬ 
tacked  is  required  in  order  to  model  the  movement  of  the 
target  and  predict  future  locations.  SRT  '  s  are  either 
time  constrained  or  time  critical  targets  (11:2).  Time 
constrained  targets  have  a  dwell  time  (time  remaining  in  one 
spot)  greater  than  the  weapon  system  cycle  time  (11:2). 
Weapon  system  cycle  time  is  the  total  time  required  from 
detection  of  a  change  in  target  location  to  weapon  system 
arrival  (10:7).  Time  critical  targets  have  a  dwell  time 
less  than  the  weapon  system  cycle  time  (11:2).  Classifying 
an  SRT  in  one  of  these  two  categories  allows  an  appropriate 
weapon  to  be  selected  for  the  attack. 

Weapon  T ype .  One  of  the  problems  facing  strategic 
planners  is  the  type  of  weapon  system  to  use  against  an  SRT. 
In  general,  either  missiles  or  bombers  can  be  used  against 


SRT ’ s  (10:7).  Missiles  can  get  to  the  target  faster  but 


lack  the  -flexibility  of  a  bomber. 


missiles  need  an 


A1  so , 

accurate  target  location  when  they  are  launched.  Therefore, 
missiles  can  be  used  more  effectively  against  time  con¬ 
strained  SRT  '  s  since  no  movement  will  take  place  before  the 
missile  arrives. 

On  the  other  hand,  bombers  are  needed  for  time  critical 
SRT ' s  because  movement  of  the  target  will  occur  even  before 
a  missile  could  arrive.  Bombers  provide  the  flexibility  to 
apply  a  constrained  search  far  the  SRT  and  attack  when  the 
target  is  found  (11:61).  Since  bombers  are  the  primary 
concern  of  this  project,  the  assumed  scenario  will  be  bom¬ 
bers  attacking  time  critical  targets  or  targets  with  very 
short  dwell  times.  But  before  a  bomber  can  attack  an  SRT, 
information  concerning  the  target  itself  is  needed. 

SRT  Characteristics.  Research  into  the  SRT  question 
shows  that  some  factor  values  about  SRT ' s  are  known  with 
certainty  but  little  is  known  about  other  factors.  First, 
even  though  SRT ' s  move  around,  they  do  need  a  home  base  for 
certain  support  requirements  (11:2).  This  will  limit  the 
full  range  of  movement.  Also,  unless  a  war  is  in  progress, 
there  is  no  reason  to  relocate  the  target  (11:2).  Movements 
from  the  home  base  may  occur  due  to  training  but  these 
distances  will  be  small.  Therefore,  it  can  be  assumed  that 
at  the  start  of  hostilities  the  target  will  be  at  its  home 
base  and  will  not  stray  more  than  a  specified  distance  from 
that  base  throughout  the  war.  Further,  target  characteris¬ 


tics  such  as  speed,  dwell  time, 


target  type  (mobile  or 


movable),  set  up  and  tear  down  times,  and  possible  location 
sites  are  all  assumed  to  be  known. 

Methodol oqy 

Si  mul  at  i  on  Model  .  Because  of  the  complexities  and 
number  of  variables  involved  in  computing  detection  proba¬ 
bilities  of  SRT's,  simulation  was  chosen  as  the  general 
approach  to  solve  the  problem.  Simulation  is  a  way  to 
represent  the  behavior  of  a  system  as  it  changes  from 
state  to  state  according  to  a  set  of  well-defined  operating 
rules  (9:6).  For  the  purposes  of  this  project,  the  state  of 
the  system  will  be  the  location  of  the  target  at  each  time 
period.  Each  type  of  SRT  in  every  possible  region  repre¬ 
sents  a  new  situation  requiring  different  values  for  the 
variables  of  the  system  (dwell  time,  setup  time,  etc.).  For 
example,  an  army  garrison  moving  in  mountainous  terrain 
would  require  different  input  values  than  a  motile  missile 
unit  in  flat  terrain.  Many  situations  such  as  these  can  be* 
modeled  with  a  computer  simulation  by  inputinq  different 
combinations  of  values  for  the  factors. 

The  flexibility  of  simulation  allows  the  desired  re¬ 
sults  to  be  computed  with  an  opportunity  to  test  other 
aspects  of  the  SF;T  mission.  Tor  instance,  the  army  garrison 
mentioned  ear  1 1 er  would  have  different  input  values  for 
speed,  detection  factors,  set  up  and  tear  down  times,  etc. , 
than  the  mobile  missile  unit.  However,  both  situations  can 
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be  modeled  by  just  changing  the  input  information  with  no 
changes  required  in  the  actual  model. 

Developing  the  simulation  model  involves  two  different 
theories.  First,  target  movement  is  modeled  as  a  stochastic 
process.  A  stochastic  process  is  a  collection  of  numbers 
each  of  which  represents  the  value  of  a  particular  variable 
for  each  specified  period  of  time  (2:1).  Second,  optimal 
search  theory  and  associated  search  patterns  are  investi¬ 
gated.  One  search  pattern  will  be  used  for  all  aircraft 
missions  to  ensure  results  are  not  biased  toward  a  particu¬ 
lar  search  pattern. 

Desi gn  of  Exper i ments.  A  statistical  design  of  experi¬ 
ments  is  the  "process  of  planning  an  experiment  so  that 
appropriate  data  will  be  collected  which  may  be  analyzed  by 
statistical  methods  resulting  in  valid  and  objective  conclu¬ 
sions"  (7:2).  When  the  experiment  is  completed,  significant 
factors  (those  that  have  significant  effects  on  the  results) 
will  be  identified.  Knowledge  of  these  factors  will  aid  in 
the  future  analysis  of  the  computer  runs. 

A  scientific  approach  was  used  to  determine  which  fac¬ 
tors  are  critical  to  the  final  results  of  the  simulation. 
To  determine  the  relationship  between  update  times  and  de¬ 
tection  probabilities,  all  critical  factors  must  be  held 
constant  as  the  update  time  is  changed  to  ensure  the  change 
in  the  resulting  detection  probability  is  only  due  to  the 
change  in  the  update  time  and  not  other  critical  factors. 
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Then  a  good  estimate  -for  the  detection  probability  at  each 
update  time  can  be  calculated. 

Estimating  the  Pr obab i 1 i ty  of  Detecti on.  An  aircraft 
mission  will  be  simulated  as  i t  attacks  an  SRT.  Movement  by 
tne  target  will  be  started  at  time  zero  of  the  simulation 
and  continue  until  one  of  the  termination  conditions  is  met. 
An  update  on  the  target's  location  will  be  simulated  at  some 
specified  time  be-fore  the  aircraft  begins  the  search.  Apply¬ 
ing  the  updated  information,  the  aircraft  will  search  for 
the  target  and  record  either  a  hit  or  a  miss  for  each  mis¬ 
sion  simulated.  Many  missions  will  be  flown  at  each  update 
time  and  the  percentage  of  hits  will  be  used  as  an  estimate 
of  the  detection  probability  for  th^t  update  time.  Each 
probability  with  the  corresponding  update  time  can  then  be 
plotted  to  display  the  desired  relationship. 

Sequence  of  Presentat i on 

Following  Chapter  1,  Chapter  2  explains  the  theory  and 
procedures  used  to  model  the  target  movement.  Chapter  3 
continues  by  describing  the  aircraft  mission,  how  the  proba¬ 
bility  distribution  for  the  target's  location  is  computed, 
and  how  the  aircraft  searches  for  the  target.  In  Chapter  4, 
the  simulation  model  is  explained  as  well  as  how  the  air¬ 
craft  and  target  models  are  integrated.  Chapter  5  dis¬ 
cusses  the  validation  and  verification  of  the  model.  An 
experimental  design  used  to  identify  key  relationships  among 
the  factors  is  explained  in  Chapter  6  along  with  the  analy- 
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I  I .  Model i ng  T arget  Movement 

Introducti on 

An  aircrew  must  have  some  information  about  an  SRT ' s 
pattern  of  movement  to  make  accurate  estimates  about  the 
target's  future  locations.  Chapter  2  begins  by  relating  an 
SRT  '  s  movement  to  a  stochastic  process  and  explains  how  the 
properties  of  a  Markov  chain  are  used  to  model  the  target ' s 
movement.  Next  is  a  discussion  on  the  meaning  and  develop¬ 
ment  of  the  one-step  transition  matrix.  Finally,  a  descrip¬ 
tion  of  the  two  types  of  SRTs  is  presented  as  well  as  how  to 
compute  the  time  unit  of  the  stochasic  process  for  each  type 
of  target. 


Stochast i c  Processes 

An  aircrew  of  a  strategic  bomber  preparing  to  attack  an 
SRT  is  interested  in  the  location  of  the  target  at  the  time 
the  aircraft  enters  the  search  area.  At  some  point  in  time 
prior  to  an  attack,  on  an  SRT,  intelligence  information  would 
report  the  exact  location  of  the  target  as  is  currently  done 
with  static  targets.  However,  since  an  SRT  moves,  the 
location  at  (future)  times  t,,  t  =  ,  t3,...  is  desired.  At 
each  time,  tA ,  a  value  is  observed  for  the  random  variable 
shown  as  Xt i ,  , . . . ( 1 : 293 ) .  A  set  of  random  varia¬ 
bles  CX*.  i  :  ti  6  T>  is  called  a  stochastic  process  and  T  is 
called  the  index  set  Cti,  t3 ,  ts,...l  of  the  stochastic 


process  ( 1 : 293) . 


To  obtain  analytical  results  when  evalua- 


ting  the  process,  assumptions  about  the  joint  distribution 
of  Xfci,  Xta,  are  necessary  (4:372). 

One  passible  assumption  is  that  the  stochastic  process 
is  a  Markov  chain.  A  Markov  chain  assumption  eases  the 
analytical  computations  in  making  inferences  about  a  sto¬ 
chastic  process  (4:372).  However,  a  -finite-state  Markov 
chain  must  satisfy  four  requirements  (4:373): 

1.  A  -finite  number  of  states 

2.  The  Markovian  property 

3.  Stationary  transition  pr obab i 1 i t i es 

4.  A  set  of  initial  probabilities  F'tX0=il  for  all  i 
To  explain  how  each  requirement  is  applied  to  an  SRT  and  to 
illustrate  how  SRT  ’  s  can  be  modeled  as  finite-state  Markov 
chain,  the  following  discussion  introduces  a  small  example. 

Finite  Number  of  States.  In  describing  an  SRT ' s  move¬ 
ment  as  a  stochastic  process,  the  state  of  the  system  at 
time  ti  is  the  location  of  the  target.  To  define  the  loca¬ 
tion,  each  SRT  is  assumed  to  stay  within  a  specified  range 
of  its  main  operating  base  (MOB)  due  to  support  requirements 
(11:2).  Therefore,  a  square  box  can  be  drawn  around  the  MOB 
to  define  the  area  of  movement.  Within  the  large  square, 
smaller  square  cells,  which  will  be  referred  to  as  grids, 
can  be  drawn.  Even  though  the  size  of  the  grids  may  vary,  a 
finite  or  countable  number  of  grids  will  always  exist.  The 


state  of  the  system  (target  location)  at  time  tt  can  be 
defined  as  the  number  of  the  grid  occupied  by  the  target  at 


t i me  t i 


As  an  example,  suppose  an  SRT  is  restricted  to  a  15 
mile  radius  of  its  HOB.  Assume  -further  that  the  grid  size 
is  10  miles  square.  A  square  drawn  around  the  possible 
movement  area  and  divided  into  appropriate  grids  would  re¬ 
flect  the  situation  depicted  in  Fig.  2.1.  Nine  grids 
(states)  make  up  the  possible  location  for  the  SRT.  The 
grids  are  numbered  from  left  to  right  and  down  the  box  as 
shown.  Whereas  maximum  travel  distance  by  the  SF:T  and  grid 
size  may  change,  each  situation  involving  an  SRT  can  be 
modeled  in  this  way.  A  finite  number  of  states  is  therefore 


established  for  the  process. 


Fig. 2.1.  Typical  Search  Area 
Markov! an  Property .  Simply  stated,  any  stochastic 
process  with  the  Markovian  property  implies  the  future  state 
of  the  system  depends  only  on  the  current  state  and  not  on 
any  state  observed  in  the  past  (1:294).  In  a  more  rigorous 
fashion,  the  conditional  probability  of  any  future  state 
given  any  past  state  is  independent  of  the  past  state  and 
depends  only  on  the  present  state  (4:372). 

Assume  the  SRT  in  the  previous  example  starts  at  its 


MOB  in  grid  5  and  moves  to  qrid  4,  1,  and  2  in  three  consec- 


utive  time  periods.  The  conditional  probability  of  moving 
■from  grid  2  to  any  other  grid  is  independent  of  how  the  SRT 
got  to  grid  2.  Being  in  grid  2  at  the  present  time  is  the 
driving  factor.  If  the  SRT  had  moved  to  grids  6,  3,  and  2, 
the  conditional  probabilities  for  the  next  move  would  be  the 
same  as  if  the  SRT  moved  through  4,  1,  and  2.  Not  only  does 
this  fact  show  the  existence  of  the  Markovian  property,  but 
it  introduces  the  idea  of  conditional  probabilities. 

St  at  ionary  Transition  F'robabi  I  i  ti  es.  Conditional  prob¬ 
abilities  in  a  Markov  chain  are  called  transition  probabili¬ 
ties  (4:372).  Mathemat i cal  1 y ,  transition  probabi 1 i ties  can 
be  defined  as  P  CXt-.i  =  j|  Xt=i  3  which  represents  the  probabil¬ 
ity  that  the  state  of  the  process  at  time  t+1  will  be  j 
given  at  time  t  the  state  of  the  process  is  i  (4:372).  To 
be  considered  stationary,  an  added  restriction  is  intro¬ 
duced  . 

To  be  stationary,  the  probability  of  moving  from  state 
l  in  time  period  t  to  state  j  in  time  period  t+1  must  be 
equal  to  the  probability  of  moving  from  state  i  in  time 
period  0  to  state  j  in  time  period  1.  This  relationship 
must  hold  for  all  possible  grids  i  and  j.  In  mathematical 
notat l on : 

f  >-X  t.  i  =  j|  X  t  =  i  )  —  F  x  X  i  =  j  |  X  o—  l  3  for  all  t  — 0,1,*..,  (1) 
The  individual  probabi 1 i t l es  are  denoted  by  pt j  (4:372). 
Another  look  at  the  previous  example  will  help  define  the 
transition  pr obab i 1 i t i es . 

Starting  from  grid  2,  assume  the  target  has  a  0.3 


probability  of  moving  to  grid  1,  a  0.4  probability  of  moving 
to  grid  4,  and  a  0.3  probability  of  moving  to  grid  5.  The 
situation  is  depicted  in  Fig.  2.2.  Assume  further  that 
three  time  periods  have  passed.  Identifying  one  case.  Fig. 
2.2  shows  the  the  F'  CX^=1  |  X3=21  =  0.3  where  t=3  is 


Fig.  2.2.  Probable  Movement  of  Target  From  One  Grid 


the  current  time  period.  An  equivalent  statement  could 
be  made  about  grid  4  or  grid  5.  To  be  considered  station¬ 
ary,  PCX^=1  |  X 3=2 >  =  P  C X i = 1 1  X o=  2 >  =  0.3.  In  other  words, 
once  in  grid  2,  the  conditional  probability  of  moving  to  any 
other  grid  must  be  the  same  whether  the  move  takes  place  at 
time  0,  at  time  1,  or  at  time  n.  There  is  no  evidence  in  the 
available  literature  (10,11)  to  suggest  that  restrictions 
are  placed  on  SRT's  that  would  change  the  conditional  proba¬ 
bilities  from  one  time  period  to  another;  therefore,  the 


probabilities  are  assumed  to  remain  the  same  for  all  time 
periods  and  can  be  considered  stationary.  Conditional  prob¬ 
abilities  give  insight  into  the  probab i 1 i t l es  of  transition¬ 
ing  from  one  state  to  another,  but  if  the  unconditional 
probability  of  a  future  location  of  the  SRT  is  desired,  the 


probability  distribution  of  the  initial  state  is  also  re¬ 


quired  (4: 376) . 

Initial  Condi t i ons .  In  order  to  specify  the  initial 
probability  distribution,  the  SRT's  location  must  be  known 
at  some  point  in  time.  Considering  the  example  again,  if 
the  SRT  is  observed  in  grid  2,  the  time  of  observation  is  t0 
and  X0=2.  The  initial  probability  distribution  would  show 
the  P£X0=23  =  1  and  PCX0=i3  =  0  for  all  i*2.  Knowledge  of 
this  distribution  satisfies  the  fourth  requirement  for  a 
Markov  chain  and  combined  with  the  conditional  (transition) 
probabi 1 i t i es ,  future  uncondi t i onal  probabilities  can  be 
calculated.  A  convenient  method  for  computing  the  probabil¬ 
ities  is  through  matrix  operations. 

What  i s  a  Transi tion  Matr i x 

Transition  probabilities  can  be  represented  by  matrix 
notation  as  follows  (4:373): 


where  p± j  =  probability  of  going  from  state  i  to  state  j  in 
one  step  (one  time  period).  Using  a  3  X  3  grid,  assume  the 


one  step  transition  matrix  is  as  follows: 


starting  point  in  the  example  is  grid  2.  Then  the  initial 
probability  distribution  is  p0  =  CO  1  0  O  0  0  0  0  03. 
Multiply  the  row  vector ,  p©,  by  the  one  step  transition 
matrix,  P<4>,  and  the  results  are  pt  =  C.3  0  0  .4  .3  0  0  0 
03  which  means  starting  in  grid  2,  after  one  time  period 
there  is  a  0.3  probability  of  being  in  grid  1,  0.4  probabil¬ 
ity  for  grid  4,  and  0.3  probability  for  grid  5.  If  p0  were 
multiplied  by  the  2  step  transition  matrix,  P<a>,  the  uncon¬ 
ditional  distribution  would  represent  probab i 1 l t i es  after 
two  time  periods. 

In  any  case,  the  i nter pr etat i on  is  the  same.  Multi¬ 
plying  the  initial  probability  distribution  vector  times  the 
desired  n-step  transition  matrix  results  in  the  uncondi¬ 
tional  probability  distribution  for  the  state  of  the  process 
after  time  n.  Because  so  much  information  is  embodied  in 
the  one  step  transition  matrix,  accurately  developing  the 
conditional  probabilities  is  critical. 

Develop l ng  the  Trans l t i on  Matr l x 

Associating  each  SRT  with  its  area  of  movement  creates 
a  unique  situation.  To  accurately  model  the  target  move¬ 
ment,  an  accurate  one  step  transition  matrix  must  be  devel¬ 
oped.  There  will  be  a  row  and  a  column  for  each  grid  in  the 
movement  area.  From  each  grid,  estimates  must  be  made  about 
the  probability  of  the  SRT  moving  into  any  other  available 
grid  in  the  next  time  period.  Estimating  the  probabi 1 i t i es 
involves  several  considerations. 


Considering  the  new  grid,  three  situations  could  exist. 
First,  the  grid  could  contain  terrain  -features  such  that  an 
SRT  could  not  move  through  the  grid.  In  this  case,  the 
conditional  probability  of  moving  there  would  be  zero.  A 
second  situation  could  be  one  where  the  terrain  is  such  that 
the  SRT  may  not  have  any  place  to  set  up  and  dwell  but  could 
move  through  the  area  enroute  to  another  grid.  Finally,  a 
third  case  is  a  grid  that  will  allow  movement  through  the 
area  and  allow  the  SRT  to  set  up.  Knowing  the  characteris¬ 
tics  of  the  individual  SRT  and  terrain  around  its  operating 
area  will  allow  each  grid  to  be  classified  as  one  of  the 
three  types. 

Another  consi derat i on  in  developing  the  conditional 
probabilities  involves  the  speed  of  the  target.  Faster 
moving  targets  can  reach  grids  farther  away  in  one  time 
period  than  slower  moving  targets.  Therefore,  this  consid¬ 
eration  is  dependent  on  the  time  period  chosen  for  the 
stochastic  process.  However,  grids  farther  from  the  current 
location  will  also  have  a  smaller  probability  of  being 
entered  than  closer  grids. 

Finally,  the  probability  of  staying  in  the  same  grid 
for  two  consecutive  time  periods  will  be  assumed  to  be  zero. 
From  the  known  character  1 st i cs  of  the  SRT,  the  length  of 
time  spent  in  one  grid  can  be  calculated.  These  are  average 
times.  If  there  is  any  chance  of  an  SRT  staying  in  one  grid 
for  more  than  one  time  period,  it  will  be  reflected  in  the 
average  "dwell"  time  of  the  SRT.  In  this  way  a  time  unit 


■for  the  stochastic  process  (time  between  changes  o-f  state) 


can  be  developed  -from  the  average  dwell  time  and  a  change  o-f 
state  must  occur  after  each  time  unit.  Development  of  the 
transition  matrix  is  a  key  to  modeling  SRT  movement  but  the 
type  of  target  being  modeled  must  also  be  known  to  develop 
the  time  unit  of  the  stochastic  process. 


Matchi ng  Target  T ype  and  T i me  Unit 

Strategic  relocatable  targets  can  be  classified  as 
either  mobile  or  movable  targets  (11:2).  A  mobile  target  is 
a  target  that  is  always  in  motion  and  assumed  to  be  at  a 
constant  speed.  Thus  the  amount  of  time  spent  in  one  grid 
can  be  calculated  as  the  grid  size  divided  by  the  target's 
speed  and  the  result  used  as  the  time  unit  of  the  stochastic 
process.  Also,  since  the  mobile  target  does  not  have  to  set 
up,  there  are  only  two  kinds  of  grids  rather  than  three 
(grids  the  SRT  can  enter  and  grids  they  can  not).  A  mobile 
target  will  move  randomly  to  an  available  adjacent  grid 
after  each  time  unit  and  the  state  of  the  system  will  change 
accordingly.  Movable  targets  require  a  little  more  informa¬ 
tion  to  compute  the  time  unit  of  the  process. 

Movable  SRT's  are  targets  that  set  up  in  one  location 
for  a  specified  time  period  then  relocate  to  another  posi¬ 
tion  (11:2).  Two  pieces  of  information  about  the  SRT  must 
be  known  to  accurately  compute  a  time  unit  for  the  stochas¬ 
tic  process.  First,  the  total  time  spent  in  one  spot  must 
be  given.  A  combination  of  set  up  time,  dwell  time,  and 


tear  down  time  (time  required  to  dismantle  equipment  to 
prepare  -for  another  move)  make  up  the  total  dwell  time  spent 
in  one  spot.  Once  the  SRT  is  ready  to  move,  another  time 
variable  is  introduced. 

To  move  -from  one  grid  to  another  ,  the  time  required 
depends  not  only  on  the  distance  to  move  but  also  the  speed 
of  the  SRT.  Distances  between  possible  set  up  points  will 
vary.  Therefore,  an  average  distance  per  move  must  be 
calculated  for  each  area.  For  a  mobile  target,  the  distance 
is  always  one  grid,  but  for  a  movable  target  the  distance 
can  be  longer  (not  shorter  since  it  is  assumed  the  grid 
number  must  change  after  each  time  unit).  Dividing  the 
average  distance  moved  by  the  speed  of  the  target  will 
result  in  the  average  time  required  to  move  from  one  loca¬ 
tion  to  another.  The  actual  distances  moved  will  depend  on 
the  grid  size  selected.  How  the  grid  size  is  determined 
will  be  discussed  in  the  aircraft  model.  For  now,  adding 
the  average  move  time  to  the  total  dwell  time  becomes  the 
time  unit  for  the  stochastic  process  for  a  movable  target. 

Cone  1 usi on 

A  model  for  moving  the  target  can  now  be  built  for  the 
simulation  using  Markov  chain  procedures.  Development  of 
the  one  step  transition  matrix  and  the  time  unit  of  the 
stochastic  process  are  two  key  elements  in  modeling  the  SRT 
movement.  Character i st i cs  of  SRT  '  s  are  available  in  classi¬ 
fied  documents  and  to  get  realistic  results,  this  classified 


information  must  be  accessed  and  used  to  develop  the  key 
input  data.  However,  verification  and  validation  of  the 
model  can  be  done  with  notional  numbers.  But  to  complete  the 
model,  an  aircraft  mission  must  be  included  to  simulate  the 
attack  on  an  SRT  so  that  the  detection  probab 1 1 i t 1 es  can  be 


computed . 


III.  Model i ng  the  Aircraft  Mi ssi on 


Introducti on 


To  estimate  the  probability  of  detection,  an  aircraft 
mission  attacking  an  SRT  must  be  simulated.  At  many  times 
planners  will  have  no  control  over  the  target  variables.  But 
in  the  aircraft  model  most  variables  can  be  set  as  the 
planners  desire.  Therefore,  a  little  more  flexibility  is 
needed  to  model  the  aircraft  mission.  Planners  can  adjust 
aircraft  parameters  such  as  speed  and  altitude  to  improve 
the  probability  of  mission  success.  Since  each  SRT  mission 
is  unique,  simulating  the  aircraft  attack  must  be  flexible 
enough  to  accept  many  values  for  the  aircraft  factors  and 
integrate  them  to  arrive  at  an  estimate  for  the  probability 
of  detection. 

Chapter  3  describes  the  model  of  the  aircraft  mission. 
First  . s  a  discussion  about  what  makes  up  the  system  cycle 
time  and  how  the  aircraft  mission  is  modeled  until  the 
search  begins.  An  explanation  of  the  search  pattern  used  by 
the  aircraft  to  look  for  the  SRT  follows  including  some 
comparisons  among  three  possible  search  patterns.  Third, 
intelligence  updates  are  introduced  with  a  discussion  on  how 
the  aircraft  makes  adjustments.  Finally,  the  fourth  section 
explains  the  conditions  for  terminating  a  single  aircraft 
mi ssi on . 
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System  Cycle  Time 

System  cycle  time  is  the  total  time  required  from 
detection  of  a  change  in  a  relocatable  target's  location  or 
status  until  the  weapon  system  arrives  at  the  target  (10:7). 
Because  the  target  location  is  not  known  with  certainty,  the 
weapon  system  arrival  time  assumed  in  this  model  will  be  the 
arrival  time  of  the  aircraft  at  the  edge  of  the  search  area. 
System  cycle  time  can  be  broken  into  five  segments  but  with 
some  overlap  involved  in  certain  situations. 

Intel  1 i gence  Cyc 1 e  T i me.  Intelligence  cycle  time  is 
the  first  part  of  the  total  time.  Time  starts  when  a  change 
in  an  SRT's  location  is  detected  and  ends  when  the  informa¬ 
tion  is  in  the  hands  of  the  mission  planners  (13).  Usually 
military  staff  planners  or  the  crew  members  flying  the 
mission  do  the  flight  and  target  planning.  However,  the 
main  point  is  that  there  is  always  a  delay  from  when  a 
change  occurs  in  an  SRT  location  and  when  plans  can  be  made 
to  attack  the  target  at  the  new  location.  This  delay  is 
input  into  the  aircraft  model  as  a  constant  to  allow  realis¬ 
tic  target  movement  to  occur  while  the  aircraft  is  preparing 
for  the  attack. 

hi  ssi  on  F’l  ann  l  nq  .  The  second  part  of  the  system  cycle 
time  is  mission  planning  time.  Once  new  information  reaches 
planners  another  delay  is  encountered  while  the  actual  air¬ 
craft  mission  is  planned  (13).  Crew  members  or  higher  level 
planners  can  accomplish  this  task.  Mission  planning  time 
may  not  be  applicable  in  some  cases  such  as  when  staff 


officers  plan  the  attack.  In  these  cases,  mission  planning 
time  runs  concurrently  with  crew  rest. 

Crew  Rest .  Crew  members  must  be  afforded  certain  rest 
time  before  -flying  the  required  mission  which  constitutes 
the  third  part  of  system  cycle  time.  As  stated  earlier,  if 
crew  members  do  not  plan  the  mission,  crew  rest  can  run 
concurrently  with  mission  planning  time  and  only  one  delay 
needs  to  be  considered.  However,  even  if  the  time  runs 
concurrently,  a  small  delay  is  still  needed  to  account  for 
briefing  the  crew  on  the  mission  before  they  takeoff. 

Fr ef 1 i ght  T i me .  As  the  fourth  part  of  system  cycle 
time,  any  ground  delays  such  as  preflight  and  taxi  time 
must  be  considered  in  the  model.  Other  delays  may  also  be 
encountered  due  to  alert  time  or  unexpected  maintenance 
problems  prior  to  takeoff  (13).  An  average  expected  delay 
can  be  computed  for  these  circumstances  and  should  always  be 
included  in  the  model . 

FI  1 gh t  T i me.  The  final  delay  in  the  system  cycle  time 
is  the  actual  flight  time  to  the  search  area.  Each  aircraft 
mission  must  account  for  the  time  from  takeoff  until  the 
aircraft  is  in  a  position  to  begin  the  search.  Target  move¬ 
ment  must  be  allowed  to  continue  according  to  its  schedule 
until  the  aircraft  arrival  at  the  search  area. 

After  all  five  delays  are  computed,  the  total  system 
cycle  time  is  calculated  as  the  sum  of  the  five  delays.  In 
order  to  accurately  estimate  the  aircrew  s  probability  of 


detecting  the  target, 


the  target  movement  and  the  aircraft 


mission  must  be  real i sti cal  1 y  integrated.  All  delays  be¬ 
tween  the  time  the  target  location  is  known  and  the  aircraft 
arrival  time  must  be  included  to  allow  the  target  to  move 
far  the  correct  amount  of  time.  If  more  or  less  time  is 
allowed  before  the  aircraft  begins  the  search,  the  detection 
probability  may  be  affected  which  will  then  bias  the  results 
away  from  the  true  probability  of  detection.  For  instance, 
if  not  enough  movement  is  allowed,  the  location  distribution 
for  the  target  will  cover  a  smaller  area  than  it  should 
which  will  incorrectly  reflect  a  higher  probability  of  de¬ 
tection.  However,  if  the  proper  delays  are  accounted  for, 
the  unconditional  probability  distribution  of  the  target's 
location  will  be  a  better  representati on  of  the  real  world 
and  the  search  can  begin  with  information  as  good  as  can  be 
expected. 

Ai r cr af t  Search 

One  search  pattern  must  be  selected  and  used  for  all 
computer  runs  so  that  any  changes  in  detection  probabilities 
can  only  be  attributed  to  changes  in  the  time  of  the  last 
update.  To  accurately  plot  the  relationship  between  detec¬ 
tion  probab i 1 i t i es  and  the  time  of  the  last  intelligence 
update,  all  the  aircraft  factors  except  the  update  time  must 
remain  constant  for  all  computer  runs.  In  so  doing,  a  change 
in  detection  pr obab i 1 i t 1 es  from  one  update  time  to  another 


will  only  be  due  to  the  change  in  the  update  time. 


Most  of 


the  aircraft  factors  have  values  that  can  be  derived  from 
technical  manuals  or  operational  procedures.  However,  the 
search  pattern  is  not  so  easily  defined.  Three  search 
patterns  were  considered  as  possible  ways  to  look  for  SRT's. 

Exhaustive  Search.  When  using  an  exhaustive  search 
pattern,  the  entire  search  area  is  scanned.  As  depicted-in 
Fig.  3.1,  an  aircraft  enters  from  the  upper  left  corner  of 
the  region  and  searches  from  left  to  right  (14:slide  8). 
When  at  the  right  edge  of  the  area  a  180  degree  turn  is  made 


Fig.  3.1.  Exhaustive  Search 


to  reverse  the  direction  and  to  search  the  next  row  of  grids 
in  the  opposite  direction.  At  the  left  edge  the  procedure 
is  repeated.  When  the  last  row  (bottom  row)  has  been 
searched  the  aircraft  exits  the  search  area.  By  defining 
the  size  of  each  square  grid  as  two  tunes  the  search  radius 
of  the  aircraft  sensors,  the  entire  area  is  searched  which 
gives  a  high  probability  of  target  detection;  however,  there 
are  some  drawbacks. 

To  use  the  exhaustive  search  requires  a  lot  of  time  and 


does  not  use  the  updated  in-formation.  Each  aircraft  will 
fly  for  a  considerable  time  to  get  to  the  search  area  due  to 
the  location  of  probable  targets.  Also,  each  aircraft  will 
usually  be  assigned  more  than  one  target  per  mission.  With 
such  long  missions,  fuel  constraints  and  other  mission  re¬ 
quirements  will  limit  the  time  available  to  search  for  each 
target.  In  addition,  the  entire  area  is  searched  regardless 
of  any  updated  information  that  may  be  received.  Therefore, 
an  exhaustive  search  is  impractical  and  a  shorter  search  is 
needed . 

Effort  Pi str i but  1  on .  Dne  search  pattern  which  requires 
less  time  is  one  which  divides  the  total  search  time  avail¬ 
able  among  the  individual  grids  to  maximize  the  probability 
of  detection  (16:5-1).  Initially  the  total  available  search 
time  is  given  and  defined  as  the  total  search  effort.  An 
assumption  is  made  that  spending  more  time  (effort)  search¬ 
ing  a  particular  grid  will  increase  the  probability  of 
detecting  the  target  if  the  target  is  in  the  grid  being 
searched.  With  this  assumption,  the  unconditional  probabili¬ 
ties  of  the  target  being  in  each  grid  are  sorted  from  high¬ 
est  to  lowest.  The  available  search  effort  is  then  divided 
among  the  grids  on  the  sorted  list  to  maximize  the  probabil¬ 
ity  of  detecting  the  target  (16:5-1).  However,  a  problem 
with  this  search  pattern  is  also  encountered. 

An  effort  distribution  algorithm  such  as  just  described 
does  not  account  for  any  movement  of  the  target  during  the 
actual  search.  One  of  the  assumptions  made  for  the  SRT 


jo 


missions  is  that  the  scenario  would  involve  bombers  attack¬ 
ing  SRT ' s  with  very  short  dwell  times.  Even  though  aircraft 
search  times  may  be  short,  targets  could  move  sometime  after 
the  search  begins  and  before  the  search  time  is  over. 
Therefore,  a  search  pattern  that  allows  target  movement 
during  the  aircraft  search  must  be  used. 

Markov  Moti on-Myop i c  Sear ch .  A  myopic  search  pattern 
can  adjust  to  target  movement  during  the  search  (16:6-1).  A 
myopic  search  is  one  which  applies  all  available  search 
effort  for  one  time  in  the  location  with  the  highest  proba¬ 
bility  of  containing  the  target.  In  so  doing,  this  search 
pattern  maximizes  the  probability  of  detection  for  a  speci¬ 
fic  point  in  time  given  the  target  has  not  been  found  in  the 
previous  time  periods  (16:6-5). 

Three  inputs  are  required  to  compute  the  unconditional 
probability  distribution  of  the  target's  location  using  the 
myopic  search  pattern.  First,  the  probability  of  the  target 
being  in  each  grid,  x,  at  a  time  t  without  being  detected  by 
any  previous  search  is  needed  (usually  this  will  represent 
the  initial  position  where  t=0) .  Second,  the  probability 
that  a  target  in  grid  x  at  time  t  will  not  be  detected  at 
time  t  must  be  known.  Finally,  the  third  input  is  the 
probability  that  a  target  in  a  grid  x  at  time  t  goes  next  to 
grid  y  (this  represents  the  transition  matrix).  Multiplying 
these  three  inputs  together  results  in  the  unconditional 
probab 1 1 1 t 1 es  of  the  target  being  in  each  possible  grid  at 


time  t+1  and  not  detected  in  the  previous  time  periods. 


From  these  pr obab 1 1 i t i es ,  the  grid  with  the  highest  proba¬ 
bility  of  containing  the  target  at  time  t+1  is  searched. 

A  myopic  search  pattern  searches  only  one  grid  each 
time  period  to  maximize  the  detection  probability.  However , 
tor  the  SRT  mission,  several  grids  will  be  searched  during 
one  time  period.  Depending  on  how  -fast  the  aircraft  flies, 
how  large  the  grids  are,  and  how  much  search  time  is  allot¬ 
ted,  any  number  of  grids  can  be  searched  each  time  period. 
The  goal  is  to  maximize  the  probability  of  detecting  the 
target  given  that  a  specified  number  of  "looks"  (grids 
searched)  are  allowed.  Therefore,  the  pattern  is  modified 
slightly. 

Only  two  of  the  previous  inputs  are  required  for  the 
aircraft  SRT  mission.  By  using  the  properties  of  a  Markov 
chain,  multiplying  the  initial  position  vector  (computed  at 
the  update  time)  times  the  n-step  transition  matrix  results 
in  the  unconditional  probability  distribution  of  the  tar¬ 
get's  location  at  the  beginning  of  the  search.  The  grids 
are  searched  from  highest  to  lowest  according  to  the  proba¬ 
bility  of  containing  the  target.  Termination  occurs  when 
the  target  is  found  or  when  the  allotted  number  of  grids 
have  been  searched.  If  a  time  period  ends  while  the  search 
is  in  progress,  new  pr obab 1 1 i t 1 es  are  computed  and  the 
search  continues  with  the  new  results. 

Using  the  modified  Markov  motion-myopic  search  pattern 
allows  the  aircraft  to  search  the  "best"  grids  first  and 


allows  adjustments  to  target  movement  while  the  search  in  is 
progress.  However,  making  the  adjustments  and  getting  the 
■first  set  of  unconditional  probabilities  requires  computing 
the  initial  position  vector  and  the  n-step  transition  ma¬ 
trix. 

Updated  I  n -for  mat  i  on 

Information  at  several  update  times  will  be  simulated 
but  no  matter  what  the  update  time,  the  required  calcula¬ 
tions  are  the  same.  Target  movement  and  the  aircraft  mis¬ 
sion  progress  independently  from  one  another.  However,  when 
an  update  time  arrives,  the  initial  position  vector  of  the 
target  as  seen  by  the  aircrew  is  reset  to  reflect  the  actual 
target  location.  The  remaining  time  until  the  aircraft 
reaches  the  search  area  is  then  computed. 

Using  the  time  unit  of  the  stochastic  process,  an 
estimate  of  the  number  of  moves  to  be  made  by  the  target  is 
calculated.  From  this  information,  the  appropriate  n-step 
transition  matrix  can  be  computed  and  subsequently  the  un¬ 
conditional  probability  distribution  of  the  target's  loca¬ 
tion  computed  for  the  aircraft's  arrival  time.  Also,  an 
estimate  is  computed  for  the  time  of  the  target's  first  move 
after  the  aircraft  begins  the  search.  If  a  move  is  antici¬ 
pated  by  the  aircrew  while  the  search  is  in  progress,  ad¬ 
justments  are  made  to  the  unconditional  probability  dis¬ 
tribution  of  the  target's  location. 
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By  ranking  the  unconditional  pr  obab  1  1  i  t  i  es  of  target 


detection  from  high  to  low  probability,  the  aircraft  can 
search  grids  with  the  highest  probability  o-f  target  presence 
-first.  Each  aircraft  mission  will  only  be  allowed  to  search 
a  limited  number  of  grids.  If  the  aircraft  searches  the 
grid  containing  the  target,  a  random  number  is  drawn  to  test 
for  target  detection  by  the  aircraft  sensors  based  on  the 
computed  conditional  target  detection  factor.  If  the  air¬ 
craft  estimates  that  the  target  is  moving  to  another  grid 
while  the  search  is  in  progress,  the  (n+l)-step  transition 
matrix  is  computed  followed  by  calculation  of  the  new  uncon¬ 
ditional  location  distribution.  After  the  new  distribution 
is  sorted,  the  search  continues  in  the  grid  at  the  top  of 
the  sorted  list. 

For  example,  assume  the  aircraft  is  initially  given  the 
following  vector  of  grids  to  search:  [4,3,6,5,1,2,7,9,83. 
In  addition,  assume  grid  4  has  been  searched  and  the  air¬ 
craft  is  starting  to  search  grid  3.  If  a  move  is  antici¬ 
pated,  new  calculations  are  completed.  Assume  further  the 
new  vector  of  grids  (sorted  from  highest  to  lowest)  is 
[9,8,6,7,5,4,3,1,23.  After  searching  grid  3  the  aircraft 
will  search  grid  9  from  the  top  of  the  new  list  rather  than 
grid  6  from  the  old  list.  In  this  way,  the  most  current 
information  on  the  target's  location  is  being  used  for  the 
search.  To  arrive  at  the  estimate  of  the  detection  proba¬ 


bility,  many  missions  are  generated  to  get  accurate  results. 


Search  Time.  Each  bomber  mission  will  be  limited  to  a 


speci-fic  amount  o-f  search  time.  As  mentioned  earlier,  cur¬ 
rently  SAC  planners  use  one  sweep  through  a  search  area  as 
the  time  limit  <5,13).  Consi derati ons  such  as  other  targets 
to  strike  and  exposure  to  enemy  defenses  come  into  play  when 
determining  the  actual  limit.  For  the  aircraft  model,  i  -f 
the  time  limit  runs  out  before  the  aircraft  finds  the  tar¬ 
get,  the  mission  records  a  miss  and  checks  are  then  done  to 
see  if  the  simulation  run  is  complete  or  if  another  mission 
should  be  generated. 

Mini  mum  Fuel .  Very  closely  related  to  the  search  time 
limit  is  the  fuel  constraint.  In  order  to  continue  the  mis¬ 
sion  and  land  safely,  an  aircraft  will  be  required  to  end  a 
search  with  no  less  than  a  specific  amount  of  fuel  reserve. 
If  more  search  time  is  available  when  the  check  is  done  but 
the  fuel  reserves  are  below  the  minimum  allowed,  the  air¬ 
craft  will  have  to  exit  the  search  area.  Again,  if  the 
target  has  not  been  found  a  miss  is  recorded  for  that  par¬ 
ticular  aircraft  mission  and  appropriate  checks  are  com¬ 
pleted  as  before. 

T ar get  Detect i on .  Finding  the  target  is  the  quickest 
way  to  terminate  the  aircraft  mission.  Two  requirements 
must  be  satisfied  before  a  hit  can  be  recorded.  Firs  ,  the 


aircraft  must  choose  to  search  a  grid  that  contains  the 
target.  If  the  correct  grid  is  selected,  the  second  re¬ 
quirement  to  be  satisfied  is  that  the  aircraft  sensors  must 
be  able  to  recognize  the  target.  To  model  the  sensors,  a 
conditional  detection  factor  is  input  as  a  probability. 
Then  a  random  number  is  drawn  to  determine  if  the  target  is 
found  (a  random  number  less  than  or  equal  to  the  detection 
factor  means  a  hit).  Appropriate  steps  are  then  taken  to 
end  the  run  or  generate  another  mission  depending  on  the 
final  checks. 

Cone 1 usi on 

Integrating  the  aircraft  model  with  the  target  model 
will  allow  an  estimate  of  the  detection  probability  to  be 
computed  over  a  range  of  intelligence  update  times.  Since 
the  target  continues  to  move  throughout  the  simulation  inde¬ 
pendent  of  the  aircraft  model,  the  timing  of  the  aircraft 
mission  is  very  important  to  ensure  the  target  is  given 
ample  time  to  move  around.  As  was  indicated  earlier,  not 
allowing  enough  movement  time  could  bias  the  detection  prob¬ 
abilities  to  the  high  side.  Also,  by  using  a  modified 
myopic  search  pattern  for  all  aircraft  missions  and  keeping 
all  other  variables  constant,  any  changes  in  the  estimates 
for  the  different  update  times  can  be  attributed  to  the 
changing  update  time  only.  To  integrate  the  two  models,  a 
computer  simulation  model  is  presented  in  the  next  chapter 
which  is  used  to  generate  the  data  to  compute  the  estimate 


■for  the  detection  probabilities.  From  these  probabi  1  i  ti  es 
the  desired  plot  o-f  detection  probabi  1  i  ti  es  versus  intelli 
gence  update  times  can  then  be  drawn. 


IV.  Simulation  Model 


Introduction 

Development  of  the  simulation  model  requires  an  accur¬ 
ate  description  of  the  real  world.  Previous  chapters  have 
explained  the  factors  affecting  both  SRT  movement  and  air¬ 
craft  SRT  missions.  Emphasis  now  turns  to  developing  a 
model  that  will  output  an  estimate  for  the  probability  of 
detecting  an  SRT  for  specified  intelligence  update  times. 
Fig.  4.1  shows  the  flow  of  the  entire  simulation  model. 

After  the  computer  languages  used  are  introduced,  Chap¬ 
ter  4  continues  in  a  chronol ogi cal  fashion  to  describe  the 
other  activities  in  the  simulation.  First,  initial  condi¬ 
tions  for  the  simulation  are  discussed.  Next,  procedures 
for  moving  the  SRT  are  explained.  Third  is  a  description  of 
how  the  aircraft  mission  is  initiated,  how  receipt  of  intel¬ 
ligence  updates  are  handled,  and  how  the  flight  to  the 
search  area  is  modeled.  Finally,  the  last  two  sections 
describe  haw  the  aircraft  search  is  modeled  and  how  the 
aircraft  mission  is  terminated. 

Computer  Languages 

SLAM.  A  Simulation  Language  for  Alternative  Modeling, 
called  SLAM,  is  the  primary  computer  language  used  for  this 
model.  SLAM  is  a  FORTRAN  based  language  that  allows  simula¬ 
tion  models  to  be  built  using  one  of  three  different  orien¬ 


tations  <9:  i  :< )  . 


In  describing  SRT  movement  and  aircraft  SRT 


low  Chart  o-f  the  Complete  Simulation 
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missions,  the  orientation  assumed  is  a  combined  network- 
discrete  event  orientation. 

A  network  orientation  or  process  orientation  models  the 
-flow  of  entities  through  a  defined  process  (9:78).  Each 
aircraft  and  each  target  represent  one  entity  or  item  in  the 
SRT  process.  The  process  is  the  sequence  of  events  and 
activities  involved  in  moving  the  target  and  flying  the 
aircraft  mission.  For  the  SRT  model,  the  network  will 
generate  the  flow  of  the  aircraft  and  target  entities  from 
the  beginning  of  the  simulation  until  one  of  the  termination 
conditions  is  met.  While  network  orientations  have  the 
advantage  of  easier  development,  not  all  systems  can  be 
accurately  represented  by  the  available  network  elements 
(9:323) . 

In  a  discrete  event  orientation,  the  model  is  built  by 
defining  the  changes  that  occur  at  discrete  points  in  time 
called  event  times  (9:229).  At  each  event  time  the  logic 
associated  with  the  occurrence  of  that  event  is  executed  in 
a  time  ordered  sequence  (9:66).  A  model  is  built  by  identi¬ 
fying  those  events  where  changes  in  the  system  state  occur 
and  developing  the  logic  to  make  the  required  changes  in  the 
system  when  the  event  does  occur.  Changing  the  state  of  an 
SRT  mission  involves  updating  the  location  of  the  target  and 
can  be  done  with  the  available  network  elements.  However, 
other  changes  require  matrix  operations  using  with  the  Har¬ 
kov  chain  procedures  and  cannot  be  done  within  the  network 
orientation.  Therefore,  to  completely  model  the  SRT  move- 


ment  and  the  aircraft  missions,  a  combined  network-discrete 
event  orientation  is  modeled. 

By  using  the  combined  -framework,  the  advantages  o-f  both 
orientations  are  enjoyed.  Portions  of  the  system  that  can 
be  described  by  a  network  are  modeled  as  a  network  with  a 
discrete  event  viewpoint  used  only  for  portions  that  require 
the  added  flexibility  (9:323).  A  reduction  in  modeling 
effort  is  realized  since  the  major  flow  of  the  system  uses 
the  network  orientation.  Only  the  basic  matrix  algebra 
requires  modeling  outside  the  network.  To  accomplish  the 
matrix  operations  requires  interfacing  the  SLAM  program  with 
some  FORTRAN  subroutines. 

FORT  RAN.  Each  discrete  event  in  the  model  is  coded  as 
a  FORTRAN  subroutine  (9:236).  As  entities  flow  through  the 
network,  the  time  will  come  for  the  occurrence  of  a  discrete 
event.  SLAM  allows  specific  coding  to  call  the  appropriate 
FORTRAN  subroutine  to  perform  the  required  functions 
(9:236).  However,  the  entire  simulation  is  controlled  by 
the  SLAM  network  including  the  simulated  clock  time. 

SLAM  advances  time  as  appropriate  and  relieves  the 
programmer  of  sequencing  events  in  their  proper  chronologi¬ 
cal  order  (9:237).  This  allows  the  programmer  to  write 
independent  units  of  code  that  can  occur  simultaneously 
while  the  simulation  is  running.  A  SLAM  executive  program 
chooses  the  "next  event"  to  execute  as  the  simulation 
progresses.  When  there  are  no  events  remaining  on  the  next 


event  "calendar",  the  simulation  is  terminated  and  a  new  run 
can  be  initiated.  However,  at  the  start  of  every  run, 
certain  actions  must  be  taken  to  initialize  the  model. 

Time  Zero  Activities 

An  SRT  is  created  at  time  zero  and  the  movement  started 
(Fig.  4.2  shows  the  flow  of  the  target  model).  The  center 
grid  of  the  search  box  is  identified  based  on  the  input  data 
and  the  target  is  placed  in  the  center  position.  An  initial 
target  location  vector  is  developed  to  indicate  the  target's 
postion.  Also,  since  no  dwell  time  is  assumed  for  the  first 
move,  a  random  direction  is  chosen  and  a  delay  time  initia¬ 
ted  for  the  time  required  to  move  the  target  to  the  boundary 
between  the  old  grid  and  the  new  grid.  As  the  target  s 
position  is  initialized,  the  aircraft  prepares  for  the  mis¬ 
sion  . 

A  decision  must  be  made  about  the  first  event  to  occur 
in  the  aircraft  model  before  any  activities  can  occur  for 
this  segment.  The  first  event  will  either  be  the  aircraft 
takeoff  or  the  intelligence  update  (the  flow  of  the  aircraft 
model  is  shown  in  Fig.  4.3).  If  the  takeoff  occurs  -first,  a 
delay  must  be  initiated  to  account  for  the  time  until  the 
scheduled  takeoff  followed  by  a  delay  for  the  flight  time  to 
the  update  point  with  a  final  delay  for  the  flight  time  to 
the  search  area.  Otherwise,  a  delay  is  initiated  for  the 
time  until  the  update  is  received  while  still  on  the  ground 
followed  by  other  delays  for  the  remaining  ground  time 


45 


and  the  -flight  time  -from  the  takeoff  point  to  the  search 


area.  As  the  aircraft  model  continues,  the  target  moves 
independently  of  any  aircraft  actions. 

Continuous  T ar get  Movement 

Before  any  move  can  be  made,  the  relative  position  of 
the  SRT  inside  the  defined  search  "box"  must  be  identified. 
A  random  direction  is  selected  for  the  next  move  to  be  made 
by  the  SRT.  However,  depending  on  the  position  of  the 
target  relative  to  the  search  box,  either  eight,  five,  or 
three  directions  will  be  available  to  choose  from.  In  the 
center  area  of  the  box,  eight  grids  surround  the  current 
location  and  therefore  one  of  eight  directions  can  be  selec¬ 
ted  (Fig.  4.4a).  But  if  the  target  is  on  the  boundary  of 
the  box,  only  one  of  five  directions  can  be  selected 


Fig.  4.4.  Three  Relative  Positions  Possible  for  the 
Target  in  the  Search  Area 
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because  movement  outside  the  box  is  prohibited  (Fig.  4.4b>.  / 

Even  more  restrictive,  if  a  corner  position  is  occupied,  • 


only  one  of  three  directions  can  be  selected  (Fia.  4.4c) 
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Determining  the  relative  position  of  the  target  must  come 
before  a  direction  can  be  chosen.  To  determine  the  posi¬ 
tion,  further  screening  must  be  accomplished. 

For  SRT ' s  located  on  the  boundary  of  the  box ,  further 
tests  must  be  done  to  determine  exactly  which  edge  or  which 
corner  is  occupied.  Imagine  the  search  box  in  a  north-south 
orientation  as  depicted  in  Fig  4.5a.  From  the  center  grid, 
eight  possible  directions  can  be  selected,  each  representin 


Fig.  4.5.  Possible  Directions  for  the  Target  to  Move 
a  compass  direction  as  labeled  in  the  figure.  Now  consider 
the  situation  in  Fig.  4.5b.  On  the  "edge"  of  the  box  only 
five  directions  are  available.  Moving  to  Fig  4.5c,  in  a 
corner  position  only  three  directions  are  available  (Note: 
To  be  consistent,  the  box  is  always  oriented  in  a  N-S  direc¬ 
tion  and  the  eight  compass  points  remain  the  same).  Each  of 
the  four  corners  and  each  of  the  four  edges  have  the  same 
number  of  possible  move  directions  but  have  a  unique  set  of 
actual  move  directions  (Fig.  4.6).  Therefore,  to  determine 
which  directions  are  possible  from  the  current  position,  the 
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Fig.  4.6.  Contrast  of  Relative  Positions  and 
Choice  of  Directions 

exact  location  -for  boundary  and  corner  grids  must  be  identi¬ 
fied. 

Once  the  exact  location  o-f  the  occupied  grid  has  been 
identified  a  new  direction  can  be  chosen.  Branching  within 
the  simulation  model  will  place  the  target  entity  in  the 
appropriate  section  of  code  associated  with  its  current 
position.  A  random  number  from  a  uniform  distribution  will 
be  chosen.  The  range  of  the  distribution  will  be  zero  to 
three,  five,  or  eight  depending  on  the  situation.  Each  of 
the  appropriate  compass  directions  is  associated  with  a  sub¬ 
range  of  the  distribution  and  a  direction  will  be  chosen 
based  on  the  value  of  the  random  number. 

For  example,  if  five  directions  are  available,  the 
direction  parameters  will  be  0  to  5.  A  random  number  less 
than  one  will  refer  to  one  direction,  between  one  and  two 
another  direction  and  so  on  until  all  five  directions  are 
covered.  Once  the  new  direction  is  chosen,  the  actual  grid 
to  be  entered  is  computed.  But  before  the  move  can  be  made, 
the  status  of  the  new  grid  must  be  checked. 

A  test  must  be  done  on  the  future  location  to  see  if  it 


is  suitable  for  the  move. 


With  the  input  data,  each  grid  is 


assigned  a  status  of  0,  1,  or  2:  0  means  no  passage,  1  means 
the  target  can  pass  through  but  not  stop,  and  2  means  the 
target  can  pass  through  and  stop.  If  the  new  direction  will 
put  the  SF:T  in  a  grid  with  a  0  status,  a  new  direction  is 
chosen  before  any  move  takes  place  (no  simulated  time  elap¬ 
ses).  If  a  status  of  1  is  discovered,  the  move  is  made  with 
the  appropriate  amount  of  time  passing  but  another  move  is 
initiated  before  any  dwell  time  is  applied  (for  a  mobile 
target,  only  0's  and  Is  apply  since  mobile  SRT ' s  have  0 
dwell  time).  Finally,  if  a  grid  with  a  status  of  2  is 
selected  the  move  can  be  made  with  a  delay  initiated  repre¬ 
senting  the  appropriate  move  time  and  dwell  time.  At  the 
completion  of  the  delay,  the  process  repeats  and  continues 
until  the  termination  of  the  mission. 

Aircraft  Update  Procedures 

While  the  SRT  movement  is  being  simulated,  an  aircraft 
mission  continues  according  to  the  scheduled  timing.  As 
mentioned  earlier,  the  aircraft  will  follow  a  certain  se¬ 
quence  depending  on  whether  the  takeoff  or  update  occurs 
first.  In  either  case,  time  delays  are  built  into  the  model 
to  allow  for  an  appropriate  amount  of  random  target  movement 
up  to  the  update  time.  But  the  update  procedures  are  the 
same  no  matter  what  time  the  update  occurs. 

Simulating  update  procedures  requires  revealing  to 


the  aircrew  the  current  location  of  the  target. 


When  the 


update  event  occurs,  the  initial  target  location  vector  is 
changed  to  reflect  the  new  position.  Then,  to  compute  the 
unconditional  probability  distribution  o-f  the  target's  loca¬ 
tion  for  the  aircraft  arrival  time,  the  appropriate  n-step 
transition  matrix  must  be  calculated.  Since  the  one-step 
transition  matrix  is  raised  to  th  nth  power,  to  compute  the 
n-step  transition  matrix  the  value  of  n  must  be  known. 

To  -find  the  value  o-f  n,  the  anticipated  number  of 
moves  the  target  will  make  before  the  aircraft  arrives  is 
calculated.  At  the  time  of  the  update,  intelligence  sources 
do  not  know  how  long  the  SRT  has  been  in  the  current  loca¬ 
tion.  To  cut  down  on  some  of  the  error,  an  assumption  is 
made  that  the  target  has  been  in  the  reported  location  for 
one  half  of  the  total  time  that  is  normally  spent  in  one 
location.  The  aircrew  will  then  anticipate  a  move  after  one 
half  of  the  normal  time  has  elapsed.  After  the  first  move, 
all  future  moves  will  be  anticipated  after  a  full  time  unit. 
Dividing  the  time  remaining  until  the  aircraft  arrives  (mi¬ 
nus  the  time  for  the  first  move)  by  the  full  time  unit  will 
give  (n-1)  which  is  needed  to  compute  the  appropriate  n-step 
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An  unconditional  location  distribution  must  be  computed 
to  guide  the  aircraft  in  the  search  for  the  uPT.  Multi¬ 
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transition  matrix  will  result  in  the  unconditional  location 
distribution  of  the  target  -for  the  aircra-ft  arrival  time. 
The  grids  are  then  sorted  -from  high  to  low  probabilities. 
Sorting  mates  programming  the  aircra-ft  search  easier  because 
when  the  grids  are  selected  from  top  to  bottom  on  the  or¬ 
dered  list  the  highest  probability  grids  are  searched  -first.. 
When  all  update  calculations  are  complete,  the  aircra-ft 
mission  continues  on  its  scheduled  time  until  the  aircra-ft 
arrives  in  the  search  area. 

Aircraft  Search 

A  search  begins  by  moving  the  aircra-ft  to  the  first 
grid  on  the  sorted  list.  Since  the  maximum  probability  of 
finding  the  target  is  desired,  searching  the  grids  with  the 
highest  probability  of  containing  the  SRT  first  will  achieve 
the  desired  goal  over  the  long  term.  A  check  is  done  to  see 
if  the  aircraft  location  (grid  number)  is  the  same  as  the 
target's  actual  location.  If  the  locations  are  not  the 
same,  the  process  repeats  for  the  next  grid  on  the  sorted 
list.  However,  if  the  locations  are  the  same,  a  check  is 
done  to  see  if  the  aircraft  sensors  identify  the  target. 

Sensor  systems  are  checked  by  comparing  a  random  number 
to  the  probability  of  detection  associated  with  the  sensors. 
Since  the  sensors  are  not  perfect,  inputs  into  the  model  can 
specify  a  probability  distribution  to  account  for  errors  in 
the  system.  A  random  number  is  drawn  and  compared  to  the 
conditional  detection  probability  to  determine  the  success 


of  the  sensors  -for  that  particular  “look"  at  the  target.  I-f 
a  success  occurs,  a  hit  is  recorded  -for  that  aircraft  mis¬ 
sion.  Failure  of  the  sensors  is  indicated  by  a  message 
telling  the  decision  maker  that  the  aircraft  found  the 
correct  grid  but  still  missed  the  target  due  to  the  sensors. 
In  this  case,  the  search  continue s. 

To  continue  the  search,  the  next  grid  on  the  sorted 
list  is  selected  by  the  aircraft.  Before  a  move  is  made 
into  the  next  grid  a  check  is  made  to  see  if  the  probability 
of  the  target's  presence  is  greater  than  0.  If  not,  the 
grid  is  not  searched  because  time  would  be  wasted  in  the 
long  term.  Instead  the  aircraft  returns  to  the  first  grid 
on  the  sorted  list  and  searches  from  the  top  of  the  list 
down  again.  But,  if  the  next  grid  on  the  sorted  list  has  a 
probability  greater  than  0  then  the  search  is  conducted  as 
before.  As  the  search  continues,  if  the  time  anticipated  by 
the  aircrew  for  the  next  target  move  occurs,  updated  infor¬ 
mation  must  be  calculated. 

New  unconditional  probabilities  must  be  computed  if  the 
target  moves  during  the  aircraft  search.  When  the  update  on 
the  target ' s  location  is  received  by  the  aircrew  while 
enroute  to  the  search  area,  a  guess  is  made  as  to  how  long 
the  target  has  been  in  the  current  location.  At  the  same 
time  an  estimate  is  made  about  the  time  of  the  target's 
first  move  after  arrival  of  the  aircraft.  The  target  may  or 
may  not  move  at  that  time  but  the  aircraft  anticipates  a 


move  based  strictly  on  timing.  If  the  time  of  the  antici¬ 
pated  move  occurs  during  the  search,  the  <n+l)-step  transi¬ 
tion  matrix  must  be  calculated  and  multiplied  by  the  initial 
position  vector  calculated  when  the  intelligence  update  was 
received.  Sorting  the  resulting  vector  will  give  a  new 
unconditional  location  distribution  with  which  to  continue 
the  search  -for  the  SRT.  Starting  at  the  top  of  the  list, 
the  aircraft  continues  as  before  until  one  of  the  termina¬ 
tion  conditions  is  met. 

Ter  mi nat i on  Procedures 

Three  conditions  will  terminate  one  mission  of  the 
simulation.  First,  if  the  aircraft  finds  the  SRT  a  hit  is 
recorded  and  the  mission  is  terminated.  Second,  after  each 
grid  is  searched,  the  allotted  search  time  is  checked  to  see 
if  it  has  been  exceeded.  If  too  much  search  time  has  been 
used,  a  miss  is  recorded  and  the  mission  is  terminated; 
otherwise  the  mission  is  continued.  Finally,  the  third 
condition  checks  the  fuel  status  of  the  aircraft.  If  the 
fuel  level  is  below  a  specified  minimum  value,  the  aircraft 
must  terminate  the  search  and  a  miss  is  recorded.  But,  the 
mission  continues  if  the  fuel  is  above  the  minimum  accept- 
ab  1  e . 

In  all  cases,  when  the  mission  is  terminated,  the 
procedures  are  the  same.  A  specified  number  of  missions  are 
run  for  each  update  time  being  tested.  When  one  mission  is 
terminated  a  check  is  made  to  see  if  it  is  the  last  mission 


■for  the  current  simulation  run.  If  not,  the  aircraft  and 
target  entities  are  initialised  to  their  starting  positions 
and  another  mission  started.  After  the  last  mission  for  one 
update  time  is  terminated,  a  new  update  time  is  computed  and 
another  set  of  missions  is  run.  Finally,  when  the  last 
update  time  has  been  tested,  the  simulation  is  terminated. 

Concl usi on 

After  integrating  the  aircraft  mission  and  target  move¬ 
ment  into  one  simulation  model,  computer  runs  can  be  made  to 
generate  the  desired  results.  Timing  is  very  important 
throughout  the  simulation  to  ensure  the  results  are  as 
realistic  as  passible.  However,  one  thing  to  keep  in  mind 
is  that  the  simulation  may  not  exactly  match  the  real  world. 
An  exact  match  is  virtually  impossible  to  achieve.  A  simu¬ 
lation  model  is  only  as  good  as  the  inputs  into  the  model 
and  the  assumptions  made  while  building  the  model.  Verifi¬ 
cation  and  validation  of  the  model  is  very  important  if  the 
output  from  the  computer  runs  is  to  be  of  any  use.  Chapter 
5  discusses  how  the  verification  and  validation  stages  were 
accomp 1 i shed . 
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V.  Verification  and  Validation 


I ntr oduc t i on 

Computers  in  general  do  not  make  mistakes;  people  (pro¬ 
grammers)  make  mistakes.  Simulation  models  are  not  worth 
the  time  and  effort  to  run  if  the  model  is  not  a  reasonable 
representation  of  the  real  world.  Verification  and  valida¬ 
tion  of  the  computer  model  are  required  if  the  model  is  to 
provide  any  useful  information  to  a  decision  maker. 

Chapter  5  discusses  the  verification  and  validation 
stages  of  this  simulation  model.  Verification  is  accom¬ 
plished  by  comparing  the  model  output  to  manual  calculations 
for  a  small  example  mission.  To  validate  the  model,  input 
values  are  varied  and  changes  in  the  final  results  are 
compared  to  expected  values  based  on  knowledge  and  past 
experience  of  the  SRT  missions. 

Ver 1 f 1  cat  1  on 

Verification  of  the  model  consists  of  determining  that 
the  model  executes  as  intended  (9:12).  For  the  target 
movement,  an  output  trace  of  all  grid  locations  which  the 
target  occupies  and  the  times  of  occupation  are  output.  A 
close  look  at  the  data  indicates  if  moves  are  being  made  in 
all  directions  acording  to  the  uniform  distribution  assigned 
and  if  the  moves  occur  at  the  proper  times.  Also,  a  check 
is  made  to  see  if  other  location  restrictions  are  being 


observed . 


To  verify  the  aircraft  mission, 


a  comparison  is  made 


between  the  results  from  manual  calculations  for  the  matrix 
operations  and  the  computer  output  to  ensure  the  program  is 
doing  the  calculations  correctly.  Moreover ,  checking  the 
simulation  time  of  certain  events  throughout  the  aircraft 


mission  shows  if  appropriate  delays  are  being  integrated 
with  the  target  movement.  Output  of  an  example  using  a 
square  search  area  divided  into  nine  total  grids  (target 
locations)  will  show  the  desired  results  to  verify  the  model 
(Fig  5.1). 
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Fig.  5.1.  Search  Area  Used  to  Verify  and 
Validate  the  Model 


T ar get  Model .  To  be  consistent,  the  target  is  modeled 
to  start  every  mission  in  the  same  grid.  As  stated  earlier, 
each  target  is  assigned  to  an  MOB  (Main  Operating  Base)  with 
the  center  grid  of  the  search  area  containing  the  MOB  (Fig. 
5.1).  Before  each  mission  begins,  the  center  grid  is  com¬ 
puted  from  the  other  input  data  and  the  SRT  is  established 
in  that  grid.  Using  the  present  example,  the  output  shows 
the  target  starting  in  grid  5  (center  grid)  for  every  mis¬ 


sion  which  confirms  the  calculation  is  being  done  correctly. 


Once  initialized  to  the  proper  starting  point,  random  move¬ 
ment  can  then  begin. 

Target  movement  is  modeled  to  be  evenly  distributed  in 
all  possible  directions.  Because  only  a  small  number  of 
missions  were  simulated  for  verification  purposes,  the  per¬ 
centage  of  time  each  direction  was  chosen  will  not  be  a  very 
accurate  estimate  of  the  distribution.  However,  other  indi¬ 
cations  show  all  possible  directions  are  being  used.  Trac¬ 
ing  the  target's  location  shows  that  each  grid  is  selected 
as  a  new  location  several  times  during  the  short  simulation. 
Over  many  runs,  each  direction  should  be  chosen  an  approxi¬ 
mately  equal  number  of  times.  However,  other  location  re¬ 
strictions  prevent  same  grids  from  being  selected  as  many 
times  as  they  otherwise  would. 

Each  grid's  status  will  determine  if  a  grid  can  be 
entered  or  not  and  could  also  prevent  surrounding  grids  from 
being  selected  as  often.  As  a  new  grid  is  selected,  spe¬ 
cific  output  statements  are  printed  to  show  the  grid's 
status.  A  trace  of  the  target  movement  indicates  the  loca¬ 
tion  restrictions  are  being  observed  in  the  present  example; 
only  grids  with  a  status  of  2  have  targets  setting  up,  grids 
with  a  status  of  1  are  being  entered  with  no  dwelling,  and 
grids  with  a  status  of  0  are  not  entered  at  all.  Since 
grids  3  and  8  (fig.  5.1)  have  a  status  of  0,  grids  6  and  9 
are  selected  less  often  than  other  grids  because  more  move¬ 
ment  can  take  place  in  the  opposite  corner  (grids  1,  2,  4, 
and  5).  However,  all  grids  are  being  selected  at  some  point 


regardless  of  the  other  location  restrictions  which  leaves 
only  the  timing  o-f  the  moves  to  be  checked. 

Based  on  the  input  values  for  target  speed  and  dwell 
times,  target  moves  can  be  manually  checked  for  correct 
timing.  SRT ' s  start  the  simulation  (time  0)  in  the  center 
grid.  For  this  example,  0.5  miles  per  minute  is  used  as  the 
target  speed,  90  minutes  as  the  dwell  time,  and  17.6  miles 
as  the  grid  size.  Initially  the  target  should  only  move 
halfway  across  the  center  grid  (no  dwell  time  for  the  first 
move).  This  is  being  done  properly  in  the  model  as  confir¬ 
med  from  the  sample  output  showing  grid  1  (the  first  new 
location,  second  grid)  being  entered  at  34. 8  simulated  min¬ 
utes.  Next,  a  move  is  made  to  grid  4  at  time  194.4.  Manual 
calculations  confirm  that  with  90  minutes  of  dwell  time  and 
34.8  minutes  to  move  halfway  across  one  grid,  194.4  is  the 
correct  time  to  enter  the  third  grid.  Repetitions  of  the 
direction  selection  and  corresponding  moves  continue 
throughout  the  simulation  run  and  all  times  are  calculated 
correctl y . 

Verification  of  the  target  model  is  now  complete.  Ran¬ 
dom  directions  are  being  chosen,  location  restrictions  are 
being  observed,  and  timing  is  being  simulated  properly. 
However,  for  the  simulation  results  to  be  of  value,  the 
aircraft  model  must  also  be  verified. 

Ai rcr af  t  Model .  Verification  of  the  aircraft  model 
must  start  by  ensuring  appropriate  delays  are  accounted  for 
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before  the  first  aircraft  event  (takeoff  or  intelligence 
update).  Depending  on  the  input  information,  the  aircrew 
can  receive  the  intelligence  update  either  before  or  after 
takeoff  but  the  search  still  begins  at  the  same  simulated 
time  (when  the  aircraft  enters  the  search  area)  every  mis¬ 
sion.  Comparing  runs  with  updates  before  takeoff  with  runs 
having  the  update  after  takeoff  reveals  the  same  time  of 
860.0  simulated  minutes  as  the  time  the  search  begins  for 
the  first  mission  of  each  run.  Since  random  numbers  play  a 
part  in  some  delays,  the  start  of  the  search  for  the  first 
mission  of  future  runs  will  not  be  the  same.  However,  by 
checking  the  time  each  subsequent  mission  starts  and  adding 
the  scheduled  delays  from  time  zero  to  the  time  the  search 
starts  (860.0  for  this  example),  the  time  the  search  begins 
can  be  confirmed. 

For  example,  mission  number  1  in  run  number  1  began  at 
time  0  and  started  the  search  at  time  860.0.  A  second 
mission  started  when  all  entities  for  the  first  mission  were 
destroyed.  This  occurred  at  time  902.4.  Therefore,  the 
search  in  the  second  mission  should  begin  at  time  1762.4 
(860.0  +  902.4).  Since  output  from  the  example  confirms  the 
time,  correct  calculations  again  are  being  done. 

For  mission  number  1  in  run  number  2,  the  times  were 
the  same  as  before  except  for  the  ending  time.  Because  of 
randomness  in  some  delays,  the  first  mission  ended  at  time 
870.0  (vs  902.4).  Adding  860.0  to  870.0  indicates  the 
search  for  the  second  mission  should  begin  at  time  1730.0 


which  is  also  con-firmed  by  the  output  information.  Even 
though  the  times  are  slightly  different,  the  appropriate 
delays  up  to  the  aircraft  search  time  are  being  handled 
correctly.  But  even  though  the  search  starts  at  the  correct 
time,  the  aircrew  must  also  have  the  correct  information  in 
order  to  search  for  the  target  properly. 

Updated  information  must  be  given  to  the  aircrew  at  the 
appropriate  time  if  future  calculations  are  to  be  correct. 
If  an  80  minute  update  is  being  tested,  the  aircrew  should 
know  the  location  of  the  SRT  80  minutes  before  the  search 
begins.  Comparing  the  target  movement  trace  to  the  output 
of  the  aircraft  model  reveals  that  in  all  cases  correct 
information  is  received  by  the  aircrew. 

As  an  example,  in  mission  number  1  of  run  1,  an  80 
minute  update  is  expected  with  the  search  to  be  pin  at  860.0. 
Manual  calculations  show  that  at  time  of  780.0  the  aircrew 
should  receive  the  intelligence  update  on  the  target's  ac¬ 
tual  location.  Target  output  shows  the  target  location  as 
grid  1  at  time  780.0.  The  aircraft  model  output  confirms 
the  aircrew  received  information  that  the  target  was  in  grid 
1  at  the  80  minute  update.  This  thought  process  can  be  done 
for  all  missions  and  the  results  are  confirmed  in  c 11  cases. 
However,  yetting  the  correct  grid  number  at  the  correct  time 
to  the  aircrew  is  only  part  of  the  updated  information 
required  to  start  the  search. 


Matrix  operations  to  compute  the  unconditional  proba¬ 
bility  distribution  of  the  target's  location  must  be  done 
before  the  aircraft  search  begins.  Computing  the  uncondi¬ 
tional  probabilities  requires  the  initial  position  (or  posi¬ 
tion  at  the  update  time)  and  the  n-step  transition  matrix. 
If  no  move  can  be  made  by  the  target  between  the  time  the 
update  is  received  by  the  aircrew  and  the  time  the  search 
begins,  the  grid  number  passed  in  the  intelligence  update 
will  be  the  only  grid  searched.  If  one  or  more  moves  are 
possible,  the  one-step  transition  matrix  must  be  raised  to  a 
power  equal  to  the  estimated  number  of  moves.  Multiplying 
the  n-step  transition  matrix  times  the  initial  position 
vector  results  in  the  desired  unconditional  probab i 1 1 t i es . 
Verification  of  these  operations  can  be  done  as  before  by 
comparing  manual  calculations  with  the  computer  output  re¬ 
sulting  from  the  aircraft  search. 

Sorting  the  grid  number  by  unconditional  probability 
from  highest  to  lowest  results  in  a  ready  made  ordered  list 
of  grids  for  the  aircraft  to  search.  After  computing  the 
unconditional  probabilities,  a  sort  routine  is  used  to  de¬ 
velop  the  ordered  list.  Aircraft  will  then  enter  the  first 
grid  number  on  the  list  to  check  for  the  target.  If  the 
target  is  not  found,  the  next  grid  on  the  list  is  checked 
only  if  the  unconditional  probability  for  target  presence  is 
greater  than  0.  The  search  continues  in  this  fashion  until 
the  target  is  found  or  the  allotted  search  time  runs  out. 
Output  information  from  the  aircraft  model  indicates  that 


Since  an  80  minute  update  is  used,  the  target  is  estimated 
to  be  able  to  make  only  one  move  a-fter  the  update  is  re¬ 
ceived  by  the  aircrew  and  before  the  search  begins  (Assuming 
the  SRT  has  been  in  place  for  half  of  its  allotted  time,  the 
target  should  move  in  159.6/2  or  79.8  minutes  which  is  less 
than  80  minutes).  Therefore,  the  one-step  transition  matrix 
should  be  raised  to  the  first  power  and  multiplied  times  the 
position  vector  for  grid  1  (the  location  of  the  target  at 
the  update  time).  Results  should  reveal  the  first  row  of 
the  one-step  transition  matrix.  After  sorting  from  highest 
to  lowest  unconditional  probablity,  the  results  are  exactly 
as  expected.  However,  this  case  shows  the  sorting  routine 
worked  but  does  not  indicate  if  the  matrix  multiplication 
routine  works  because  the  one-step  transition  matrix  re¬ 


quired  no  changes  from  the  input  values. 
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Using  a  240  minute  update  in  run  number  3, 
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are  possible  by  the  target  from  the  time  of  the  update  until 
the  search  begins.  When  the  update  is  received,  the  assump¬ 
tion  that  the  target  has  spent  half  the  total  time  in  the 
current  grid  equates  to  a  grid  change  at  time  79.8.  After 
dwelling  and  moving  through  the  second  grid  (using  159.6 
minutes),  the  third  grid  will  be  entered  at  an  estimated 
time  of  239.4.  Since  the  two  moves  will  be  made  in  239.4 
minutes  and  the  aircraft  will  not  arrive  until  240  minutes 
after  the  update,  the  one-step  transition  matrix  must  be 
squared . 

Results  from  mission  number  1  on  run  number  3  are  as 
foil ows: 


Grid  Number 


Frob  of  Target  Presence 
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(6) 


For  this  mission,  the  target  was  observed  in  grid  5  at  the 
update  time.  Therefore,  the  expectation  would  be  that  the 
fifth  row  of  the  two-step  transition  matrix  (sorted)  would 
be  output.  Manual  calculations  result  in  the  following  for 
the  two  step  transition  matrix: 
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Row  number  -five,  after  being  sorted,  produces  the  same 
results  as  the  program's  calculations.  Checks  of  other 
missions  also  show  that  the  computations  and  sorting  are 
being  performed  correctly.  However,  given  the  aircrew  has 
the  correct  information ,  a  check  must  be  done  to  see  it  the 
search  is  being  done  as  intended. 

Comparing  the  sorted  unconditional  probability  list 
with  the  actual  grids  searched  by  the  aircraft  verifies  that 
the  aircraft  is  searching  the  correct  grids  in  the  right 
order.  From  mission  10  of  run  number  2  (a  random  choice) , 
results  are  the  following: 

T ar get  Locat 1  on  Pi str 1  but  1  on : 

Grid  Number  Frob.  of  Target  Presence 


Aircraft  Search: 


Grid  Searched 
4 

1  (9) 
9 

Since  three  grids  can  be  searched  in  the  allotted  time,  only 
the  -first  three  grids  on  the  sorted  list  were  searched  and 
in  the  proper  order.  A  question  that  comes  to  mind  is  which 
grid  would  be  searched  next  if  four  grids  could  be  searched. 

Rather  than  search  a  grid  that  is  estimated  to  have  a  0 
probability  of  containing  the  target,  the  search  reverts  to 
the  top  of  the  list  when  a  probability  of  zero  is  encoun¬ 
tered.  In  the  mission  described  above,  since  the  next  grid 
(grid  2)  has  a  probability  of  zero  of  target  presence,  grid 
4,  at  the  top  of  the  list,  should  be  revisited.  Another 
mission  illustrates  the  handling  of  this  type  of  situation. 

Returning  to  mission  number  1  of  run  number  1,  the 
search  returned  to  the  first  grid  on  the  sorted  list  rather 
than  search  a  grid  with  a  0  probability  of  target  presence. 
Grids  4  and  5  were  searched  first  because  each  had  a  proba¬ 
bility  of  0.5  of  containing  the  target.  However,  all  other 
grids  were  estimated  to  have  a  probability  of  zero  of  con¬ 
taining  the  target.  Therefore,  after  searching  grid  5  with 
no  detection,  the  third  grid  searched  was  4  again  rather 
than  1.  Output  information  from  the  run  is  as  follows: 


tion  must  be  calculated  using  the  location  con-firmed  at  the 
update  time.  Finally,  a-fter  sorting  the  probability  distri¬ 
bution,  the  search  continues  in  the  -first  grid  on  the  new 
sorted  list.  However,  one  critical  point  at  this  stage  is 
that  the  search  time  used  by  the  aircraft  is  cumulative.  A 
new  search  does  not  start.  The  same  search  continues  using 
the  new  sorted  list.  If  two  grids  were  already  searched  and 
time  only  allows  a  total  of  three  to  be  searched,  only  one 
grid  can  be  searched  from  the  new  list.  Verification  of  the 
model  in  these  circumstances  is  accomplished  using  the  same 
comparisons  as  before. 

Another  example  run  with  an  180  minute  update  gave  the 
following  results  (target  was  in  grid  5  at  the  update  time): 

Initial  Location  Distribution: 


Grid  Number 


F'rob.  of  Target  Presence 


After  the  move  the  new  location  distribution  was: 


Grid  Number 


Frob.  of  Target  Presence 
0 . 65 
0.  2 
0.  15 
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comparison  of  the  model  and  system  is  done  as  a  test  -for 


reasonableness  (9:13).  Past  system  outputs  and  knowledge  o-f 
system  performance  behavior  based  on  experience  should  be 
used  as  yardsticks  for  the  comparison  (9:13). 

Data  Inputs.  Information  concerning  strategic  reloca¬ 
table  targets  is  classified.  Char acter i st 1 cs  such  as  dwell 
time,  set  up  time,  tear  down  time,  and  possible  setup  loca¬ 
tions  are  all  available  but  only  in  classified  form  and 
cannot  be  addressed  in  an  unclassified  project.  For  verifi¬ 
cation,  notional  or  sample  numbers  in  the  range  of  actual 
values  can  be  used.  However,  results  that  are  output  from 
the  model  using  the  notional  numbers  cannot  be  used  for 
validity  checks  since  experience  with  SRT  missions  is  lim¬ 
ited.  But,  the  results  can  be  used  by  doing  some  sensi¬ 
tivity  analysis  with  various  input  values  in  conjunction 
with  the  validation  of  the  model  elements. 

Model  El ements .  As  with  the  data  inputs,  the  output 
from  the  computer  program  will  not  validate  the  model  ele¬ 
ments  but  can  indicate  a  reasonable  representation  of  the 
real  world  by  way  of  some  sensitivity  analysis.  Outputs 
from  past  systems  are  not  available  because  the  area  of 
SRT's  is  still  quite  new.  However,  assessments  can  be  made 
about  the  reasonableness  of  the  model  by  varying  the  input 
values  and  observing  the  changes  in  the  results  (9:13). 
Variations  in  four  factors  were  observed. 

First,  the  grid  size  was  varied.  In  all  cases,  the 
probability  of  finding  the  target  was  lower  or  unchanged 
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when  the  grid  size  was  increased.  Since  more  area  has  to  be 
covered  with  the  larger  grid  size,  this  is  a  logical  result. 
Aircraft  search  time  is  limited  and  i  -f  more  time  is  spent 
searching  one  grid,  less  time  remains  to  search  other  grids. 
Also,  a  larger  grid  size  allows  more  room  for  the  target  to 
hide.  Therefore,  responses  to  changes  in  grid  size  appear 
valid. 

Search  time  is  the  second  factor  varied.  As  would  be 
expected,  given  more  search  time  an  increase  in  the  proba¬ 
bility  of  finding  the  target  is  observed.  If  more  grids  can 
be  searched,  chances  are  better  of  finding  the  SRT. 

Third,  the  time  unit  of  the  stochastic  process  is 
changed.  The  time  unit  is  the  total  time  an  SRT  spends  in 
one  grid  (including  dwell  time,  set  up  time,  tear  down  time 
and  moving  time).  For  this  variation,  the  pattern  was  not  as 
obvious.  However,  the  overall  indication  is  that  as  the 
time  unit  is  increased  (SRT  spends  more  time  in  one  loca¬ 
tion)  ,  the  probability  of  detection  also  increases.  Again  a 
logical  conclusion. 

Finally,  the  fourth  factor  to  be  varied  is  the  aircraft 
sensor  system  detection  factor.  Logically,  as  sensor  sys¬ 
tems  improve,  the  chance  of  finding  the  target  with  them 
should  be  better.  In  each  case,  the  model  output  showed  an 
increased  or  unchanged  probability  of  detection  indicating 


again  the  correct  sensitivity  of  the  model  to  various 


in- 


Subsystems  and  Interface  Points.  Validation  of  the 
subsystems  and  interface  paints  has  really  already  been 
addressed  as  they  pertain  to  this  model.  Each  subsystem  of 
the  model  involved  specific  calculations  which  have  already 
been  verified.  Also,  in  verfifying  the  model,  the  timing  of 
events  was  shown  to  be  correct.  Each  event  in  the  SLAM 
model  called  the  appropriate  FORTRAN  subroutines  at  the 
proper  time  to  perform  the  required  actions.  Validation  of 
the  'nterface  between  SLAM  and  FORTRAN,  as  well  as  between 
the  target  and  aircraft  models  is  complete  due  to  the  veri¬ 
fication  of  the  timing  and  validation  of  the  model  elements. 

Cone  1 usi on 

Verification  and  validation  of  a  model  are  necessary 
stages  of  model  building.  Verification  proves  the  model 
does  what  is  intended.  Validation  goes  one  step  farther  to 
ensure  what  is  intended  is  a  reasonable  representation  of 
the  real  world  system  being  simulated.  Verification  was 
accomplished  in  this  model  by  comparisons  with  manual  calcu¬ 
lations  using  a  small  example.  Validation  was  difficult 
because  of  no  past  data.  However ,  a  sensitivity  analysis  on 
the  inputs  into  the  model  indicated  the  results  shifted  in 
logical  directions  as  specific  factors  were  varied.  The 
numbers  themselves  were  not  important  since  only  notional 
numbers  were  input.  But  what  is  important  is  that  the 
di r ect i ore 


the  probability  of  detection  moved  as  the  inputs 


varied  are  logical  and  give  a  good  indication  of  the  valid¬ 
ity  of  the  model. 

A  primary  purpose  of  this  project  was  to  build  a  simu¬ 
lation  model  for  an  SRT  mission;  but  this  is  not  the  end. 
With  a  valid  model,  large  examples  can  be  run  and  further 
analysis  can  be  accomplished  on  the  results.  Even  more,  the 
value  of  the  output  information  and  how  the  information  can 


be  used  will  be  investigated  in  the  final  chapters. 


VI.  Experimentation  and  Resul ts 

Introduct i on 

A  representati ve  SRT  mission  was  simulated  to  demon¬ 
strate  how  this  model  can  be  used.  To  keep  the  thesis 
unclassified,  sample  numbers  within  the  range  of  actual 
values  were  input.  Subsequent  analysis  on  the  output  of  the 
simulation  runs  using  the  sample  inputs  illustrates  how  the 
model  could  be  adequately  used  in  a  classified  environment. 

Chapter  6  first  describes  the  factorial  design  used  for 
initial  screening  of  the  factors.  Analysis  of  the  results 
follows  to  indicate  which  factors  or  interaction  effects  are 
considered  significant.  Using  the  results  of  the  factorial 
design,  final  "production"  runs  are  conducted  and  their 
results  analysed. 

Exper i mental  Desi gn 

Many  experiments  involve  the  study  of  the  effects  of 
two  or  more  factors  on  a  desired  result.  A  factorial  design 
has  been  found  to  be  most  efficient  in  analysing  these 
effects  (7:189).  All  combinations  of  the  factors  involved 
are  tested  at  each  of  their  possible?  values  and  the  effect 
on  the  overall  result  is  measured  in  each  case. 

Deter mi nat i on  of  Factors .  Before  the  design  can  be 
run,  the  factors  to  be  tested  must  be  determined.  Any 
factor  (variable,  input  value)  that  can  change  the  output 


value  should  be  considered. 


For  the  SRT  mission.  all  the 


variables  were  reduced  to  -four  factors  that  can  affect  the 


desired  result:  grid  size,  allotted  search  time,  time  unit 
of  the  stochastic  process,  and  the  conditional  detection 
probabi 1 i ty . 

Grid  size  is  a  combination  of  the  search  radius  of  the 
aircraft  and  the  altitude  flown.  The  search  radius  is  how 
far  the  aircraft  sensors  can  "see"  ahead  of  the  aircraft 
(anywhere  from  90°  to  the  left  through  90°  to  the  right). 
Changing  altitude  changes  the  distance  the  aircraft  can  see. 
To  ensure  the  aircraft  can  scan  an  entire  grid  as  it  flies 
straight  through  the  grid's  center,  a  grid  size  twice  that 
of  the  search  radius  is  used.  Therefore,  from  the  center- 
line  of  the  grid,  looking  left  or  right  will  enable  the 
sensors  to  see  all  the  area  included  in  the  grid.  Because 
grid  size  and  search  radius  are  proportional  to  one  another, 
only  grid  size  will  be  considered  as  a  factor. 

Imbedded  in  the  second  factor,  allotted  search  time,  is 
the  speed  of  the  aircraft.  Search  time  is  how  much  time  the 
aircraft  has  to  search  for  the  SRT.  Faster  flying  aircraft 
can  search  more  area  in  the  allotted  time  than  slower  air¬ 
craft.  Changes  in  the  model  s  result  due  to  i ncreased /de¬ 
creased  search  time  can  indicate  an  effect  due  to  a  change 
in  the  aircraft's  speed  or  a  change  in  search  time.  Again 
both  factors  are  proportional  and  can  be  tested  for  signifi¬ 
cance  by  considering  just  the  search  time. 
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A  third  -factor  is  the  time  unit  of  the  stochastic 
process  (time  between  changes  in  the  SRT ' s  location).  As  in 
the  other  two  factors,  the  time  unit  is  determined  by  a 
combination  of  other  variables  such  as  dwell  time,  speed  of 
the  SRT,  and  the  average  distance  per  move  by  the  SRT. 
Changing  any  one  or  a  combination  of  these  variables  changes 
the  time  between  new  locations  for  the  SRT.  For  instance, 
since  the  grid  size  is  constant,  a  faster  speed  by  the  SRT 
means  less  time  is  spent  in  one  grid  whereas  a  longer  dwell 
time  means  a  longer  time.  Because  these  variables  directly 
affect  the  time  unit  used,  only  the  time  unit  needs  to  be 
considered  over  its  possible  range  of  values.  Changes  in 
the  time  unit  can  be  made  by  changing  any  one  of  the  imbed¬ 
ded  variables  with  the  subsequent  analysis  indicating 
whether  the  factor  caused  a  significant  change  in  the  final 
r esul t . 

Finally,  the  fourth  factor  is  the  conditonal  detection 
probability.  As  referred  to  here,  the  detection  probability 
is  the  probability  the  aircraft  detects  the  target  given  the 
sensors  are  looking  directly  at  the  SRT  and  is  computed  as 
the  product  of  two  other  factors.  Multiplying  the  target 
detection  factor  by  the  aircraft  sensor  system  detection 
factor  results  in  the  conditional  detection  probability 
(14:slide  10).  Therefore,  this  probability  is  computed  and 
input  into  the  model  as  a  constant  to  be  used  for  all  runs. 
So  while  it  is  a  probability,  it  remains  a  constant  factor 
for  one  set  of  input  data.  Varying  the  conditional  detec- 
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tion  probability  over  its  possible  range  will  indicate  its 
significance  as  well  as  the  significance  of  the  other  two 
given  factors. 

Screeni ng  Runs .  Usually  initial  screening  runs  in  a 
design  will  indicate  which  factors  have  a  significant  effect 
on  the  results  and  which  factors  can  be  disregarded  due  to 
their  lack  of  si gni f i cance.  However,  the  structure  of  the 
SRT  problem  is  slightly  different.  Rather  than  determine 
which  factors  or  interactions  are  significant  at  one  update 
time,  a  relationship  between  each  update  time  is  desired  to 
determine  if  any  factors  or  interactions  change  in  signifi¬ 
cance  as  the  update  time  changes. 

Three  update  times  were  chosen.  Since  significant 
changes  in  the  overall  detection  probability  will  more  than 
likely  occur  for  shorter  update  times  than  for  longer  update 
times,  SO  minutes,  160  minutes,  and  240  minutes  were  chosen 
for  the  screening  runs.  These  times  should  cover  the  range 
of  time  where  significant  changes  in  the  results  should 
occur . 

To  complete  the  screening  runs,  a  2*  factorial  design 
with  three  replicates  at  each  treatment  level  is  accom¬ 
plished.  Each  of  the  four  factors  is  considered  at  two 
levels  as  defined  in  Table  6.1.  Sixteen  total  runs  is 
required  to  complete  one  replicate  of  the  design.  For  this 
model,  a  complete  design  was  accomplished  three  times  for 
each  of  the  three  updates  resulting  in  a  total  of  48  runs  at 


78 


TABLE  6. 1 


List  of  Factors 

Letter _ Factor _ Low _ Hi  gh _ 

A  Grid  Size  17.4  miles  38.8  miles 

B  Search  Time  21  min  46  min 

C  Time  Unit  34.8  min  276.24  min 

_ D _ Detection  Prob. _ 0. 5 _ 1  ■  0 _ 

each  update  time.  But  to  determine  the  significance  of 

these  factors  on  the  probability  of  detecting  the  SRT ,  a 

model  to  estimate  the  probability  must  also  be  developed. 

To  estimate  the  probability  of  the  aircraft  finding  an 
SRT,  the  percent  of  successes  for  all  missions  flown  is 
used.  Estimating  the  probability  of  detection  is  a  binomial 
experiment.  A  binom-rl  experiment  is  a  sequence  of  identi¬ 
cal  and  independent  runs  which  can  result  in  only  one  of  two 
possible  outcomes  (6:78).  Each  mission  flown  results  in  a 
hit  (1)  or  a  miss  (0).  The  target  parameter,  p,  can  be 
estimated  by  Y/n  where  Y  is  the  number  of  hits  and  n  is  the 
total  number  of  missions  flown  at  one  update  time  (6:297). 
For  each  run,  the  number  of  hits,  Y,  is  computed  but  the 
value  of  n  is  still  needed  to  ensure  the  estimate  is  within 
a  specified  accuracy. 

By  specifying  a  maximum  variance  and  significance 
level,  the  value  of  n  can  be  computed.  To  be  957.  sure  that 
the  estimate  computed  is  within  .1  of  the  actual  probabil¬ 
ity,  957.  of  the  estimates  must  lie  within  two  standard 
deviations  of  the  actual  probability  with  repeated  computer 
runs  (6:317).  In  other  words 


TABLE  6.2 


Treatment  Levels  and  Results 


Probability 

of  Detection 

Est l mate 

BO  min  update 

160  min  u 

pdate 

240 

min  update 

Treatment  1 

2 

1 

o 

T 

■_> 

1 

<7 

3 

(  1 ) 

.  30 

.  30 

.  19 

.21 

.24 

.  19 

.  10 

.  19 

a  25 

a 

nn 

.  31 

.  20 

*7 

.  12 

.  13 

.  21 

.  13 

.  07 

b 

.41 

.  49 

.  36 

.31 

.41 

.39 

.31 

~T  *7 
a  JJ 

.37 

ab 

.44 

.  34 

.  31 

.28 

.  26 

.  24 

.  26 

a  25 

a  25 

c 

.  08 

.  06 

.07 

.  27 

.  28 

.21 

.32 

.  38 

.  36 

ac 

.  01 

.  00 

.  00 

.20 

.  16 

.  1  1 

.  20 

.  23 

.  15 

be 

.  10 

.  06 

.  10 

.  28 

» -1*6 

.  28 

.38 

.  44 

.42 

abc 

.  09 

.  07 

.  08 

.  29 

.  27 

.  27 

-T  r-V 

■ 

.31 

a  35 

d 

•  52 

.  51 

.  47 

.  42 

.  48 

.  46 

T  T 

a 

.  37 

.  43 

ad 

.  36 

.  38 

.38 

.  30 

.  37 

a  33 

.  24 

bd 

.77 

.  76 

.  73 

.  59 

.69 

.75 

.  45 

.52 

.57 

abd 

.  56 

.  63 

.58 

.  44 

.  40 

.57 

.  47 

.47 

.  47 

cd 

.  1  1 

.  08 

.  1 1 

.42 

.  35 

.28 

.  43 

.  48 

.41 

acd 

.  0  1 

.  00 

.  00 

.  37 

.32 

.  26 

.  29 

.  39 

.  36 

bed 

.  1  1 

.  08 

.  1 1 

.42 

•  35 

.28 

.  43 

.48 

.41 

abed 

.  1  1 

.  08 

.  1 1 

.42 

•  35 

.28 

.  43 

.  48 

.41 

The 

al  gebr ai c 

signs 

used  to 

calculate 

the 

effects 

of 

the  four  factors  is  displayed  in  Table  6.3  with  the  respec¬ 
tive  ANOVA  Tables  for  the  three  update  times  shown  in  Tables 
6.4,  6.5,  and  6.6.  While  definite  conclusions  may  be  hard  to 


draw  from  the  results,  some  general  inferences  can  be  made 


TABLE  6.3 


TABLE  6.4 


ANOVA  Table  -for  the  80  Minute  Update 


Source  of  Variation 


Grid  Size  <A) 

.  054 

Search  Time  (B) 

.  1645 

Time  Unit  (C) 

1 . 6465 

Detect  Prob  (D) 

.  1838 

AB 

.  00 1 3 

AC 

.  0073 

AD 

.  0105 

BC 

.  058 1 

BD 

.  0073 

CD 

.  1398 

ABC 

.  0099 

ABD 

.  0002 

ACD 

.  0068 

BCD 

.  0068 

ABCD 

.  0003 

Error 

.0413 

MS 


.  054 
.  1645 
1 . 6465 


87  47 


41.84 
127.46 
1275. 74 


1.012 


8.  138 


5.  62to 


8. 


7.677 


TABLE  6 


ANOVA  Table 

■for  the 

160  Minute  Update 

Source  of  Variation 

SS 

df 

MS 

F° 

Gr  iu  Size  (A) 

«“*n  -7  / 

•  t  O 

\ 

X 

r»0“ 7  L. 

•  U  f  W 

27 . 324* 

Search  Time  (B) 

.  1376 

1 

.  1376 

42. 94- 

Time  Unit  (C) 

.  0554 

1 

.  0554 

17.265- 

Detect  Prob  (D) 

.  3024 

1 

.  3024 

94.383- 

AB 

. 00002 

1 

. 00002 

.  0065 

AC 

.  0248 

1 

.  0248 

7.725- 

AD 

.  0020 

1 

.  0020 

.625 

BC 

.  0369 

1 

.  0369 

1 1 . 50- 

BD 

•  000004 

1 

. 000004 

. 00128 

CD 

.  0514 

1 

.  05 1 4 

16.03- 

ABC 

.  0083 

1 

.  0083 

2.58 

ABD 

.  0002 

1 

.  0002 

.  052 

ACD 

.  0158 

1 

.  0 1 58 

4 . 92 lto 

BCD 

.0144 

1 

.0144 

4. 43* 

ABCD 

Err  or 

.  0003 

.  1 025 

1 

.  0003 

.  0032 

.  079 

Total 

.  8395 

47 

TABLE  6.6 
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ANQVA  Table  -for  the  240  Minute  Update 


Source  of  Variation 


bria  6ize  u-w 

.  u  /  ZV 

Search  Time  (B) 

.  1530 

Time  Unit  (C) 

.  0336 

Detect  Prob  (D) 

.  2255 

AB 

.  0073 

AC 

.  0006 

AD 

.  0029 

BC 

.  0204 

BD 

.  0002 

CD 

.  0213 

ABC 

.  0035 

ABD 

.  0032 

ACD 

.  0063 

BCD 

.  0054 

ABCD 

.  00 1 5 

Err  or 

.  0633 

T  otal 

.  6207 

*  Significant  at  1  percent 
b  Significant  at  10  percent 


77 . 35 
16.98 
113. 997- 
3. 668b 
.  3065 
1 . 444 


10.75- 
1 . 767 
1 . 599 
.  183^ 
.736 
.  772 


Analysis  of  Screen  1 ng  Runs .  For  each  of  the  three 
update  times  selected,  all  four  main  effects  (effects  due 
only  to  one  factor)  caused  significant  changes  in  the  result 
(7:189).  Using  the  magnitude  of  the  resulting  F  test  as  an 
indication  of  significance,  the  order  of  the  four  factors 
from  most  significant  to  least  significant  is  the  condi¬ 
tional  detection  probability,  allotted  search  time,  grid 
size,  and  the  time  unit  of  the  stochastic  process.  Although 
all  four  factors  are  significant,  if  tradeoffs  are  required, 
a  decision  maker  (nay  be  able  to  use  the  general  inferences 
just  stated  to  make  a  future  decision.  Even  more  helpful 
however ,  may  be  the  inferences  that  can  be  made  from  the 
interaction  effects. 

Only  interaction  effects  involving  two  factors  indi¬ 
cated  a  significant  change  in  the  resulting  probability 
estimate.  Any  interaction  effect  involving  the  grid  size 
(A)  was  not  found  to  be  significant  On  the  other  hand,  the 
time  unit  of  the  process  (C>  was  found  to  be  significant 
when  combined  with  either  the  allotted  search  time  (B)  or 
the  conditional  detection  probability  (D).  Since  B  and  D 
individually  seemed  to  be  the  most  significant,  the  combined 
significance  was  probably  due  more  to  just  the  effects  of  B 
and  D  wi th  C  s  effects  relatively  insignificant. 

Results  of  thus  design  do  not  indicate  which  factors 
are  necessarily  more  significant  than  the  others  but  the 
results  can  indicate  which  factors  may  need  emphasis  if 


tr adeof  f  s 


must  be  made. 
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are  two  -fold.  First,  all  main  effects  are  significant. 
Second,  each  main  effect  is  significant  at  all  update  times 

\ 

1  in  the  range  tested.  If  a  decision  maker  has  to  decide 

I 

'  between  a  change  in  two  factors  such  as  aircraft  altitude 

I 

I  <grid  size)  and  sensors  to  use  for  the  search  (conditional 

detection  probability),  results  from  this  design  can  steer 
the  decision  maker  in  the  right  direction.  Knowledge  of  the 
I  significant  factors  and  interaction  effects  allows  an  input 

data  set  to  be  built  for  use  in  the  "production"  runs  to 
demonstrate  the  usefulness  of  the  model. 

t 

I 

Product  1  on  Runs 

To  demonstrate  the  final  output,  an  example  SRT  mission 
is  simulated.  All  variables  are  set  to  realistic  values  to 
simulate  an  actual  mission.  The  values  remain  the  same  for 
all  computer  runs.  Initially,  100  missions  are  simulated 
with  the  aircrew  receiving  an  intelligence  update  30  minutes 
prior  to  entering  the  search  area.  The  number  of  hits 
(divided  by  100)  is  used  as  the  estimate  for  the  probability 
of  detection  for  that  update  time.  A  plot  is  started  using 
this  result  as  the  first  data  point. 

Continuing,  another  100  missions  are  flown  with  a  60 
minute  update  and  the  process  repeated  as  before.  Each 
subsequent  run  has  the  aircrews  receiving  an  update  30 
minutes  earlier  than  the  previous  set  of  aircrews.  After 
the  16fc^  run  representing  an  eight  hour  update,  a  one  hour 


increment  is  used  between  runs  for  the  final  eight  runs. 


A 


total  of  24  runs  are  made  covering  a  range  of  update  times 
from  30  minutes  to  960  minutes  (16  hours).  The  table  of 
results  is  shown  in  Table  6.7  with  the  plot  of  all  update 
times  and  the  respective  detection  probabilities  shown  in 
Fig.  6.1.  At  this  point  a  decision  maker  would  have  some 
valuable  information  to  help  utilize  his  intelligence  re¬ 
sources.  However,  since  this  plot  is  only  good  for  this  set 
of  input  data,  a  new  plot  has  to  be  generated  for  each 
different  mission. 


TABLE  6.7 

Table  of  Example  Results 


Run 

Update  (nun) 

Detect  Frob 

1 

30 

.  92 

o 

60 

.  90 

•.!> 

90 

.  46 

4 

120 

•  52 

5 

150 

.  58 

6 

180 

.47 

7 

210 

.  44 

a 

240 

.  49 

9 

270 

.  49 

10 

300 

.  47 

1 1 

330 

.  42 

12 

360 

.  40 

13 

390 

.  40 

14 

420 

.  47 

15 

450 

.  43 

16 

480 

.  43 

17 

540 

•  36 

18 

600 

.  31 

19 

660 

20 

720 

.  27 

21 

780 

.  28 

nn 

840 

•  J  -J 

O  T 

900 

T  n 

24 

960 

.  36 

Developing  the  Confidence  Bands 


For  each  run,  the  percentage  of  hits  out  of  100  mis¬ 
sions  was  plotted  against  the  update  time  cor respondi ng  to 
the  run.  However,  as  a  further  step,  a  confidence  band  was 
drawn  around  the  estimated  probabilities  to  better  indicate 
the  desired  relationship. 

To  obtain  a  confidence  band  around  the  computed  curve, 
the  Bonferroni  approach  was  used.  A  confidence  interval  for 
each  estimate  is  constructed  and  each  point  plotted  on  the 
same  graph  with  the  estimated  pr obabi 1 i t i es .  By  connecting 
the  points  both  above  and  below  the  original  plot,  a  "confi¬ 
dence  band"  will  result.  However,  the  accuracy  of  the  band 
and  inferences  that  can  be  made  depend  on  the  confidence 
coefficient  chosen  and  the  number  of  confidence  intervals  in 
the  "family"  (how  many  estimates  are  involved). 

Each  confidence  interval  was  built  using  a  confidence 
coefficient  (  )  of  .01  which  means  997.  of  the  intervals 
computed  at  each  update  time  will  include  the  actual  parame¬ 
ter  value  (actual  detection  probability)  for  that  update 
time  (8:707).  Qveral 1 ,  since  each  run  is  independent,  the 
confidence  coefficient  is  (  1  -  )  «  where  <x.  is  the  confidence 
coefficient  for  one  interval  and  g  is  the  number  of  esti¬ 
mates  (8:150).  Therefore,  the  computed  confidence  coeffi¬ 
cient  is: 

(!-<*.)«  -  (l-.Ol)3*  ^  .  79  (20) 
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Bonf  erron  1 


confidence  intervals  are  constructed  using 


the  following  formula (8: 158) : 

+  Bs<?„)  (21) 

where:  yv,  is  the  estimated  probability 

B  is  the  Bonferroni  constant  (B=t ( l-o< /2g; n-2) ) 
g  is  the  number  of  confidence  intervals 
h  is  the  sample  size 

s(yv,)  is  the  standard  deviation  of  the  estimate 
(  <S  =  p  ( 1-p  >  /n  ) 

The  results  of  each  calculation  is  shown  in  Table  6.8  along 
with  the  previous  results  and  the  confidence  band  is  plotted 
on  the  graph  with  the  estimated  probabilities  in  Fig.  6.2. 
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TABLE  6.8 

Example  Problem  Results  With  Confidence 
Interval  Calculations 


B=3. 405 

Update  (min)  Detect  Prob  s(y^)  B(s(Vr,)) 

30  . 92  . 027  . 092 

60  .90  .030  .102 

90  .46  .050  .170 

1 20  . 52  . 050  .  1 70 

150  .58  .049  .167 

180  .47  .050  .170 

210  .44  .050  .170 

240  .49  .050  .170 

270  .49  .050  .170 

300  .47  .050  .170 

330  .42  .049  .167 

360  .40  .049  .167 

390  .40  .049  .167 

420  .47  .050  .170 

450  .43  .050  .170 

480  .43  .050  .170 

540  .36  .048  . 163 

600  .31  .046  .157 

660  .32  .047  .160 

720  .27  .044  .150 

780  .28  .045  .153 

840  .33  .047  .160 

900  .32  .047  .160 

960  .36  .048  .163 


Limits 
Lower  Upc 


828 

1.012 

798 

1 . 002 

290 

.  630 

350 

.  690 

413 

.747 

300 

.  640 

270 

.610 

320 

.  660 

320 

.  660 

300 

.  640 

253 

.587 

— i77 

ZOO 

.  567 

nr77 
Z.  o 

.567 

300 

.  640 

260 

.  600 

260 

.  600 

197 

«=:'“?  t 

153 

.  467 

1 60 
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Before  the  project  was  undertaken,  a  generally  accepted 
thesis  was  that  the  closer  the  aircraft  is  to  the  target 
when  the  update  is  received,  the  better  the  chance  o-f  find¬ 
ing  the  target.  Results  from  the  example  run  are  in  line 
with  this  thinking.  However,  as  questioned  earlier,  how 
"close"  should  the  update  be  to  the  aircraft  arrival  time 
to  produce  an  acceptable  detection  probability. 

Updates  at  30  minutes  and  60  minutes  resulted  in  much 
higher  detection  probabilities  (by  more  than  .30)  than  any 
other  point  on  the  plot.  As  it  turns  out  the  input  data 
used  required  69.6  minutes  for  the  SRT  to  move  half  way 
through  one  grid.  Also,  besides  the  move  time,  90  minutes 
of  dwell  time  was  required.  Therefore,  an  SRT  spent  159.8 
minutes  in  grids  that  allowed  stops.  Since  the  aircraft 
does  not  know  how  long  the  SRT  has  been  in  the  grid  location 
reported  at  the  update  time,  the  assumption  made  earlier  in 
Chapter  4  is  that  half  of  the  total  possible  time  has  been 
spent  in  the  reported  location.  Based  on  this  assumption, 
79.8  minutes  was  used  in  every  aircraft  mission  throughout 
the  example.  Both  the  30  and  60  minute  updates  are  less 
than  the  assumed  time  which  means  the  aircraft  assumes  the 
target  cannot  move  before  the  search  begins.  When  a  90 
minute  update  (or  earlier)  is  received,  one  or  more  moves 
will  be  assumed  by  the  aircraft  and  the  results  show  a 
significant  drop  in  the  probability  of  detection  for  that 


particular  update  time.  Something  significant  happens 

between  the  60  and  90  minute  updates. 

From  this  example,  the  conclusion  reached  is  that  the 
best  detection  probabilities  can  be  achieved  i  -f  the  air¬ 
craft  receives  an  update  at  a  time  that  is  less  than  half 
the  assumed  dwell  time  of  the  target  (where  the  dwell  time 
must  include  the  move  time  through  a  grid,  set  up  time,  tear 
down  time,  and  actual  sitting  time  in  one  spot).  Actual 
values  from  this  example  data  may  not  be  accurate,  but 
significant  changes  in  the  detection  probabi 1 i t 1 es  indicate 
the  value  of  a  possible  rule  of  thumb  to  carry  over  to 
future  planning. 

Conclusi on 

From  the  design  of  experiments,  inferences  were  made 
concerning  the  significance  of  the  selected  factors.  All 
main  effects  were  found  to  be  significant  at  all  update 
times.  Search  time  (B)  and  the  conditional  detection  proba¬ 
bility  (D)  were  thought  to  affect  the  overall  probability  of 
detection  slightly  more  because  of  the  magnitude  of  the  F 
test  values.  Interactions  involving  two  factors  seem  to  lend 
credence  to  this  inference.  When  the  time  unit  (C)  was 
combined  with  B  or  B,  the  interaction  effect  was  also  sig¬ 
nificant  for  all  update  times.  With  these  inferences  in 
mind,  simulation  runs  were  made  to  develop  the  desired  plot. 

Aircraft  missions  were  simulated  at  different  update 
times  to  estimate  the  probability  of  detection  at  each 
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selected  time.  From  the  estimate,  a  plot  was  developed  from 
which  further  inferences  could  be  made  concerning  the  "best" 
time  to  send  an  intelligence  update  to  the  attacking  air¬ 
crew.  A  general  conclusion  drawn  was  that  the  update  time 
should  be  planned  for  a  time  less  than  one  half  the  total 
time  an  SRT  spends  in  one  grid.  While  this  is  a  general 
result  from  one  set  of  sample  data,  further  testing  could 
prove  this  inference  to  be  a  valuable  rule  of  thumb. 


VII.  Recommend at i ons  fo r  Further  Stud 


Introduction 

Computer  simulations  are  only  as  good  as  the  assump¬ 
tions  they  are  built  upon  and  the  information  put  into  the 
model.  Building  a  model  to  simulate  SR'T  missions  has  so 
much  uncertainty  surrounding  the  target  movement  and  air¬ 
craft  search,  that  the  assumptions  and  input  data  are  criti¬ 
cal  in  determining  the  value  of  the  results.  More  specifi¬ 
cally,  two  areas  require  more  research  to  improve  the 
model's  accuracy. 

Ai r cr af t  Sear ch 

To  form  a  realistic  search  pattern,  the  grids  searched 
by  the  aircraft  must  all  be  connected.  In  this  model,  the 
search  pattern  that  was  developed  borrowed  some  ideas  from 
articles  written  about  search  theory.  However,  the  frame  of 
reference  in  the  literature  did  not  consider  feasible  air¬ 
craft  flight  paths.  As  stated  in  Chapter  1,  most  research 
addressed  the  problem  of  maximising  the  probability  of  de¬ 
tection  for  a  given  amount  of  search  effort  or  time. 
Searching  in  this  manner  may  result  in  the  highest  possible 
detection  probability  (or  at  least  an  upper  bound  on  the 
probability)  but  does  not  ensure  a  feasibe  search  path. 

More  research  is  required  to  ensure  the  grids  generated 
+ or  the  aircraft  to  search  form  a  feasible  search  "path". 


T  wo 


approaches  came  to  mind  while  studying  the  problem  but 


time  did  not  allow  -further  implementation.  One  approach 
could  involve  restrictions  placed  on  the  one-step  transition 
matrix.  If  the  transition  pr obab i 1 1 t i es  were  input  such 
that  the  highest  probabi 1 1 t 1 es  centered  on  a  certain  "path" 
of  grids  (such  as  fallowing  a  road) ,  the  resulting  uncondi¬ 
tional  pr obab 1 1 1 t 1 es  would  reflect  an  appropriate  path  (Fig. 
7.1).  However,  one  must  be  careful  in  this  situation  so  as 
not  to  restrict  the  search  too  much  so  that  it  is 
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Fig.  7.1.  Sample  Search  Area  Where  the 
F'r obab  i  1  i  t  i  es  Form  a  F'ath 

A  second  approach  takes  the  current  search  pattern  one 
step  farther.  After  computing  the  search  pattern  that  gives 
the  upper  bound  on  the  probability,  a  check  should  be  done 
on  the  grids  to  be  searched  to  ensure  a  path  is  formed.  If 
not,  new  grids  can  be  selected  to  make  the  pattern  a  path  by 
continuing  down  the  sorted  unconditional  probability  list. 
In  thus  way,  the  new  grids  selected  still  have  the  highest 
probability  of  containing  the  target  but  with  the  added 
constraint  of  having  to  form  a  path  with  the  other  selected 


gr  ids.  A  check  such  as  this  could  be  incorporated  very 


e  a  s  l  1  y 


into  the  present  model  as 


FORTRAN  subroutine. 


While  the  actual  values  tor  the  probability  of  detection 
estimates  may  change,  the  shape  of  the  curve  generated  from 
the  estimates  versus  the  update  times  would  not  degrade  much 
and  would  more  realistically  reflect  a  real  aircraft  mission 
profile. 

T arget  Time  i n  Current  Locati on 

Another  area  that  could  improve  the  model  deals  with 
how  long  the  target  has  been  in  the  current  location  when  an 
update  is  received  by  the  aircrew.  An  assumption  made  in 
this  model  is  that  when  the  target  is  spotted  by  intelli¬ 
gence  or  by  whatever  means,  the  SRT  has  been  there  half  of 
the  total  time  normally  spent  in  one  location.  By  dividing 
the  time  in  half,  the  error  from  not  knowing  the  exact  time 
should  be  reduced.  For  targets  with  short  dwell  times  the 
error  will  be  small.  However,  as  dwell  times  get  larger, 
errors  also  increase.  Since  another  assumption  for  the  model 
is  that  only  short  dwell  time  targets  are  being  considered, 
more  flexibility  can  be  built  into  the  model  by  developing  a 
routine  to  better  estimate  how  long  the  target  has  been  in 
the  current  location  when  the  update  is  observed. 

As  an  example,  79.8  minutes  was  half  the  normal  time 
spent  in  one  grid  for  the  SRT  s  in  the  previous  example. 
Actual  time  already  spent  in  the  reported  grid  at  the  update 
time  averaged  about  62  minutes  as  calculated  from  the  sample 
output.  In  fact,  only  one  run  had  an  average  time  greater 
Therefore,  the  target,  on  the  average,  had  not 


than  79.8. 


been  in  the  reported  location  as  long  as  the  aircrew  as¬ 
sumed.  Because  of  this,  the  aircrew  predicted  a  move  by  the 
target  before  the  target  actually  moved. 

On  the  average,  the  searches  were  taking  place  in  grids 
too  -far  away  -from  the  reported  position  given  at  the  update 
time.  A  routine  that  could  better  estimate  how  long  a 
target  has  been  in  its  current  location  could  improve  the 
accuracy  of  the  model.  However,  more  intelligence  updates 
may  be  required  to  compute  this  estimate  and  it  may  prove  to 
be  too  much  effort  for  the  benefit  gained.  On  the  other 
hand,  the  extra  research  might  turn  up  an  estimating  proce¬ 
dure  that  could  improve  the  accuracy  of  the  model. 

Conclusion 

Strategic  relocatable  targets  are  very  uncertain  crea¬ 
tures.  Much  uncertainty  surrounds  their  pattern  of  movement 
and  attempts  to  find  them.  HQ  Strategic  Air  Command  and 
Pentagon  planners  are  continually  working  to  decrease  this 
uncertainty  and  to  plan  strikes  against  SRT ' s  with  a  high 
probability  of  success.  However,  much  more  work,  is  still 
needed . 

From  this  model  some  of  the  uncertainty  about  the 
target  movement  has  been  put  to  rest.  So  many  problems  still 
exist  because  every  target  in  every  section  of  enemy  terri¬ 
tory  poses  a  unique  threat  with  a  unique  set  of  information 
required  to  find  and  destroy  each  SRT.  To  build  a  general 
model  to  cover  all  contingencies  based  on  the  appropriate 


FORTRAN  Program 


PROGRAM  MAIN 

************************************************************ 


*  * 

*  THIS  PROGRAM  WORKS  IN  CONJUNCTION  WITH  THE  SLAM  PROGRAM  * 

*  "RELOC. TARGET"  TO  SIMULATE  STRATEGIC  BOMBERS  ATTACKING  * 

*  STRATEGIC  RELOCATABLE  TARGET  ( SRT  S )  .  IT  IS  USED  TO  DO  * 

*  THE  OPERATIONS  THAT  SLAM  CANNOT  DO.  FOR  THIS  TYPE  OF  * 

*  PROGRAM,  THE  PRIMARY  OPERATIONS  ARE  MATRIX  OPERATIONS  * 

*  AND  SORTING.  ALSO,  SOME  SPECIAL  OUTPUT  INFORMATION  IS  * 

*  EXTRACTED  AND  PRINTED  TO  AN  EXTERNAL  FILE  TO  AID  IN  * 

*  FOLLOW  ON  ANALYSIS.  ALL  TIMES  ARE  IN  MINUTES  AND  RATES  * 

*  ARE  PER  MINUTE.  * 

*  * 


************************************************************ 

♦DECLARATIONS 

D I MENS I ON  NSET ( 1 OOOO ) 

CHARACTER*?  TYPE 

REAL  SPEED, AVGD IS, TGTDET 

I NTEGER  DWELL , SETUP , TERDWN , SEARCH , VELCT Y 
I  NTEGER  I NTTME  ,  MPT  I  ME  ,  CRWRST  ,  RDUN I T  ,  WTUN I T  ,  MAX  BOX 
INTEGER  RUNCTR,UPINCR 
♦COMMON  STATEMENTS 

COMMON/SCOM 1 / ATRI B  < 100)  ,DD (100)  , DDL (100)  ,DTNOW, 1 1 , MFA , 

MSTQF’ ,  NCLNR 

1  ,  NCRDR  ,  NF'RNT  ,  NNRUN  ,  NNSET  ,  NT  APE  ,  SS  (  1 00 )  ,  SSL  (100)  ,  TNEXT  , 

TNOW, XX (100) 

* 

COMMON/  I NF'TS/T  YPE ,  SPEED ,  DWELL ,  SETUP  ,  TERDWN  ,  SEARCH , 

VELCTY 

1 , AVGD I S , I NTTME , MPT IME , CRWRST , TGTDET , SERRAD , LSTGRD , 

RUNCTR , UP I NCR 

* 

COMMON/ ARAYS/ONESTF'  (  169,169)  , STATUS  (0:  169)  ,  ACFTDT  (169) 
1  ,SRTPRB(  169,2)  , DETECT  (169)  ,  ALLSTF'  (169,169) 

1  ,  POST  IN  ( 169)  ,  TMF'MAT  (169,169)  ,TMPLOC  (169) 

1  ,  TMPSRT  (169)  ,  TF'ST AT  (0:169) 

* 

COMMON/RDWRT/RDUNIT, WTUN IT 

* 

COMMON/ ARYSZE/MAXBOX 

* 

COMMON  QSET ( 1 0000 ) 

EC!U  I  VALENCE  (NSET  ( 1 )  ,QSET(1)  ) 

♦OF'EN  INPUT  AND  OUTPUT  FILES 

OPEN (UNI T=20, FI LE= 'RELOC. SEVENV ' ,STATUS=  OLD ' ) 

OPEN (UNI T=21 ,FILE= ' RELOC . RESLTV ' ,STATUS=  NEW ' ) 
♦INITIALIZE  SLAM  VARIABLES 
NNSET = 10000 


iv  v"  nv1*  v 


NCRDR=5 
NPRNT=6 
NTAPE=7 
NPLOT =2 

♦SET  INPUT  AND  OUTPUT  FILE  UNIT  NUMBERS 
RDUNIT=20 
WTUNIT=21 

♦SET  SIZE  OF  SEARCH  AREA  BOX  FOR  THIS  RUN 
MAXB0X=49 
♦READ  IN  DATA 

REW I ND ( UN  I T =RDUN I T ) 

♦QNESTP ( I  ,  J  )  IS  THE  ONE  STEP  TRANSITION  MATRIX 
DO  10,  1=1, MAXBOX 

READ  < RDUN I T  ,  ♦ )  <ONESTP(I  ,  J)  ,  J  =  1  , MAXBOX) 

10  CONTINUE 

♦STATUS(I)  IS  THE  VECTOR  CONTAINING  THE  STATUS  OF  EACH 
♦GRID  (0,1, OR  2) 

READ (RDUN IT ,♦) (TPSTAT(I) ,1=0, MAXBOX) 

♦ACFTDT  < I )  IS  THE  VECTOR  OF  AIRCRAFT  SENSOR  DETECTION 
♦FACTORS  PER  GRID 

READ  < RDUN IT,*)  (ACFTDT ( I ) , 1  =  1 , MAXBOX) 

♦SET  TYPE  OF  TARGET 
TYPE=' MOVABLE' 

♦SET  SPEED  OF  TARGET 
SPEED=0. 25 

♦SET  DWELL  TIME  OF  TARGET  (O  FOR  MOBILE  TARGETS) 

DW£LL=60 

♦SET  TIME  REQUIRED  FOR  SETTING  UP  IN  A  NEW  LOCATION 
♦AFTER  A  MOVE 
SETUP= 15 

♦SET  TIME  REQUIRED  TO  DISMANTLE  EQUIPMENT  FOR  A  NEW  MOVE 
TERDWN=15 

♦SET  SEARCH  TIME  AIRCRAFT  IS  ALLOWED 
SEARCH=21 

♦SET  AIRSPEED  OF  AIRCRAFT 
VELCTY=6 

♦SET  AVERAGE  DISTANCE  TARGET  MOVES  EACH  TIME  IT  RELOCATES 
AVGD I S= 17.4 

♦SET  INTELLIGENCE  CYCLE  TIME 
INTTME= 120 

♦SET  TIME  ALLOWED  FOR  MISSION  PLANNING 
MPT IME=60 

♦SET  TIME  REQUIRED  FOR  CREW  REST 
CRWRST =480 

♦SET  THE  TARGET  DETECTION  FACTOR 
TGTDET  =1.0 

♦INITIALIZE  THE  RUN  COUNTER 
RUNCTR= 10 

♦INITIALIZE  THE  INCREMENT  OF  TIME  FOR  CHANGING  THE 
♦UPDATE  TIME 

UP  I NCR=60 

♦RUN  THE  SLAM  PROGRAM 
CALL  SLAM 


END 

****************** ***** ****** ******************************* 
********* *********** *********************** *********** ****** 
SUBROUTINE  INTLC 

************************************************************ 
♦SUBROUTINE  INTLC  IS  CALLED  AUTOMATICALLY  BEFORE  EACH  RUN  TO 
♦INITIALIZE  APPROPRIATE  VARIABLES  TO  SPECIFIED  VALUES 
************************************************************ 
♦DECLARATIONS 

CHARACTER*?  TYPE 
REAL  SERRAD, TBTDET, SPEED, AVGDIS 
INTEGER  I NTTME , MPT I ME , CRWRST , DWELL 
I NTEGER  SETUP , TERDWN , VELCTY , SEARCH 
INTEGER  MAXBOX , RDUN I T , WTUN I T , LSTGRD 
INTEGER  RUNCTR , UP  I NCR 
♦COMMON  STATEMENTS 

COMMON/SCOM1 / ATRIB ( ICO)  , DD ( 100)  , DDL (100)  ,DTNOW, 1 1 , MFA , 

MSTOP , NCLNR 

1 , NCRDR , NPRNT , NNRUN , NNSET , NT APE , SS ( 1 00 ) , SSL ( 1 00 ) , TNEXT , 

TNOW,XX (100) 

* 

COMMON/  I NPTS/TYF'E ,  SPEED  ,  DWELL  ,  SETUP  ,  TERDWN ,  SEARCH , 

VELCTY 

1  ,  AVGD I S  ,  INTTME ,  MPT  IME  , CRWRST  ,  TGTDET  ,  SERRAD ,  LSTGRD  , 

RUNCTR, UP  I NCR 

* 

COMMON/ ARAYS/ONESTP ( 169, 169) , STATUS (0: 169) ,ACFTDT(169) 
1 ,SRTPRB( 169,2) ,DETECT(169> , ALLSTP ( 1 69 , 1 69 ) 

1  , POSTIN  (169)  ,TMPMAT  (  169, 169)  ,TMF'L0C(169) 

1  , TMF'SRT  (169)  , TPSTAT  (0:169) 

* 

COMMON /RDWRT /RDUNI T , WTUN IT 

* 

COMMON/ ARYSZE/ MAX BOX 
♦SET  THE  SEARCH  RADIUS  OF  THE  AIRCRAFT 
5ERRAD=8. 7 

♦ESTABLISH  MINIMUM  FUEL  FOR  THE  AIRCRAFT  TO  STOP  SEARCH 
X  X ( 7 ) =50 . 0 

*SET  FLIGHT  TIME  FROM  TAKEOFF  TO  BEGIN  SEARCH 
XX (3) =360. 0 

♦SET  FUEL  LOAD  FOR  THE  AIRCRAFT  (IN  1000  LBS) 

XX (6) =400. 0 

♦SET  MAXIMUM  NUMBER  OF  AIRCRAFT  TO  FLY  PER  RUN 
X  X ( 1 4 ) = 1 00 . O 

♦SET  TIME  BETWEEN  FUEL  DECREMENTS 
XX (36) =30. 0 

♦SET  AMOUNT  OF  FUEL  TO  DECREASE  EACH  DECREMENT  (IN  1000  LBS) 
XX (37) =0. 0 

♦COMPUTE  SYSTEM  CYCLE  TIME 

XX (31 ) =INTTME+MPTIME+CRWRST+XX (3) 

♦COMPUTE  GRID  SIZE 

XX ( 39 ) =SERRAD*2 . 0 


IN  ONE  GRID) 


♦COMPUTE  TIME  IT  TAKES  AIRCRAFT  TQ  SEARCH  ONE  GRID 
XX (5) =XX (39) /VELCTY 
♦SAVE  VALUE  QF  FUEL  LOAD 
XX  <40) =XX (6) 

♦COMPUTE  NUMBER  OF  GRIDS  AIRCRAFT  CAN  SEARCH  IN 
♦ALLOTTED  TIME 

I I=SEARCH/XX (5) 

XX (9) =1 1 

♦COMPUTE  THE  CONDITIONAL  DETECTION  PROBABILITY  (PROBABILITY 
♦AIRCRAFT  SENSES  THE  TARGET  GIVEN  THEY  ARE  BOTH  IN  THE  SAME 
♦GRID) 

DO  10,  1=1 , MAX BOX 

DETECT ( I ) =TGTDET ♦ACFTDT  ( I ) 

10  CONTINUE 

♦COMPUTE  GROUND  TIME  OF  AIRCRAFT  (TIME  BEFORE  TAKEOFF) 
XX(1)=XX(31)-XX(3) 

♦COMPUTE  NUMBER  OF  GRIDS  ON  EACH  SIDE  OF  SEARCH  AREA  BOX 
I  I=MAXBOX^O.  5 
XX (25) =1 1 

♦SET  TIME  OF  ARRIVAL  IN  FIRST  GRID  TO  ZERO 
XX (49) =0. 0 

♦COMPUTE  TOTAL  DWELL  TIME  (TOTAL  TIME  IN  ONE  GRID)  FOR 
♦THE  TARGET 

X  X ( 2 1 ) =DWELL+SETUF'+TERDWN 
♦INCREMENT  RUN  COUNTER 
RUNCTR=RUNCTR+ 1 
XX (47) =RUNCTR 

♦COMPUTE  THE  INTELLIGENCE  UPDATE  TIME  FOR  NEXT  RUN 
XX (2) =RUNCTR^UPINCR 
♦FOR  MOVABLE  TARGETS 

IF (TYPE. EQ.  MOVABLE' )  THEN 

♦  COMPUTE  TIME  TO  MAKE  ONE  MOVE  TO  A  NEW  LOCATION 
XX (38) =AVGDIS/SPEED 

♦  COMPUTE  TIME  TO  MAKE  ONE  HALF  MOVE  TO  A 

♦  NEW  LOCATION 
XX (20) =XX (38) / 2 . 0 

♦  COMPUTE  TIME  UNIT  OF  STOCHASTIC  PROCESS 

XX (41 ) =XX (21 ) +XX (38) 

END  IF 

♦FOR  MOBILE  TARGETS  (CONSTANTLY  IN  MOTION) 

IF (TYPE. EQ. 'MOBILE ' )  THEN 

♦  COMPUTE  TIME  TO  MAKE  ONE  MOVE  TO  A  NEW  LOCATION 
XX (38) =XX (39) /SPEED 

♦  COMPUTE  TIME  TO  MAKE  ONE  HALF  MOVE  TO  A 

♦  NEW  LOCATION 
XX (20) =XX (38) / 2 . 0 

♦  COMPUTE  TIME  UNIT  OF  STOCHASTIC  PROCESS 

XX (41 ) =XX (38) 

END  IF 

♦SAVE  INITIAL  STATUS  OF  ALL  GRIDS 
DO  20,  I=0,MAXB0X 

STATUS ( I ) =TPST  AT ( I ) 

20  CONTINUE 


IN  MOTION) 


NEW  LOCATION 


♦  SAVE 


RETURN 

END 

****** ************ ************** ******************* ********* 
************************************************************ 
SUBROUTINE  EVENT ( I > 

************************************************************ 
♦SUBROUTINE  EVENT  IS  CALLED  BY  EVENT  NODES  IN  THE  SLAN 
♦PROGRAM  TO  DO  MATRIX  OPERATIONS  AND  OUTPUT  ROUTINES  THAT 
♦SLAM  CANNOT  DO 

************************************************************ 

♦DECLARATIONS 

DIMENSION  GDSRCH< 169,2) 

REAL  TEMP  ,  TIME  , PRETK 

I NTEGER  NOMOVE , SRTCTR , RUNC TR , UP I NCR , GR I DUP 
INTEGER  ACFCRT , RDUNI T , WTUNIT , MAX BOX , LSTGRD 
♦COMMON  STATEMENTS 

COMMON/SCOM1 / ATRIB ( 100)  ,DD  < 100)  , DDL (100)  ,DTNOW, 1 1 , MFA , 

MSTOP , NCLNR 

1 , NCRDR , NPRNT , NNRUN , NNSET , NT APE , SS ( 1 00 ) , SSL ( 1 00 ) , TNEXT , 

TNOW , X  X ( 1 00 ) 

* 

COMMON/ INPTS/TYPE , SPEED , DWELL , SETUP , TERDWN , SEARCH , 

VELCTY 

1 , AVGD I S , I NTTME , MPT  I ME , CRWRST , TGTDET , SERRAD , LSTGRD , 

RUNCTR , UP I NCR 

* 

COMMON /AR A YS/ONESTP ( 169 ,169) , STATUS (0: 169) ,ACFTDT(169) 
1 ,SRTPRB( 169,2) ,DETECT(169) , ALLSTP ( 169 , 169) 

1 , POST  IN (169)  , TMPMAT ( 1 69 , 1 69 )  ,TMPL0C(169) 

1 , TMPSRT (169), TPSTAT (0:169) 

* 

COMMON / RDWRT / RDUN IT , WTUNIT 

* 

COMMON/ ARYSZE/MAXBOX 
♦BRANCH  TO  THE  CORRECT  EVENT 

GO  TO  ( 1 , 2 , 3 , 4 , 5 , 6 , 7 , 6 , 9 , 1 0 , 1 1 , 1 2 , 1 3  >  , I 
************************************************************ 
♦EVENT  1  COMPUTES  THE  TARGET  LOCATION  DISTRIBUTION  PREDICTED 
♦FOR  THE  AIRCRAFT  ARRIVAL  TIME  IF  THE  UPDATE  IS  PRIOR  TO 
♦AIRCRAFT  TAKEOFF  TIME.  THE  DISTRIBUTION  IS  SORTED  FROM 
♦HIGH  TO  LOW  PROBABILITIES  SO  THE  AIRCRAFT  SEARCHES  THE 
♦HIGHEST  PROBABILITY  AREAS  FIRST 
************************************************************ 

*  COMPUTE  THE  TIME  IT  TAKES  THE  TARGET  TO  MOVE  OUT  OF 

*  THE  FIRST  GRID 

1  TEMF'=  (XX  (21)  /2.0)+XX  (20) 

IF (TEMP. LE. XX (2) )  THEN 

*  TGT  MOVES  OUT  OF  FIRST  GRID  BEFORE  THE  AIRCRAFT 

*  ARRIVES 

I  1= (XX (2) -TEMP) /XX  (41 ) 

*  COMPUTE  NUMBER  OF  MOVES  THE  TGT  WILL  MAKE  BEFORE 

*  ACFT  ARRIVES 
XX (42) =1  I 
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CALL  NXTMVE 
CALL  NSTEP 
CALL  NEWPOS 
CALL  LOCATN 
CALL  SORT 
ELSE 

*  ACFT  ARRIVES  BEFORE  TGT  COMPLETES  FIRST  MOVE 
NOMOVE=l 

CALL  NEWPOS 
DO  90,  1=1, MAXBOX 

TMF'LOC  ( I )  =F‘OST  I N  ( I ) 

90  CONTINUE 

CALL  SORT 

*  SET  THE  TIME  OF  THE  TGT'S  NEXT  MOVE  PREDICTED 

*  BY  THE  ACFT 

XX (22)= (XX (38) /2. 0) +TNOW 
ENDIF 

♦SET  VALUES  FOR  OUTPUT  STATEMENTS 
I=XX(11)+1 
PRETK=XX (2) -XX (3) 

T I ME=TNOW-X  X (49) 

GRIDUF-XX  (15) 

WRITE (WTUNIT, 100) I , PRETK , GR I DUP , T I ME 
100  F ORMAT ( / '  ','FOR  MISSION  NO  ,13,'  AN  UPDATE  WAS 

RECEIVED  , F  9 . 3 . 

1  MINUTES  BEFORE  TAKEOFF'/'  ', ' THE  TARGET  WAS  DETECTED 

IN  GRID  NO  ' 

1,13,'  AND  HAD  BEEN  THERE  ,F9.3,'  MINUTES'/) 

RETURN 

♦♦♦♦♦♦♦♦♦♦♦*♦*♦♦♦♦*♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦«♦♦♦♦♦♦»♦♦♦♦♦ 

♦EVENT  2  COMPUTES  THE  TARGET  LOCATION  DISTRIBUTION  PREDICTED 
♦FOR  THE  AIRCRAFT  ARRIVAL  TIME  IF  THE  UPDATE  15  AFTER  THE 
♦ACFT  TAKEOFF  TIME.  THE  DISTRIBUTION  IS  SORTED  FROM  HIGH  TO 
♦LOW  PROBABILITIES  SO  THE  ACFT  SEARCHES  THE  HIGHEST 
♦PROBABILITY  AREAS  FIRST 

*  ***♦♦♦♦*♦♦*♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦■#♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

♦COMPUTE  THE  TIME  IT  TAKES  THE  TGT  TO  MOVE  OUT  OF  THE  FIRST 

♦  GR  I D 

2  TEMF- (XX (21 ) /2.0) +XX (20) 

IF (TEMP. LE. XX (2) )  THEN 

♦  TGT  MOVES  OUT  OF  THE  FIRST  GRID  BEFORE  THE 

♦  AIRCRAFT  ARRIVES 

I  1= (XX (2) -TEMP) /XX (4  1 ) 

♦  COMPUTE  NUMBER  OF  MOVES  THE  TGT  WILL  MAKE  BEFORE 

♦  ACFT  ARRIVES 
XX (42) =1  I 
CALL  NXTMVE 
CALL  NSTEP 
CALL  NEWPOS 
CALL  LOCATN 
CALL  SORT 

ELSE 

♦  AIRCRAFT  ARRIVES  BEFORE  THE  TGT  MOVES  OUT  OF 


*  FIRST  GRID 
NOMOVE= 1 
CALL  NEWF'OS 

DO  95,  1=1 ,MAXBOX 

TliPLOC  ( I )  =F'OST  I N  ( I ) 

95  CONTINUE 

CALL  SORT 

*  COMPUTE  TIME  OF  NEXT  MOVE  PREDICTED  BY  THE 

*  AIRCRAFT 

XX (22) = (XX (38) / 2 . O ) +  TNOW 
END  IF 

♦  SET  VALUES  FOR  THE  OUTPUT  STATEMENT 

I  =  X  X (11) +1 
T1ME=TN0W— XX (49; 

GRIDUP=XX ( 15) 

WRITE  (WTUNIT  ,  110)  I  ,  GR I  DUF' ,  T  I  ME 
110  FORMAT (/'  , ' FOR  MISSION  NO  ',13,'  AT  THE  UPDATE 

TIME  THE  TARGET  W 

IAS  IN'/'  ‘,'GRID  NO  ',13,'  AND  HAD  BEEN  THERE  FOR  ', 

F9. 3 , '  MINUTES 

1  '  /) 

RETURN 

************************************************************ 
♦EVENT  3  RETURNS  TO  THE  SLAM  PROGRAM  WITH  THE  NEXT  GRID  FOR 
♦THE  AIRCRAFT  TO  SEARCH  AND  ALSO  RETURNS  THE  PROBABILITY  OF 
♦TARGET  DETECTION  GIVEN  THE  TARGET  IS  IN  THAT  GRID 
******♦**♦♦♦**♦♦♦♦****  ^*^ ♦♦♦♦** ***************************** 
♦UPDATE  THE  SORTED  VECTOR  POSITION  COUNTER 
3  SRTCTR=SRTCTR+1 

♦UPDATE  THE  NUMBER  OF  GRIDS  SEARCHED 
XX (10)=XX (10) +1.0 

IF (SRTPRB (SRTCTR, 1) .GT.O.O)  THEN 

♦  SEARCH  THE  NEXT  GRID  ON  THE  LIST 
XX (13) =DETECT (SRTPRB (SRTCTR, 2) > 

XX  (8)  =SRTF'RB  (SRTCTR,  2) 

ELSE 

*  NEXT  GRID  HAS  A  F'ROBAB I L I TY=0  OF  TARGET  PRESENCE 

♦  GO  BACK  TO  THE  FIRST  ONE  ON  THE  LIST 
SRTCTR= 1 

XX (13) =DETECT (SRTPRB (SRTCTR, 2) ) 

XX  (8)  =SRTF'RB  (SRTCTR  ,  2) 

END  I F 
T IME=TNQW 
I  =  X  X  (  1 0 ) 

♦SAVE  GRID  BEING  SEARCHED  TO  BE  PR I NTED  OUT  LATER 
GDSRCH ( I , 1 ) =XX (8) 

GDSRCH (1,2) =T IME 
RETURN 

************************************************************ 
♦EVENT  4  RECOMPUTES  THE  UNCONDITIONAL  LOCATION  DISTRIBUTION 
♦WHEN  THE  TARGET  IS  ASSUMED  TO  HAVE  MOVED  DURING  THE  SEARCH 
************************************************************ 

♦  UPDATE  THE  NEXT  PLED  I CTED  MOVE  TIME 
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4 


XX (22) =XX (22) +X X  (41 ) 

IF (NOMOVE. EQ. 1 )  THEN 

*  COMPUTE  APPROPRIATE  N-STEP  TRANSITION  MATRIX 

*  FOR  ONE  MORE  MOVE 
XX (42) =0. O 

CALL  NSTEP 
N0M0VE=0 
ELSE 

XX (42) =1 . 0 
CALL  NSTEP 
ENDIF 

CALL  LOCATN 
CALL  SORT 

*  WRITE (WTUNIT, 120) 

*120  FORMAT  (//'  '  ,  ' THE  TARGET  MOVED  DURING  THE  SEARCH'/ 

'  ' , ' THE  NEW  TAR 

*  1  GET  LOCATION  DISTRIBUTION  IS:  '//'  ',12X,'3RID 

NUMBER '  , 1  OX ,  ' PROB  0 

*  IF  TARGET  PRESENCE'/) 

*  DO  20,  1=1, MAXBOX 

*  IF  (SRTPRB (1,1). GT. 0)  THEN 

*  WRITE (WTUNIT, 130) SRTPRB (1,2) , SRTPRB (1,1) 

*130  FORMAT ( 12X , F5. 1 , 23X , F5. 3) 

*  ENDIF 

*20  CONTINUE 

♦RESET  NEXT  POSITION  IN  SORTED  VECTOR  TO  SEARCH  THE  TOP  OF 
*THE  LIST 

SRTCTR=0 

RETURN 

************************************************************ 
♦EVENT  5  RESETS  APPROPRIATE  VARIABLES  TO  GET  READY  FOR 
♦ANOTHER  MISSION 

*  ♦-a-***********-********************************************** 
♦RESET  POSITION  ON  SORTED  VECTOR  TO  TOP  OF  LIST 

5  SRTCTR=0 

♦RESET  NUMBER  OF  GRIDS  SEARCHED  FOR  THE  CURRENT  MISSION 
X  X ( 1 O ) =0 . 0 

♦RESET  FUEL  LOAD  ON  THE  AIRCRAFT 
XX (6) =XX (40) 

* RE SET  N-STEP  TRANSITION  MATRIX  TO  THE  ONE-STEP  TRANSITION 
♦MATRIX 

DO  40,  1=1, MAXBOX 

DO  30,  J=l, MAXBOX 

ALLSTP  ( I , J ) =QNESTP <  I  ,  J  ) 

30  CONTINUE 

40  CONTINUE 

♦RESET  STATUS  OF  ALL  GRIDS  TO  THE  ORIGINAL  VALUES 
Du  50,  1=0, MAXBOX 

STATUS  (  I  )  =TF’ST AT  (  I  > 

50  CONTINUE 

♦RESET  FIRST  MOVE  COUNTER  TO  INDICATE  FIRST  MOVE 
XX ( 19) =0. O 

♦SET  STATUS  OF  INITIAL  GRID  TO  3 


XX (24) =3. 0 

♦RESET  TARGET  DESTRUCTION  MARKER  TO  "DO  NOT  DESTROY" 

X  X ( 26 ) =0 . 0 

♦RESET  TARGET'S  LAST  POSITION  AS  GIRD  0 
LSTGRD=0 

♦COMPUTE  ROW  OF  CENTER  GRID  FOR  INITIAL  POSITION  OF  TARGET 
1 1= (XX (25) /2. 0) +1 . 0 
XX ( 17) =1  I 

♦SET  INITIAL  COLUMN  POSITION  EQUAL  TO  THE  INITIAL  ROW  FOR 
♦TARGET 

XX ( 18) =XX (17) 

♦COMPUTE  CENTER  GRID  NUMBER  FOR  TARGET  INITIAL  POSITION 
1 1= (MAXBOX/2) +1 . 0 
XX ( 15) =1 1 

♦SET  INITIAL  POSITION  VECTOR  TO  THE  INITIAL  GRID  LOCATION 
DO  60,  1=1, MAX BOX 

IF ( I . EQ. XX  ( 15)  )  THEN 
POST  IN ( I ) = 1 . O 
ELSE 

POST  I N  ( I )  =0.  O 
ENDIF 

60  CONTINUE 

I  I  =  X  X ( 1 1 ) +1 

♦PRINT  MESSAGE  FOR  MISSION  NUMBER  THAT  IS  STARTING 
WRITE (NPRNT, 140) 

140  FORMAT ('  ',72 ('♦')/) 

IF (XX ( 10) . EQ.O)  WRITE (NPRNT , 150) 1 1 
150  FORMAT ('  ',  AIRCRAFT  NUMBER  ',13,'  SEARCHING  NOW:'/) 

RETURN 

************************************************************ 
♦EVENT  6  RETURNS  THE  STATUS  OF  THE  GRID  WHICH  THE  TARGET  IS 
♦TRYING  TO  MOVE  INTO  NEXT 

************************************************************ 

6  I I=XX (16) 

XX (24) =STATUS (II) 

RETURN 

************************************************************ 
♦EVENT  7  CHANGES  THE  STATUS  OF  THE  NEW  TARGET  LOCATION  TO  1 
♦AND  THE  STATUS  OF  THE  OLD  GRID  BACK  TO  2 
************************************************************ 

7  I  I  =  X  X  (  1 5 ) 

STATUS (I  I) =1.0 
STATUS (LSTGRD) =2.0 
LSTGRD= 1 1 

RETURN 

************************************************************ 
♦EVENT  8  PRINTS  A  MESSAGE  WITH  INFORMATION  ABOUT  THE  NEW  RUN 
************************************************************ 

8  I  I  =  X  X  ( 2 ) 

WRITE (WTUNIT, 160) 

160  FORMAT ( '  ', 72 ('*')/ ) 

♦PRINT  MESSAGE  SHOWING  THE  NEW  UPDATE  TIME  FOR  THE  CURRENT 


170 


WRITE (WTUNIT, 170)11 

FORMAT (/'  ',  UPDATE  RECEIVED  ,14,'  MINUTES  BEFORE 

SEARCH  BEGINS 

1:  '  /) 

I  I =X X (47) 

WRITE (NPRNT, 180) 

180  FORMAT ( '  ', 72 ('*')/ ) 

♦PRINT  MESSAGE  SHOWING  THE  NEW  RUN  NUMBER 
WRITE (NPRNT , 190)  1 1 

190  FORMAT ('  '  ,  '  FOR  RUN  NUMBER  ,13,':'/) 

RETURN 

**♦***-***■»(•-•«■*•»♦♦♦♦**♦♦♦♦*♦♦♦♦♦♦♦***♦♦♦♦**♦♦**♦*♦♦♦♦♦♦♦♦**♦♦♦♦ 
♦EVENT  9  PRINTS  OUT  MESSAGES  TO  TRACK  THE  MOVEMENT  OF  THE 
♦TARGET  AS  WELL  AS  THE  TIME  OF  THE  MOVE 
♦♦♦**♦**♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦*♦ 
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200 


2  1 0 


22  V 


230 


I  I  =  XX (24) 

T I ME=TNOW 
XX  (49) =TNOW 
IF (II. EQ. O)  THEN 
I=XX (16) 

WR I TE ( NPRNT , 200 ) I , T I ME 
FORMAT ( '  , ' THE  TARGET 

1  F  1  1 . 3  ) 

ELSEIF ( I  I . EQ.  1 )  THEN 
I=XX ( 15) 

WRITER  NPRNT , 2 1 0 )  I , T I  ME 
FORMAT ('  ,  THE  TARGET 

1  I  ME  '  , F 1 1 . 3 ) 

ELSEIF ( I  I . EQ. 2)  THEN 
I  =  X  X ( 15) 

WRITE (NPRNT, 220) I , TIME 
FORMAT (  ' ,  THE  TARGET 

11.3,  TO  SET  UP  '  ) 

ELSEIF (II. EQ. 3)  THEN 
I=XX < 15) 

WRITE ( NPRNT , 230 )  I , T I  ME 
FORMAT (  ' , ' THE  TARGET 

1 ,F1 1 .3) 

END  I F 


TRIED  TO  MOVE  TO 
13,  ' 


GRID  '  , 

AT  TIME  ' , 


STARTED  MOVING  THROUGH 

GRID  ' , 13, '  AT  T 


ENTERED  GRID  NUMBER  ',13, 

AT  TIME  ' ,F1 


STARTED  IN  GRID  NUMBER  ', 
13, '  AT  TIME 


RETURN 


♦EVENT  10  PRINTS  OUT  A  MESSAGE  TO  INDICATE  WHY  THE  PREVIOUS 
♦MISSION  ENDED 

*<•»***■»<*■»•«**♦■«•«**■»**■»***  +  ■»♦******#■>*****♦*  +  «■»■♦*•«♦***♦*♦*■«■+* 


10  I  I - X X ( 46 ) 

1  =  XX  ( 1 1 ) +  1 

GOTO  ( 50 1 , 502 , 503 , 504 )  ,  I  I 
501  WRITE (WTUNIT ,240) I 

240  FORMAT ( /  'MISSION  NO  .13.'  ENDED  DUE  TO  NO  SEARCH 


TIME  ALLOWED 


1  '  /) 

RETURN 

502  WRITE (WTUNIT ,250) I 

250  FORMAT (/'  ', 'MISSION  NO  ',13,'  ENDED  DUE  TO  ALLOTTED 

SEARCH  TIME  E 

1XPIRED '  /) 

RETURN 

503  WR I TE ( WTUN I T , 260 ) I 

260.  FORMAT ( / '  ', 'MISSION  NO  ',13,'  ENDED  DUE  TO  AIRCRAFT 

BELOW  MINIMUM 

1  FUEL ' / > 

RETURN 

504  WR I TE ( WTUN I T , 270 ) I 

270  FORMAT  (/'  ', 'MISSION  NO  ',13,'  ENDED  DUE  TO  TARGET 

DETECTION '  /  ) 

RETURN 

************************************************************ 
♦EVENT  11  PRINTS  OUT  THE  SORTED  UNCONDITIONAL  LOCATION 
♦DISTRIBUTION  COMPUTED  FROM  THE  UPDATED  INFORMATION  THAT  THE 
♦AIRCRAFT  WILL  USE  TO  SEARCH  FOR  THE  TARGET  INITIALLY 
************************************************************ 

11  I JK= 1 

*11  WRITE (WTUNIT, 280) 

*280  FORMAT ( / '  ' , ' I N I T I ALLY ,  THE  TARGET  LOCATION 

DISTRIBUTION  ASSUMED  B 

*  1Y'/'  ,  '  THE  AIRCRAFT  IS  (FROM  HIGH  TO  LOW 

PROBABILITIES):'//'  ',10 

*  IX, 'GRID  NUMBER' , 10X, ' PROB  OF  TARGET  PRESENCE') 

*  DO  70,  1=1, MAXBOX 

*  IF  ( SRTPRB (1,1) . GT . O )  THEN 

*  WRITE (WTUNIT, 290) SRTPRB (I, 2) , SRTPRB (1,1) 

*290  FORMAT (14X,F5. 1 ,23X,F5.3) 

*  END IF 

*70  CONTINUE 

RETURN 

************************************************************ 
♦EVENT  12  PRINTS  OUT  THE  LIST  OF  GRIDS  SEARCHED  B"  THE 
♦AIRCRAFT  ON  HIS  MISSION  IN  THE  ORDER  THEY  WERE  SEARCHED  AND 
♦ALSO  PRINTS  OUT  THE  TIME  THE  SEARCH  IN  EACH  GRID  STARTED 
************************************************************ 

12  WRITE (WTUN IT, 300) 

300  FORMAT ('  ' , ' THE  AIRCRAFT  SEARCHED  THE  FOLLOWING  GRIDS 

IN  THE  LISTE 

ID  ORDER:  '//'  ','ACFT  NO '  , 10X ,  ' GRID  SEARCHED 1 1 X ,  ' 

TIME  SEARCH  BEG 

IAN  '  ) 

I I=GDSRCH (1,1) 

I=XX  (  1 1 ) 

WRITE (WTUNIT, 3 10) I , I  I , GDSRCH (1,2) 

310  FORMAT ('  '  ,2X, 13, 17X, 1 3 , 20 X , F 1 1 . 3 ) 

K.=  X X  (10) 

DO  80,  J  =  2 , K 


1 1 =GDSRCH ( J , 1) 

WRITE (WTUNIT , 320 ) II ,GDSRCH(J,2) 

320  FORMAT ('  ' , 22 X ,  13 , 20 X , F 1 1 . 3) 

80  CONTINUE 

WRITE (WTUNIT , 330) 

330  FORMAT ( '  ' , 72 ('*')/ ) 

RETURN 

*******************  *********  *******  *************  ************ 
♦EVENT  13  PRINTS  A  MESSAGE  INDICATING  THE  AIRCRAFT  FOUND  THE 
♦CORRECT  GRID  BUT  THE  SENSORS  DID  NOT  SEE  THE  TARGET 
************************************************************ 
13  I I  =  XX (15) 

WRITE (WTUNIT, 340) II 

340  FORMAT (/'  ' , ' THE  AIRCRAFT  SEARCHED  THE  GRID  THAT 

CONTAINED '  /  '  '  ,  'T 

1HE  TARGET  BUT  STILL  MISSED'/'  ','(IN  GRID  NUMBER  ', 

13,  '  )  '  /  ) 

RETURN 

************************************************************ 

END 

************************************************************ 
*********************** ****************** ******************* 
SUBROUTINE  NSTEP 

************************************************************ 
♦SUBROUTINE  NSTEP  COMPUTES  THE  APPROPRIATE  N-STEP  TRANSITION 
♦MATRIX  CORRESPONDING  TO  THE  NUMBER  OF  MOVES  THE  TARGET  CAN 
♦MAKE  FROM  THE  TIME  OF  THE  LAST  UPDATE  UNTIL  THE  TIME  THE 
♦SEARCH  BEGINS 

************************************************************ 

♦DECLARATIONS 

REAL  TEMP 

INTEGER  NBRMVE , MAXBOX 
♦COMMON  STATEMENTS 

COMMON /SCOM 1 / ATR I B  < 1 00 )  , DD ( 1 00 )  , DDL ( 1 00 )  , DTNOW ,11, MFA , 

MSTOP , NCLNR 

1  ,  NCRDR , NF’RNT  ,  NNRUN  ,  NNSET  ,  NTAF’E ,  SS  ( 100 )  ,SSL  (  100)  ,  TNEXT  , 

TNOW, XX (100) 

* 

COMMON/ ARAYS/ONESTP ( 1 69 , 1 69 ) , STATUS (0: 169) ,ACFTDT(169) 
1 ,SRTPRB( 169,2) , DETECT (169) , ALLSTP ( 169 , 169 ) 

1 , POST  I N ( 1 69 )  , TMPMAT (169,169), TMPLOC (169) 

1  ,  TMF'SRT  (169)  ,  TPSTAT  (0:169) 

* 

COMMON/ ARYSZE/MAXBDX 

♦NUMBER  OF  MOVES  THE  TARGET  MAKES  AFTER  INITIAL  MOVE 
1 I  =  X  X (42) 

DO  60,  NBRMVE= 1,11 

*  MULTIPLY  THE  N-STEP  TRANSITION  MATRIX  BY  THE 

*  ONE-STEP 

DO  30,  1=1, MAXBOX 

DO  20,  J=l, MAXBOX 
TEMP=0. 0 

DO  10,  K=l, MAXBOX 


TEMP=TEMF'+  (ONESTP  ( I  ,  K)  *ALLSTF‘  (K ,  J )  ) 

10  CONTINUE 

TMPMATd  ,  J)=TEMP 
20  CONTINUE 

30  CONTINUE 

DO  50,  1=1, MAXBOX 

DO  40,  J=1 , MAXBOX 

ALLSTP ( I , J ) =  TMF'MAT ( I , J ) 

40  CONTINUE 

50  CONTINUE 

60  CONTINUE 

RETURN 
END 

************************************************************ 
************************************************************ 
SUBROUTINE  NEWPOS 

******************  ****************************************** 
♦SUBROUTINE  NEWPOS  RESETS  THE  INITIAL  TARSET  POSITION  VECTOR 
*T0  THE  NEW  POSITION  ESTABLISHED  FROM  THE  UPDATED 
♦INTELLIGENCE  INFORMATION 

************************************************************ 

♦DECLARATIONS 

INTEGER  MAXBOX 
♦COMMON  STATEMENTS 

COMMON /SCOM 1 / ATR I B ( 1 00 )  ,DD(100> , DDL (100)  , DTNOW , 1 1 , MFA , 

MS TOP , NCLNR 

1  ,  NCRDR  ,  NF'RNT  ,  NNRUN  ,  NNSET  ,  NT  APE ,  SS  ( 1 00 )  ,  SSL  ( 1 00 )  ,  TNEXT  , 

TNOW,XX (100) 

* 

COMMON/ ARAYS/ONESTF'  <  1 69 ,169)  , STATUS  <0:  169)  ,ACFTDT<169> 
1 ,SRTPRB(169,2) , DETECT (169) , ALLSTP ( 169 , 169 ) 

1 ,F0STIN(169) ,TMPMAT< 169, 169) ,TMPLDC(169) 

1  ,  TMF'SRT  (169)  ,  TPSTAT  (0:  169) 

* 

COMMON/ ARYSZE/MAXBOX 

♦SET  THE  VECTOR  POSITION  TO  1  FDR  THE  CURRENT  T6T  LOCATION 
♦AND  O  FOR  OTHERS 

DO  10,  1=1, MAXBOX 

IF ( I . EQ. XX ( 15) )  THEN 
POST IN ( I ) =1 . 0 
ELSE 

POST  IN ( I ) =0. 0 
END  IF 

10  CONTINUE 

RETURN 
END 

************************************************************ 
************************************************************ 
SUBROUTINE  LOCATN 

************************************************************ 
♦SUBROUTINE  LOCATN  COMPUTES  THE  NEW  UNCONDITIONAL  LOCATION 
♦DISTRIBUTION  OF  THE  TARGET 

************************************************************ 


1  15 


♦DECLARATIONS 

REAL  TEMP 
INTEGER  MAXBOX 
♦COMMON  STATEMENTS 

COMMON/ ARAYS/ONESTP ( 169 , 169 ) , STATUS (0: 169) ,ACFTDT(169) 
1 , SRTPRB (169, 2) , DETECT (169) , ALLSTP (169,169) 

1 , POST  IN (169)  , TMPMAT  < 1 69 , 1 69 )  ,TMPL0C(169) 

1 , TMPSRT (169)  , TPST  AT (0:169) 

* 

COMMON/ ARYSZE/MAXBOX 

♦MULTIPLY  THE  UPDATED  POSITION  VECTOR  BY  THE  APPROPRIATE 
♦N-STEP  TRANSITION  MATRIX 
DO  20,  1=1, MAXBOX 
TEMP=0. 0 

DO  10  J=l, MAXBOX 

TEMP=TEMP+ (POST  IN (J ) ♦ALLSTP ( J  ,  I  >  ) 

10  CONTINUE 

TMPLOC ( I ) =TEMP 
20  CONTINUE 

RETURN 
END 

♦  ♦♦♦♦♦♦♦♦♦*♦♦♦♦**♦**♦*****■#♦♦**♦  **************************** 
♦  ♦*♦♦*♦♦*♦♦♦*♦*♦****•#■#**■*•*•#*♦*♦****#•♦#■#*■*■#■*•#******■*■*•*#•*■*••***# 
SUBROUTINE  SORT 

*♦*♦♦*♦**■*•*■*♦***■*****•*** ****^***************-**************** 
♦SUBROUTINE  SORT  SORTS  THE  UNCONDITIONAL  LOCATION 
♦DISTRIBUTION  FROM  HIGH  TO  LOW  PROBABILITIES  OF  TARGET 
♦PRESENCE  IN  EACH  GRID 

*****•*♦♦*♦  *************************************************** 
♦DECLARATIONS 

INTEGER  MAXBOX, RDUN I T , WTUNI T 
♦COMMON  STATEMENTS 

COMMON/ ARAYS/ONESTP ( 169, 169) , STATUS (O: 169) ,ACFTDT(169> 
1 ,SRTPRB< 169,2) , DETECT (169) , ALLSTP ( 169, 169) 

1 ,P0STIN(169) , TMPMAT ( 169, 169) ,TMPL0C(169) 

1 , TMPSRT (169) , TPSTAT ( 0 : 1 69 ) 

* 

COMMON /RDWRT /RDUN I T , WTUNI T 

♦ 

COMMON /ARYSZE/MAXBOX 

♦ 

DO  10,  1=1, MAXBOX 

TMPSRT ( I ) =TMPLOC ( I ) 

10  CONTINUE 

DO  30,  1=1, MAXBOX 

♦  INITIALLY  PUT  THE  FIRST  GRID  INFORMATION  AT  THE 

♦  TOP  OF  THE  LIST 
SRTPRB (1,1) =TMPSRT ( 1 ) 

SRTPRB ( I ,2) =1.0 

DO  20,  J=2, MAXBOX 

♦  IF  THE  PROBABILITY  OF  THE  TARGET  BEING  IN  THE 
NEXT  GRID  ON  THE  LIST  IS  GREATER  THAN  THE 
PROBABILITY  ALREADY  ON  THE  LIST,  SWITCH  THEM 


* 

# 


I F  ( SRTPRB  (1,1).  LT .  TMF'SRT  ( J  )  )  THEN 
SRTPRB  (1,1)  =TMF'SRT  ( J  ) 

SRTPRB (I, 2) =J 
END  IF 

20  CONTINUE 

*  MARK  THE  SPOTS  ALREADY  PUT  ON  THE  SORTED  LIST 
TMF'SRT  (SRTPRB  ( I  ,2)  )=-1.0 

30  CONTINUE 

RETURN 
END 

***  ******************  ******  ****************  ***************** 
************************************************************ 
SUBROUTINE  NXTMVE 

************************************************************ 
♦SUBROUTINE  NXTMVE  COMPUTES  THE  TIME  THE  AIRCRAFT  ASSUMES 
*THE  TARGET  WILL  MOVE  AGAIN 

************************************************************ 
COMMON /SCOM 1  /  ATR I B ( 100) ,DD(100> ,DDL(100) , DTNOW ,11, MFA , 

MSTOF' ,  NCLNR 

1  ,  NCRDR  ,  NF'RNT  ,  NNRUN  ,  NNSET ,  NT  APE  ,  SS  ( 1 00 )  ,  SSL  ( 1 00 )  ,  TNEXT  , 

TNOW , X  X ( 1 00 ) 

* 

FSTMVE=XX (20) + (XX (21 ) /2. 0) 

I I=XX (42) 

XX (22) =FSTMVE+ (  ( X X ( 2 1 ) +X X ( 38 )  ) * ( 1 1  +  1 )  ) +  TNOW 

RETURN 

END 

************************************************************ 
************************************************************ 
SUBROUTINE  OTF'UT 

***************************  ********************************* 

♦  SUBROUTINE  OTF'UT  PRINTS  OUT  ANY  INFORMATION  DESIRED  OTHER 
♦THAN  THE  STANDARD  SLAM  OUTPUT 

************************************************************ 
COMMON,  SCOM 1/ ATR IB ( 100)  , DD ( 100)  , DDL ( 100)  , DTNOW, 1 1 , MFA , 

MSTOF, NCLNR 

1  ,  NCRDR,  NF'RNT,  NNRUN,  NNSET  , NTAFE  ,  SS  (  100)  ,SSL  (100)  ,  TNEXT, 

TNOW , X  X ( 1 00 ) 

* 

RETURN 

END 

************************************************************ 

************************************************************ 


SLAM  Program 

GEN, QLVNICK, RELOCATABLE  TARGETS,  1 2/ 1 3/85 , 3 , N , N; 

LIMITS, 0,0, 5; 

; *************************** ********************** ********** 
■  *********************************************************** 
;*  THIS  PROGRAM  COMPUTES  THE  PROBABILITY  OF  DETECTING  * 

;*  A  STRATEGIC  RELOCATABLE  TARGET  (SRT )  BY  A  STRATEGIC  * 
;*  BOMBER.  INPUT  VARIABLES  AND  THEIR  DESCRIPTIONS  ARE  * 
;*  IN  THE  FORTRAN  "PROGRAM  MAIN".  AFTER  SOME  VARIABLES  * 
;*  ARE  INPUT  AND  OTHER  VARIABLE  VALUES  ARE  INITIALIZED  * 
;*  WITH  THE  INPUT  VALUES,  THE  COMBINATION  NETWORK-  * 

;*  DISCRETE  SIMULATION  STARTS.  A  TARGET  MOVES  ACCORDING  * 
;*  TO  THE  DEFINED  STOCHASTIC  PROCESS  WHILE  SIMULTANEOUSLY  * 
;*  AN  AIRCRAFT  MISSION  IS  SIMULATED  FROM  THE  TIME  AN  SRT  * 
;*  LEAVES  ITS  MAIN  OPERATING  BASE  (MOB)  UNTIL  THE  AIRCRAFT* 
;*  EITHER  FINDS  THE  TARGET  OR  REACHES  HIS  TIME  OR  FUEL  * 
;*  RESTRICTION.  DURING  EACH  MISSION,  THE  AIRCRAFT  WILL  * 
;*  RECEIVE  AN  INTELLIGENCE  UPDATE  ON  THE  TARGET'S  LOCATION* 
;*  A  SERIES  OF  MISSIONS  ARE  FLOWN  FOR  EACH  OF  A  SET  OF  * 
;*  UPDATE  TIMES  WITH  A  SUBSEQUENT  COMPARISON  OF  THE  PROB-  * 
;*  ABILITY  OF  DETECTING  THE  TARGET  GIVEN  A  PARTICULAR  * 

;*  UPDATE  TIME.  FROM  THIS  COMPARISON,  BASED  ON  OTHER  * 

;*  INPUT  VALUES,  A  CHOICE  CAN  BE  MADE  FOR  AN  INTELLIGENCE  * 

;*  UPDATE  TIME  THAT  CAN  ENSURE  A  HIGH  PROBABILITY  OF  * 

5*  DETECTION  BY  THE  ATTACKING  AIRCREW.  * 

************************************************************ 

; *********************************************************** 


; *****  THE  NETWORK  MODEL 


****************************** 


NETWORK; 


XX  (  1) 

XX  (2) 

X  X  <  3 ) 
XX  (4) 

X  X  <  5 ) 
XX  (6) 
XX  (7) 
XX  (8) 
XX  (9) 

X  X  ( 1 0 ) 

X  X  <  1  1  ) 


DEFINITION  OF  VARIABLES 


********************* 


TIME  FROM  WHEN  TARGET  LEAVES  MOB  (TIME  0)  UNTIL 
AIRCRAFT  T/0 

TIME  REMAINING  UNTIL  BEGIN  SEARCH  AFTER  UPDATE 
RECEIVED 

TIME  FROM  TAKEOFF (T/0)  UNTIL  BEGIN  SEARCH 

TIME  FROM  INTELL  INFO  REC ( T I  ME  O)  UNTIL  RECEIVE 

UPDATE 

TIME  IT  TAKES  AIRCRAFT  TO  SEARCH  ONE  GRID 
FUEL  LOAD  ON  THE  AIRCRAFT  AT  T/0  TIME 
MINIMUM  FUEL  WHEN  AIRCRAFT  MUST  STOP  SEARCH 
GRID  NUMBER  AIRCRAFT  IS  SEARCHING 
NUMBER  OF  GRIDS  THE  AIRCRAFT  CAN  SEARCH  IN 
ALLOTTED  TIME 

THE  NUMBER  OF  GRIDS  SEARCHED  SO  FAR  BY  THE 
CURRENT  AIRCRAFT 

NUMBER  OF  AIRCRAFT  FLOWN  WITH  CURRENT  UPDATE 


& 


TIME 

XX (12)  NUMBER  OF  AIRCRAFT  THAT  MISS  THE  TARGET  FOR  ONE 

UPDATE  TIME 

XX (13)  PROBABILITY  OF  AIRCRAFT  DETECTING  THE  TARGET  IN 

THE  GRID  BEING  SEARCHED  GIVEN  THE  TARGET  IS 
THERE 

XX (14)  MAXIMUM  NUMBER  OF  AIRCRAFT  TO  FLY  PER  UPDATE 

TIME 

XX (15)  GRID  NUMBER  OF  TARGET  LOCATION 

XX (16)  GRID  NUMBER  TO  MOVE  THE  TARGET  TO  NEXT  WHILE 

THE  STATUS  OF  THE  GRID  IS  BEING  CHECKED  FOR  AN 
ACCEPTABLE  MOVE 

XX  (17)  ROW  NUMBER  OF  THE  TARGET  LOCATION 

XX (18)  COLUMN  NUMBER  OF  THE  TARGET  LOCATION 

XX  (19)  MARKER:  O-FIRST  MOVE  OF  TARGET  1 -SECOND  OR 

LATER  MOVE 

XX  (20)  TIME  IT  TAKES  TARGET  TO  MOVE  HALF  WAY  ACROSS 

ONE  GRID 

XX  (21)  TOTAL  TIME  TARGET  STAYS  IN  ONE  LOCATION 

XX (22)  TIME  PREDICTED  BY  AIRCRAFT  FOR  THE  TARGET'S 

NEXT  MOVE 

XX (23)  RANDOM  NUMBER  FOR  DIRECTION  OF  TARGET'S  NEXT 

MOVE 

XX (24)  STATUS  OF  GRID  CHOSEN  FOR  NEXT  LOCATION 

0  -  STAY  OUT 

1  -  PASS  THROUGH  ONLY 

2  -  RELOCATE 

XX  (25)  NUMBER  OF  GRIDS  IN  ONE  ROW/COLUMN  OF  REGION 

XX  (26)  MARKER:  O-CONTINUE  MISSION  1-END  MISSION 

XX  (27)  RANDOM  NUMBER  USED  IN  TARGET  IDENTIFICATION 

XX (28)  NUMBER  OF  AIRCRAFT  THAT  DETECT  THE  TARGET  PER 

UPDATE  TIME 

XX  (29)  PROBABILITY  OF  DETECTION  PER  UPDATE  TIME 

XX (30)  DIFFERENCE  BETWEEN  THE  SYSTEM  CYCLE  TIME  AND 

LATEST  UPDATE 

XX  (31)  SYSTEM  CYCLE  TIME  (TIME  FROM  INTELLIGENCE 

INFORMATION  COLLECTED  UNTIL  WEAPON  OVER  TARGET) 
XX  (32)  TIME  INCREMENT  USED  BETWEEN  UPDATE  TIMES  BEING 

CHECKED 

XX  (33)  MARKER:  O-NEXT  RUN  NOT  LAST  RUN 

1 -NEXT  RUN  IS  LAST  RUN 

XX  (35)  TIME  BETWEEN  ACTUAL  END  OF  MISSION  AND  LATEST 

END  POSSIBLE 

XX  (36)  TIME  BETWEEN  DECREASES  IN  FUEL  ON  BOARD  THE 

AIRCRAFT 

XX  (37)  AMOUNT  OF  FUEL  TO  DECREASE  THE  AIRCRAFT  FUEL 

EACH  INCREMENT 

XX (38)  TIME  FOR  SRT  TO  MOVE  ONE  GRID 

XX (39)  GRID  SIZE  IN  MILES 

XX (40)  SAVES  INITIAL  FUEL  LOAD  BETWEEN  RUNS 

XX  (41)  TIME  UNIT  OF  STOCHASTIC  PROCESS 

XX  (42)  NUMBER  OF  MOVES  SRT  MAKES  AFTER  FIRST  MOVE 

BEFORE  AIRCRAFT  ARRIVES  TO  BEGIN  SEARCH 


rj)njn’j^r*\r*r,y'y  *y  '.v  uTir'm'Tw.TT  "  yj.”  i1'  v  i  n  ^  ■  i 


;  XX (46)  MARKER  FOR  REASON  FOR  MISSION  TERMINATION 

XX  (47)  KEEPS  TRACK  OF  THE  RUN  NUMBER 

XX (48)  TIME  BEFORE  TAKEOFF  THAT  AN  UPDATE  MAY  BE 

RECEIVED 

XX (49)  SIMULATION  TIME  TARGET  MOVES  INTO  A  NEW  GRID 

*****  CREATE  AIRCRAFT  AND  FLY  IT  TO  THE  SEARCH  AREA  *** 

CREATE;  CREATE  AIRCRAFT 

ACT; 

DQAGN  ASSIGN, XX ( 1 1 ) =0, 

XX (12) =0, 

XX (28) =0, 

ACT; 

EVENT, 8, 1; 

ACT; 

ACFT  EVENT, 5, 2; 


l; 


RESET  ACFT  COUNTERS 
FOR  NEW  RUN 

RESET  OTHER  VAR¬ 
IABLES  FOR  NEW  RUN 


UPDATE  VARIABLES 
FOR  A  NEW 
AIRCRAFT  MISSION 
CREATE  TARGET 
CONTINUE  WITH 
AIRCRAFT 

COMPUTE  TIME  UNTIL 
INTEL  UPDATE 
UPDATE  COMES  AFTER 
TAKEOFF 

UPDATE  COMES  BEFORE 
T/0,  DELAY 

;  UNTIL  UPDATE  TIME 

;****  COMPUTATIONS  FOR  UPDATE  PRIOR  TO  T/0  *************** 


ACT , , , MOVE; 

ACT; 

ASSIGN, XX (4) =XX (31 ) -XX (2) ,  1; 

ACT  ,  ,  X  X  ( 4 )  .  GT .  X  X  ( 1 )  ,  UF'DAT ; 
ACT, XX  (4)  ; 


EVENT, 1 ,1; 


ACT; 

EVENT, 1 1,1; 

ACT; 

ASSIGN, XX (48) =XX (2) -XX (3) , 
ACT, XX (48) , , FULDN ; 

ACT  ,  XX (48)  ; 


COMPUTE  NEW  TARGET 

LOCATION 

DISTRIBUTION 

PRINT  TGT  LOCATION 
DISTRIBUTION 

COMPUTE  T/0  TIME 
START  DECREASING  FUEL 
DELAY  UNTIL  T/0 
TIME 


FLY  TO  SEARCH  AREA 

****************** 


GOON, 1 ; 

ACT, XX (3) , , SERCH; 

*****  COMPUTATIONS  FOR  UPDATE  AFTER  T/0 
UP DAT  GOON, 2; 

ACT, XX < 1) ,, FULDN;  START  FUEL 

DECREASING  AFTER 
T/0 

ACT, XX (4);  FLY  TO  UPDATE  POINT 

EVENT ,2,1;  COMPUTE  TGT  LOCA¬ 

TION  DISTRIBUTION 
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ACT; 

EVENT,  11,1; 

ACT, XX (2) , , SERCH; 

;***♦  BEGIN  SEARCH  FOR  TARGET 

SERCH  EVENT, 3,1; 

ACT, , XX (9) . EQ. 0, El MSG; 

ACT, XX (5)  ,  XX  (15)  .EQ. XX  (8)  , 


PRINT  TGT  LOCATION 
DISTRIBUTION 
FLY  TO  SEARCH  AREA 
********************** 
PUT  ACFT  IN  FIRST 
GRID  to  SEARCH 
STOP,  NO  SEARCH 
TINE  ALLOWED 
CKFND;  ACFT  AND  TGT  IN 
SAME  GRID,  CHECK 
FOR  DETECTION 


ACT, XX  (5) ; 

GOON, 1 ; 

ACT, , XX ( 10) .GE. XX (9) ,E2MSG; 


STOP,  SEARCHED  ALL 
GRIDS  ALLOWED 
STOP,  ACFT  AT 
MINIMUM  FUEL 
TGT  MOVED  DURING 
SEARCH 

RECOMPUTE  LOCATION 
DISTRIBUTION 
NO  DETECTION, 

SEARCH  NEXT  GRID 

;*****  ACFT  IN  CORRECT  GRID,  CHECK  IF  DETECT  TGT  ********* 


ACT , , X X  < 6 )  . LE . XX ( 7 )  ,E3MSG; 
ACT , , TNOW . GE . X  X ( 22 )  , MVAGN ; 

ACT, , , SERCH; 


CKFND  ASSIGN, XX (27 ) =UNFRM (0. , 1 . ,4) , 
ACT , , XX (27) -LE. XX ( 13) ,E4MSG; 

? 

ACT; 

EVENT, 13,1; 

ACT, , XX ( 10) .GE. XX (9) ,E2M3G; 
ACT, , XX (6) . LE. XX (7) ,E3MSG; 
ACT , , TNOW. GE. XX (22) , MVAGN; 

5 

ACT, , , SERCH; 


1 ;  CHOOSE  RANDOM 
NUMBER 

RANDOM  NUMBER  LESS 
THAN  DETECTION 
PROBABILITY,  ACFT 
FOUND  TGT 

SENSORS  DID  NOT  SEE 
TGT,  CONTINUE 
PRINT  MESSAGE  OF 
DETECT  FAILURE 
STOP,  SEARCHED  ALL 
GRIDS  ALLOWED 
STOP,  ACFT  AT 
MINIMUM  FUEL 
TGT  MOVED  DURING 
SEARCH 

RECOMPUTE  LOCATION 
DISTRIBUTION 
NO  DETECTION, 

SEARCH  NEXT  GRID 

;****  RECOMPUTE  TARGET  LOCATION  DISTRIBUTION  ************* 
MVAGN  EVENT, 4, 1 ; 

ACT,,, SERCH;  CONTINUE  SEARCH 

WITH  NEW 

;  DISTRIBUTION 

; ****  START  TARGET  MOVEMENT  ****************************** 
MOVE  EVENT, 7,1;  RESET  STATUS  OF 

CURRENT  TGT 
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ACT, , XX (19) . EQ.O. ,FSMVE; 
ACT; 


LOCATION  TO  ONE  FOR 

NEXT  MOVE 

TGT'S  FIRST  MOVE, 

DO  NOT  DWELL 
SECOND  OR  LATER 
MOVE,  DWELL 
BEFORE  MOVE  AGAIN 


GOON, 2; 

ACT , , , DLAGN; 


ACT , STOPA  ( 1 ) ; 
FSMVE  ASSIGN, XX ( 19) =1 , 


1; 


STAY  IN  CURRENT 
GRID  FOR 

ASSIGNED  DWELL  TIME 
RESET  FIRST  MOVE 
MARKER 


ACT, , XX (26) . EQ. 1 , WTTGT ; 


STOF,  TARGET 


DETECTED 


** 

i 

** 

; 

** 

;  ** 


ACT; 

CHECK  TARGET'S  CURRENT  POSITION  BEFORE  CHOOSING  NEW  *** 
DIRECTION;  DIRECTIONS  ARE  ORIENTED  AS  UP=NORTH,  *** 
DOWN=SOUTH ,  LEFT =WEST ,  AND  RIGHT  =EAST ;  FOR  THE  5EARCH  * 


AREA  BOX,  THE 
EAST  SIDE,  THE 
ROW=SOUTH  SIDE 
TOPRW  GOON , 1 ; 

ACT , , XX (24) 
XX (24) 


LEFT 

TOP 


EDGE=WEST 

RQW=NORTH 


SIDE, 

SIDE, 


THE  RIGHT  EDGE=  *** 
AND  THE  BOTTOM  *** 


EQ . 0 . OR . 

EQ. 3 , TELME; 


ACT, , , NOGO; 
TELME  EVENT ,9,1; 


STATUS  OF 
GRID  IS  0 
STATUS  OF 
GRID  IS  1 
PRINT  MESSAGE  FOR 
TGT  MOVEMENT 


CURRENT 
OR  3 
CURRENT 
OR  2 


NOGO 


ACT; 
GOON, 1 ; 


ACT, , XX ( 17) 


GT. 

GT, 


1 , BTMRW; 
1 , URTCR; 


ACT , , XX  < 18) 

ACT , , , ULTCR; 

ULTCR  ASSIGN, XX (23; =UNFRM(0. ,3. ,4) , 


ACT, , XX (23) . LT, 


ACT , ,XX (23) . LT , 
ACT, , , SMOVE ; 
URTCR  GOON, 1 ; 

ACT , , XX ( 18) . LT 


1 , EMOVE; 
2 , 5EMVE ; 


TGT  NOT  IN  TOP  ROW 
IN  TOP  ROW,  BUT  NOT 
ON  LEFT  EDGE 
TGT  IN  UPPER  LEFT 
CORNER 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  EAST 
MOVE  SOUTH  EAST 
MOVE  SOUTH 


XX (25)  , MI  DTP; 


ACT; 

ASSIGN, XX (23) =UNFRM (O. 


3.  ,4) 


MI  DTP 


ACT, , XX (23) .LT. 1 ,WMOVE; 

ACT , , XX (23) . LT. 2 , SWMVE ; 
ACT, , , SMOVE; 

ASSIGN, XX (23) =UNFRM (0. ,5. ,4) 


IN  TOP  ROW,  NOT 
ON  RIGHT  EDGE 
IN  UPPER  RIGHT 
CORNER 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  WEST 
MOVE  SOUTH  WEST 
MOVE  SOUTH 
IN  MIDDLE  OF  TOP 
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ACT, ,XX (23) .LT. 1 ,EMOVE; 

ACT, , XX  (23) . LT. 2,SEMVE; 

ACT, , XX (23) . LT . 3 , SMOVE ; 

ACT,  , XX  (23)  . LT . 4 , SWMVE ; 

ACT , , , WMOVE ; 

BTMRW  GOGN,l; 

ACT ,  , XX  (  17)  . LT. XX (25)  , LFEGE ; 

ACT , , XX ( 18) .GT. 1 , LRTCR; 

ACT , , , LLTCR ; 

LLTCR  ASSIGN, XX (23) =UNFRM(0. ,3. ,4) ,  1; 

ACT ,  , XX  (23)  . LT.  1 , NMOVE ; 

ACT, , XX (23) . LT . 2 , NEMVE ; 

ACT , , , EMOVE ; 

LRTCR  GOON , 1 ; 

ACT , , XX ( 18) . LT. XX (25) , MDBTM ; 
ACT; 

ASSIGN, XX (23) =UNFRM(0. ,3. ,4) ,  1 ; 

ACT , , XX (23) .LT. 1 , NMOVE; 

ACT, , XX  (23) . LT . 2 , NWMVE; 

ACT, , , WMOVE ; 

MDBTM  ASSIGN, XX  (23) =UNFRM (0.  ,5.  ,4)  ,  1; 


ACT , , X  X  ( 23 )  . LT .  1 , WMOVE ; 

ACT, , XX (23) . LT. 2, NWMVE; 

ACT , , X  X ( 23 ) . LT . 3 , NMOVE ; 

ACT, , XX (23) .LT. 4, NEMVE; 

ACT, , , EMQVE; 

LFEGE  GOON , 1 ; 

ACT ,  , XX  ( 18)  . GT.  1 , RTEGE ; 

ACT , , , MIDLF ; 

MIDLF  ASSIGN, XX (23) =UNFRM (0. ,5. ,4) ,  1; 

ACT , , XX  (23)  . LT.  1 , NMOVE; 

ACT, , XX (23) .LT. 2, NEMVE; 

ACT , , XX (23) . LT . 3 , EMOVE ; 

ACT , , XX (23) . LT. 4,SEMVE; 

ACT , , , SMOVE; 

RTEGE  GOON, 1 ; 

ACT , , XX ( 18) .LT. XX (25) .M1DBX: 


ROW 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  EAST 
MOVE  SOUTH  EAST 
MOVE  SOUTH 
MOVE  SOUTH  WEST 
MOVE  WEST 

T6T  NOT  IN  BOTTOM 
ROW 

IN  BOTTOM  ROW,  NOT 
ON  LEFT  EDGE 
IN  LOWER  LEFT 
CORNER 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  NORTH 
MOVE  NORTH  EAST 
MOVE  EAST 

IN  BOTTOM  ROW,  NOT 
ON  RIGHT  EDGE 
IN  LOWER  RIGHT 
CORNER 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  NORTH 
MOVE  NORTH  WEST 
MOVE  WEST 

IN  MIDDLE  OF  BOTTOM 
ROW 

CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  WEST 
MOVE  NORTH  WEST 
MOVE  NORTH 
MOVE  NORTH  EAST 
MOVE  EAST 

TGT  NOT  ON  LEFT 
EDGE 

TGT  IN  MIDDLE  OF 
LEFT  EDGE 
CHOOSE  RANDOM 
NUMBER  FOR  MOVE 
MOVE  NORTH 
MOVE  NORTH  EAST 
MOVE  EAST 
MOVE  SOUTH  EAST 
MOVE  SOUTH 


NOT  ON  RIGHT  EDGE 


MIDRT 


MIOBX 


;  **** 

NWMVE 


ACT, , , MIDRT; 

ASSIGN, XX (23) =UNFRM(0.  ,5.  ,4)  , 

ACT , , XX (23) . LT. 1 , NMOVE ; 

ACT, , XX (23) .LT. 2, NWMVE; 

ACT , , XX (23) . LT. 3, WMOVE; 

ACT , , XX (23) . LT. 4 , SWMVE ; 

ACT , , , SMOVE ; 

ASSIGN, XX (23) — UNFRM (O.  ,6.  ,4»  , 


ACT, , XX (23) .LT. 1 , NWMVE; 

ACT ,  , X  X ( 23 )  . LT . 2 , NMOVE ; 

ACT, ,XX (23) . LT . 3 , NEMVE ; 

ACT , , XX (23) . LT. 4 , EMOVE ; 

ACT, , XX (23) . LT . 5 , SEMVE ; 

ACT , ,XX(23) .LT. 6, SMOVE; 

ACT, ,XX (23) .LT. 7, SWMVE; 

ACT, , , WMOVE; 

MOVE  TARGET  ******-******-*** 


TGT  IN  MIDDLE  OF 
RIGHT  EDGE 
1 ;  CHOOSE  RANDOM 
NUMBER 
MOVE  NORTH 
MOVE  NORTH  WEST 
MOVE  WEST 
MOVE  SOUTH  WEST 
MOVE  SOUTH 
1 .  TGT  IN  MIDDLE  OF 
SEARCH  AREA, 

CHOOSE  RANDOM 
NUMBER 

MOVE  NORTH  WEST 
MOVE  NORTH 
MOVE  NORTH  EAST 
MOVE  EAST 
MOVE  SOUTH  EAST 
MOVE  SOUTH 
MOVE  SOUTH  WEST 
MOVE  WEST 

************************ 


ASSIGN, XX ( 16) =XX ( 15) -XX (25) -1  , 
ACT; 

EVENT, 6, 1; 

ACT , , X  X ( 24 ) . EQ . 0 , TOPRW ; 

ACT; 

GOON, 2; 

ACT, , , HFMVE ; 

ACT  ,  STOF'A  (5)  ; 

GOON, 1 ; 

ACT , ,XX  (26) .EO.  1 , WTTGT ; 

ACT; 

ASSIGN, XX ( 15) =X X ( 16)  , 

XX(17)=XX(17)-1  , 

XX ( 13) =XX ( 18) -1 ,  1; 

ACT; 

EVENT, 9, 1 ; 

ACT; 

GOON, 2; 

ACT , , , HFMVE; 

ACT  ,  STOF'A  (5)  ; 

GOON, 1 ; 

ACT , , XX (26) .EQ. 1 , WTTGT; 


COMFUTE  NEW  GRID 
NUMBER 

GET  STATUS  OF  NEW 
GRID 

STATUS=0 ,  DO  NOT 
MOVE 

MOVE  TO  NEW  GRID 

MOVE  TARGET  OUT  OF 
OLD  GRID 
TO  THE  EDGE  OF 
THE  NEW  GRID 

STOP,  TARGET 
DETECTED 


CHANGE  TGT  LOCATION 
TO  NEW  GRID 

PRINT  MESSAGE  ABOUT 
TGT  MOVE 


MOVE  TARGET  TO  THE 
CENTER  OF 
THE  NEW  GRID 


STOP,  TARGET 


DETECTED 


ACT; 

GOON ,  1  ; 

ACT, ,XX  (24)  .EQ.  1 ,TOPRW;  STATUS=1 ,  MOVE 

AGAIN  BEFORE  DWELL 

ACT,,, MOVE;  STATUS=2,  DWELL 

BEFORE  MOVE  AGAIN 

;  *******■**************************•*■***************.#..**.******* 

NMOVE  ASSIGN, XX ( 16) =XX ( 15) -XX (25) ,  1;  COMPUTE  NEW  GRID 

NUMBER 

ACT; 

EVENT, 6,1;  GET  STATUS  OF  NEW 

GRID 

ACT, ,XX (24) .EQ.O,TOPRW;  STATUS=0,  DO  NOT 

MOVE 

ACT; 

GOON, 2; 

ACT , , , HFMVE;  MOVE  TARGET  OUT  OF 

OLD  GRID 

ACT,STQPA(5) ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON , 1 ; 

ACT , , XX (26) . EO. 1 , WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX ( 15) =XX ( 16) , 

XX  (17) =XX ( 17) -1 ,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT, 9,1;  PRINT  MESSAGE  ABOUT 

TGT  MOVE 

ACT; 

GOON, 2; 

ACT,,, HFMVE;  MOVE  TARGET  TO  THE 

CENTER 

ACT , STOP A ( 5 > ;  OF  THE  NEW  GRID 

GOON, 1 ; 

ACT, , XX (26) . EQ. 1 , WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

GOON, 1 ; 

ACT, , XX  (24)  . EQ.  1 , TOPRW;  STATUS=1,  MOVE 

AGAIN  BEFORE  DWELL 

ACT,,, MOVE;  STATUS=0,  DWELL 

BEFORE  MOVE  AGAIN 

j*****************************************-***********-**.**.**.* 

NEMVE  ASSIGN, XX ( 16) =XX ( 15) -XX (25) +1 ,  1;  COMPUTE  NEW  GRID 

NUMBER 

ACT; 

EVENT ,6,1;  GET  STATUS  OF  GRID 

ACT ,  , XX  (24)  . EQ. 0 , TOPRW;  STATUS=0,  DO  NOT 

MOVE 

ACT; 


GOON, 2; 

ACT , , , HFMVE;  MOVE  TARGET  OUT  OF 

OLD  GRID 

ACT ,ST0PA(5) ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON ,  1  ; 

ACT, ,XX(26> .EQ. 1,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX ( 15) =XX (16) , 

XX (17) =XX (17) -1  , 

XX(18)=XX(18)+1,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT ,9,1;  PRINT  MESSAGE  ABOUT 

TGT  MOVE 

ACT; 

GOON ,2; 

ACT,,, HFMVE;  MOVE  TARGET  TO  THE 

CENTER 

ACT , STOPA (5) ;  OF  THE  NEW  GRID 

GOON, 1; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

GOON, 1; 

ACT, ,XX (24) .EQ. 1 ,TOPRW;  STATUS=1,  MOVE 

AGAIN  BEFORE  DWELL 

ACT,,, MOVE;  STATUS=2,  DWELL 

BEFORE  MOVE  AGAIN 

; ******************************* **************************** 
EMOVE  ASSIGN ,XX(16)=XX(15)+1,  1;  COMPUTE  NEW  GRID 

NUMBER 

ACT; 

EVENT, 6,1;  GET  STATUS  OF  NEW 

GRID 

ACT, ,XX (24) .EQ.O,TOPRW;  STATUS=0 ,  DO  NOT 

MOVE 

ACT; 

GOON, 2; 

ACT,,, HFMVE;  MOVE  TARGET  OUT  OF 

OLD  GRID 

ACT , STOPA (5 ) ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON, 1 ; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX ( 15) =XX ( 16) , 

XX ( 18) =XX < 18) +1 ,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT, 9,1;  PRINT  MESSAGE  ABOUT 
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ACT; 

GOON, 2; 

ACT, , , HFMVE; 


ACT 


,  ,XX 


ACT; 

GOON,  1 ; 

ACT, , XX (24) . EQ. 1  ,TOPRW; 
ACT, , , MOVE; 


TGT  MOVE 


MOVE  TARGET  TO  THE 
CENTER 

OF  THE  NEW  GRID 

STOP,  TARGET 
DETECTED 


ST  ATUS= 1 ,  MOVE 
AGAIN  BEFORE  DWELL 
STATUS=2 ,  DWELL 
BEFORE  MOVE  AGAIN 


; ********************* ************************************** 
SEMVE  ASSIGN, XX ( 16) =XX ( 15) +XX (25) +1 ,  1;  COMPUTE  NEW  GRID 

NUMBER 

ACT; 

EVENT, 6,1;  GET  STATUS  OF  NEW 

GRID 

ACT, ,XX (24) .EQ.O,TOPRW;  STATUS=0,  DO  NOT 

MOVE 


ACT; 

GOON, 2; 

ACT,,, HFMVE;  MOVE  TARGET  OUT  OF 

OLD  GRID 

ACT  ,  STOF'A  (5)  ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON, 1; 

ACT, , XX (26) . EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 


ACT; 

ASSIGN, XX (15)=XX ( 16) , 

XX ( 17) =XX ( 17) +1 , 

XX ( 18) =XX ( 18) +1 ,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 


ACT; 

EVENT, 9, 1; 

ACT; 

GOON, 2; 

ACT, , , HFMVE; 

ACT , STOPA (5) ; 

GOON, 1 ; 

ACT, , XX (26) . EQ. 1 , WTTGT ; 

ACT; 

GOON, 1 ; 

ACT,  ,XX  (24)  .EQ.  l,TOF'RW; 
ACT, , , MOVE ; 


PRINT  MESSAGE  ABOUT 
TGT  MOVE 


MOVE  TARGET  TO  THE 
CENTER 

OF  THE  NEW  GRID 

STOP,  TARGET 
DETECTED 


STATUS=1 ,  MOVE 
AGAIN  BEFORE  DWELL 
STATUS=2 ,  DWELL 
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BEFORE  HOVE  AGAIN 

; *********************************************************** 
SHOVE  ASSIGN, XX(16)=XX(15)+XX (25) ,  1;  CONFUTE  NEW  GRID 

NUHBER 

ACT; 

EVENT ,6,1;  GET  STATUS  OF  NEW 

GRID 

ACT, ,XX (24) .EQ.O,TOFRW;  STATUS=0,  DO  NOT 

HOVE 

ACT; 

GOON, 2; 

ACT , , , HFHVE;  HOVE  TARGET  OUT  OF 

OLD  GRID 

ACT  ,  STOF'A  (5)  ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON , 1 ; 

ACT, ,XX(26> .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX(15)=XX  (16)  , 

XX(17)=XX(17)+1,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT ,9,1;  PRINT  HESS AGE  ABOUT 

TGT  HOVE 

ACT; 

GOON, 2; 

ACT,,, HFHVE;  HOVE  TARGET  TO  THE 

CENTER 

ACT  ,  STOF'A  (5)  ;  OF  THE  NEW  GRID 

GOON, 1; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

GOON, 1 ; 

ACT, ,XX (24) .EQ. 1 ,TOPRW;  STATUS=1 ,  HOVE 

AGAIN  BEFORE  DWELL 

ACT,,, HOVE;  STATUS=2 ,  DWELL 

BEFORE  HOVE  AGAIN 

; ******  *************** ************** *************** ********* 
SWHVE  ASSIGN, XX ( 16) =XX ( 15) +XX (25) -1 ,  1;  COHPUTE  NEW  GRID 

NUHBER 

ACT; 

EVENT, 6,1;  GET  STATUS  OF  NEW 

GRID 

ACT, ,XX (24) .EQ.O,TOPRW;  STATUS=0,  DO  NOT 

HOVE 

ACT; 

GOON, 2; 

ACT,,, HFHVE;  HOVE  TARGET  OUT  OF 

OLD  GRID 

TO  THE  EDGE  OF  THE 
NEW  GRID 


ACT ,STQPA (5) 


GOON , 1 ; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX ( 15) =XX ( 16) , 

XX (17) =XX ( 17) +1 , 

XX ( 18) =XX ( 18) -1 ,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT ,9,1;  PRINT  MESSAGE  ABOUT 

TGT  MOVE 

ACT; 

GOON, 2; 

ACT , , , HFMVE ;  MOVE  TARGET  TO  THE 

CENTER 

ACT , STOPA (5 ) ;  OF  THE  NEW  GRID 

GOON , 1 ; 

ACT,  ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

GOON, 1 ; 

ACT, , XX (24) . EQ. 1 ,TOPRW;  STATUS=1,  MOVE 

AGAIN  BEFORE  DWELL 

ACT,,, MOVE;  STATUS=2 ,  DWELL 

BEFORE  MOVE  AGAIN 

; *********************************************************** 

WMOVE  ASSIGN, XX ( 16) =XX < 15) -1 ,  1;  COMPUTE  NEW  GRID 

NUMBER 

ACT; 

EVENT, 6,1;  GET  STATUS  OF  NEW 

GRID 

ACT , , X X ( 24) . EQ. 0 , TOPRW;  STATUS=0,  DO  NOT 

MOVE 

ACT; 

GOON, 2; 

ACT,,, HFMVE;  MOVE  TARGET  OUT  OF 

OLD  GRID 

ACT, STOPA (5) ;  TO  THE  EDGE  OF  THE 

NEW  GRID 

GOON , 1 ; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN, XX (15) =X  X ( 16)  , 

XX ( 18) =XX < 18) -1 ,  1;  CHANGE  TGT  LOCATION 

TO  NEW  GRID 

ACT; 

EVENT ,9,1;  PRINT  MESSAGE  ABOUT 

TGT  MOVE 

ACT; 

GOON, 2; 

ACT,,, HFMVE;  MOVE  TARGET  TO  THE 

CENTER 


OF  THE  NEW  GRID 


ACT , STQPA (5) ; 

GOON , 1 ; 

ACT, ,XX (26) .EQ. 1 ,WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

GOON, 1; 

ACT, ,XX (24) .EQ. 1 ,TOPRW;  STATUS=1,  MOVE 

AGAIN  BEFORE  DWELL 

ACT,,, MOVE;  STATUS=2,  DWELL 

BEFORE  MOVE  AGAIN 

• ****  DELAY  FOR  DWELL  TIME  ******************************* 
DLAGN  GOON , 1 ; 

ACT, XX (21);  DELAY  FOR  DWELL 

TIME 

GOON , 1 ; 

ACT, , XX (26) . EQ. 1 , WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASSIGN,STOPA=l ,  1;  RELEASE  TGT  ENTITY 

ACT; 

TERM; 

■ ****  FUEL  DECREASE  ************************************** 
FULDN  GOON , 1 ; 

ACT, XX (36);  FLY  ONE  MORE 

MISSION  SEGMENT 

ASSIGN, XX (6)=XX (6)-XX (37) ,  1;  DECREASE  FUEL  USED 

LAST  SEGMENT 

ACT, , XX (26) .EQ. 1, WTTGT;  STOP,  TARGET 

DETECTED 

ACT  ,  , , FULDN; 

;*♦**  MOVE  THE  TARGET  ONE  HALF  OF  A  GRID  **************** 
HFMVE  GOON , 1 ; 

ACT, XX (20);  DELAY  FOR  ONE  HALF 

OF  A  MOVE 

GOON, 1 ; 

ACT, ,XX (26) .EQ. 1 , WTTGT;  STOP,  TARGET 

DETECTED 

ACT; 

ASS  I GN ,  STQF'A=5 ,  1;  RELEASE  TARGET 

ENTITY 

ACT; 

TERM; 

;♦***  MISSION  COMPLETION  MESSAGE  ************************* 
E1MSG  ASSIGN, XX (46) =1. ,  1; 

ACT; 

EVENT, 10,1;  PRINT  MESSAGE 

ABOUT  NO  SEARCH 

;  TIME  ALLOWED 

ACT , , , ENFLT ; 

E2MSG  ASSIGN, XX (46) =2. ,  1; 

ACT; 

EVENT, 10,1; 


PRINT  MESSAGE 
ABOUT  ALL  SEARCH 


TIME  USED 


ACT, , , ENFLT ; 

E3MSG  ASSIGN, XX (46) =3. ,  1; 

ACT; 

EVENT, 10, 1 ; 


ACT, , , ENFLT; 

E4MSG  ASSIGN, XX (46) =4. ,  1; 

ACT; 

EVENT, 10,1; 

ACT, , , IDENT; 

;***♦  END  FLIGHT, NO  DETECTION 
ENFLT  ASSIGN, XX ( 1 1) =XX ( 1 1 ) +1 , 
XX (12)=XX ( 12) +1 , 
XX (26) =1 , 

STQPA= 1 , 

ST0PA=5 ,  1 ; 


PRINT  MESSAGE 
ABOUT  ACFT  BELOW 
MINIMUM  FUEL 


PRINT  MESSAGE 
ABOUT  TGT  DETECTED 

***************** *********** 


COUNT  AIRCRAFT 
MISSES  AND 


MARK  TARGET  FOR 
DESTRUCTION 


ACT; 

COLCT, XX (12) , NUMBER  MISSES; 
ACT, , , WTTGT ; 


;#***  END  FLIGHT,  FOUND  TARGET 

* 

IDENT  ASSIGN, XX ( 1 1 ) =XX ( 1 1 ) +1 , 

XX (28) =XX ( 28 ) + 1 , 

XX (26) =1 , 

ST0PA=1 , 

ST0PA=5,  1; 


ACT; 

COLCT, XX (28) , NUMBER  HITS; 

ACT, , , WTTGT; 

;*«**  AIRCRAFT  WAITS  FOR  ALL  OTHER 
WTTGT  ACCUM ,4,4, ,1; 


ACT; 

ASSIGN, XX (29) =XX (28) /XX  < 1 1 )  , 
ACT; 

EVENT, 12,1; 


[■************************ 


COUNT  AIRCRAFT 
HITS  AND 
MARK  TARGET  FOR 
DESTRUCTION 


STRAY  ENTITIES  ******** 
WAIT  FOR  ALL 
ENTITIES (AIRCRAFT, 
TARGET, FUEL,  AND 
DWELL/MOVE) 

1;  COMPUTE  FRACTION  OF 
TGT'S  FOUND 

PRINT  THE  GRIDS  THE 
AIRCRAFT  SEARCHED 
AND  THE  TIME  THE 
SEARCH  BEGAN 


CQLCT  ,  X  X  ( 29 )  ,  Pr  DETECT  I  ON ; 

ACT, ,XX(11) .EQ.XX(14) , DONE;  RUN  COMPLETE,  FLEW 

MAX  NUMBER 

;  OF  MISSIONS  ALLOWED 

ACT , , , ACFT ;  FLY  NEXT  MISSION 

;  *******************************************  **************** 
DONE  TERM , 1 ;  END  RUN 

ENDNETWORK; 

;  *********************************************************** 
I N I T , 0 . ,1 OOOOOO . , NO / 1 ; 

SEEDS ,0(1) / Y , 0 ( 2 ) / Y , 0 ( 3  > / Y , O ( 4 ) / Y , 0 ( 5 ) / Y , 0  <  6 ) / Y , 0 ( 7 ) / Y , 

0(8) / Y , 0 ( 9 ) /Y; 

SIMULATE; 

FIN; 


Barr,  Donald  R.  and  Peter  W.  Zehna.  Probability 
Model  i ng  Uncertainty.  Menlo  Park,  Cali-fornia:  Addison- 
Wesley  Publishing  Company,  1983. 

Bhat ,  U.  Narayan.  El  ements  of  App  1 1  ed  Stochasti  c 
Processes .  New  York:  John  Wiley  and  Sons,  Inc.,  1972. 

Brown,  S.  S.  "Optimal  Search  -for  a  Moving  Target  in 
Discrete  Time  and  Space,"  Operations  Research ,  28: 

1275-1289  (Nov-Dee  1980). 

Hillier,  Frederick  S.  and  Gerald  J.  Lieberman. 
Introduction  to  Operations  Research  (Third  Edition). 
Oakland,  Cali-fornia:  Holden-Day,  Inc.,  1980. 

Meagher ,  James  J.,  Project  Manager  -for  SRAM  II  Integra¬ 
tion.  Personal  interview,  Wr l ght-Patter son  AFB,  Ohio, 
2  September  1985. 

Mendenhall,  William  and  others.  Mathematical  Statistics 
Wi th  App 1 ications  (Second  Edition).  Boston, 
Massachusetts:  Duxbury  Press,  1981. 

Montgomery,  Douglas  C.  Desi gn  and  Anal ysi s  of 
Experiments  (Second  Edition).  New  York:  John  Wiley  and 
Sons,  1984. 

Neter,  John  and  others.  Appl l ed  Linear  Statistical 
Model s  (Second  Edition).  Homewood,  Illinois:  Richard 
D.  Irwin,  Inc.,  1985. 

Pritsker,  A.  Alan  B.  Introducti on  to  Si mul at l on  and 
SLAM  1 1  (Second  Edition).  West  La-fayette,  Indiana: 
Systems  Publishing  Corporation,  1984. 

Classi-fied  Document:  Qualified  requestors  may  obtain 

this  reference  from  AFIT/ENS,  Wr i ght-Patter son  AFB, 
Ohio. 

Classified  Document:  Qualified  requestors  may  obtain 

this  reference  from  AFIT/ENS,  Wr  i  ght-F’atterson  AFB, 
Ohio. 

Stone,  L.  D.  Theor y  of  Opt i mal  Search .  New  York: 

Academic  Press,  1975. 

Sutton,  Jerry  P. ,  Deputy  Director  of  Projects  for  the 
B-l.  Personal  interviews,  Wr i ght-Patterson  AFB,  Ohio, 
28  March  through  13  December  1985. 


Trees,  Ronald  J.,  Bomber  Division  Directorate  -for 
Strategic  Force  Analyses  Assistant  Chief  o-f  Staff, 
Studies  and  Analyses.  "Modeling  Strategic  Relocatable 
Targets  Missions."  Brie-fing  presented  to  MORS  A-l 
Working  Group.  HO  USAF/SSAB,  F'entagon,  5  June  1984. 

Tierney,  Luke  and  Joseph  B.  Kadane.  "Surveillance 
Search  tor  a  Moving  Target,"  Qper at i ons  Resear ch ,  5 1 ; 
720-738  (July-August  1983). 

Washburn,  Alan  R.  Search  and  Detect i on .  Monterey, 
California:  Military  Applications  Section  Operations 
Research  Society  of  America,  May  1981. 


Captain  Donald  B.  Qlynick  was  born  on  15  April  1954  at 


Eglin  AFB ,  Florida.  He  graduated  -from  high  school  in  Can¬ 
ton,  Ohio  in  1972  and  attended  the  United  States  Air  Force 
Academy  -from  which  he  received  a  Bachelor  o-f  Science  degree 
in  Computer  Science  in  June  1976.  Upon  receiving  his  com¬ 
mission  in  the  USAF  -following  graduation  -from  the  Academy, 
he  started  undergraduate  pilot  training  at  Craig  AFB,  Ala¬ 
bama.  Pilot  training  was  completed  at  Williams  AFB,  Arizona 
due  to  the  closure  o-f  Craig  AFB.  Graduation  -from  pilot 
training  occurred  in  October  1977  where  he  then  served  as  a 
copilot  and  pilot  -for  the  23rd  Bomber  Maintenance  Squadron 
at  Minot  AFB,  North  Dakota  until  January  1983.  Following  a 
tour  as  a  command  post  Emergency  Actions  Controller  -for  the 
5th  Bomb  Wing  at  Minot  AFB  -from  January  1983  until  April 
1984,  he  entered  the  School  o-f  Engineering,  Air  Force  Insti¬ 
tute  o-f  Technology  in  May  1984. 


Permanent  address:  339  Saratoga  Avenue  N.  W. 

Canton,  Ohio  44708 


1 36 


UNCLASSIFIED 


security  classification  of  this  page 


1*  REPORT  SECURITY  CLASSIFICATION 


2».  security  classification  authority 


7b.  oeclassification/downgraoing  schedule 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AFIT/GOR/ENS/8  5D- 1 5 


6a.  NAME  OF  PERFORMING  ORGANIZATION  |6b.  OFFICE  SYMBOL 

V  Sch°o1  of  Engineering  | 


REPORT  DOCUMENTATION  PAGE 


1b.  RESTRICTIVE  MARKINGS 


6c.  ADORESS  (City,  State  and  ZIP  Code) 

Air  Force  Institute  of  Technology 
Wright-Pat ter son  AFB,  Ohio  45433 


l.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

B-1  Project  Office 


8b.  OFFICE  SYMBOL 
(If  applicable) 

ASD/B-1M 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  N UMBER (S) 


7*.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  State  and  ZIP  Code) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Be.  ADDRESS  (City,  State  and  ZIP  Code) 

10.  SOURCE  OF  FUNDING  NOS. 

Wright-Patterson  AFB,  Ohio  45433 

PROGRAM 
ELEMENT  NO.  1 

PROJECT 

NO. 

TASK 

NO. 

WORK  UNIT 
NO. 

11.  TITLE  (Include  Security  Classification) 

See  Box  19 

12.  PERSONAL  AUTHOR(S) 

Donald  B.  Olynick,  B.S.,  Capt,  USAF 


13a  type  OF  REPORT 

MS  Thesis 


IB.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM  _  TO 


14.  DATE  OF  REPORT  (Yr  ,  Mo..  Day) 

1985  December 


15.  PAGE  COUNT 


COSATI  CODES 


18.  SUBJECT  TERMS  f Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Strategic  Bombing,  Markov  Processes,  Moving  Targets, 
Simulation,  Mathematical  Models 


10.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


FIELD 

GROUP 

09 

02 

Title:  HOW  THE  TIMING  OF  UPDATES  ON  THE  LOCATION  OF  STRATEGIC 
RELOCATABLE  TARGETS  CHANGES  THE  PROBABILITY  OF  DETECTION 

Thesis  Advisor:  William  F.  Rowell,  Major,  USAF 

Assistant  Professor,  Department  of  Operational  Science 


.  u  -iv  him*  M  M  mj. 


N>  Tom 


20.  OISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIFIEO/UN LIMITE  O  X  SAME  AS  RPT.  □  OTIC  USERS  □ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

William  F.  Rowell,  Major,  USAF 


DD  FORM  1473,  83  APR  EOITION  OF  1  JAN 


t  et  TcokAoJorr 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TE  LEPHONE  NUMBER  22c  OFFICE  SYMBOL 

(Include  Area  Code i 

513-255-3362  AFIT/ENS 


EOITION  OF  1  JAN  73  IS  OBSOLETE.  ITNOT  A^TFTF.n 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


SECURITY  classification  of  this  page 


This  project  simulates  a  bomber  aircraft  mission  attacking  a  stra¬ 
tegic  relocatable  target  to  develop  the  relationship  between  the  timing 
of  intelligence  updates  and  the  probability  of  detection.  After  one 
mission  is  simulated,  the  result  is  recorded  as  a  hit  or  a  miss.  For 
one  run  using  one  update  time,  the  percentage  of  hits  from  the  total 
number  simulated  beoomes  the  estimate  for  the  probability  of  detection. 

Target  movement  is  modeled  as  Markov  chain  with  aircraft  missions 
simulated  independently.  The  simulation  is  a  combined  network  discrete- 
event  orientation  using  SLAM  and  FORTRAN  as  the  computer  languages.  Out¬ 
put  routines  track  the  target  movement  through  a  specified  grid  area  as 
well  as  tracking  the  aircraft  search  patterns. 

Results  are  displayed  graphically  by  plotting  the  estimated  detec¬ 
tion  probability  against  the  corresponding  update  time  used  to  arrive  at 
the  estimate.  Procedures  for  developing  confidence  bands  around  the  re¬ 
sulting  plot  are  also  discussed.  More  research  is  recommended  in  two  areas: 
the  search  pattern  used  by  the  aircraft  and  procedures  for  estimating  how 
long  the  aircraft  has  been  in  the  observed  position  *dien  an  intelligence 
update  is  taken.  The  conclusion  is  that  the  time  remaining  before  the  air¬ 
craft  arrives  in  the  search  area  when  the  update  is  received  should  be  less 
than  half  the  total  dwell  time  of  the  target  in  order  to  realize  sig¬ 
nificant  improvements  in  the  detection  probabilities. 
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